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Chapter 1 

Introduction 



Following part of the work iniciated by A. Agrachev and A. Sarychev in 8 
we study controllability, by means of low modes forcing, of incompressible 2D 
Navier-Stokes (NS) Equations on the two dimensional rectangle with tangent 
velocity on the boundary. 

In the present paper we deal with the 2D NS system 

u t + (u- V)u + Vp = vAu + F(xi,x 2 ) + v(t,xi,x 2 ) (1.0.1) 

(1.0.2) 
(1.0.3) 
(1.0.4) 



— I dx 2 



V 


■ u = 





on 


R 


u • 


n = 





on 


dR 


v x 


■ u = 





on 


OR 



Where R :— {(xi,x 2 ) G K 2 | a± < x\ < a 2 ; b\ < x 2 < b 2 } and V 

and n is the unit normal to the boundary. In the equation l|1.0.1|) u is the velocity 
of the fluid "particle"; p is the pressure; the only nonlinear term of the equation 
— (u ■ V)it — is called the inertial term; vAu is called the viscosity term, v > 
is the coeficient of viscosity; F is an external force and; v is a control at our 
disposal. We are interested in the case where v is a degenerate forcing, i.e., v is 
a finite sum of harmonics — v = J^keK 1 v k(t)Ek, where are eigenf unctions 
of the Stokes operator. So the components ffc(t), k G JC 1 , t € [0, T] are our 
controls which are measurable essentially bounded functions. We shall study 
Galerkin approximations, say, big enough to contain a set of modes we want to 
observe. 

A natural way to study the NSE is to study its evolution on subspaces of 
Sobolev spaces; such subspaces depend on the boundary conditions. 

We shall note by L 2 (R) the space of Lebesgue measurable square integrable 
real functions defined on R and by L 2 (i?) the product space L 2 (R) 2 . Similarly 
H 1 := {/ S L 2 {R) | §L G L 2 (R), j = 1,2} and, H*(R) := H^R) 2 . 

The most studied boundary conditions are the full Dirichlet (u = on the 
boundary, also called no-slip) and periodic (tti and u 2 are periodic). For the 
full Dirichlet conditions the study is done in the spaces 

H D := {u G L 2 (R) | V • u = & u-n = 0ondR} 
V D := {u G H\R) | V ■ u = & u = 0ondR} 
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CHAPTER 1. INTRODUCTION 



and for the periodic boundary conditions the study is done in 

Hp := {u e L 2 er (i?) | V • u = 0} 
V P := {u e Hp er (i?) | V • u = 0}. 

In the periodic case sometimes the special case of zero average space — J R udx = 
— is considered, in this case we have to take this new condition into account in 
the definition of the spaces Hp and Vp . For more information on the boundary 
conditions above see for example ^5] or |14j . 

For the boundary conditions 1)1.0.3(1 . 1(1. 0.4(1 the spaces 

H := {u e L 2 (i?) | V • u = & u-n = 0ondi?}; 
V := closure of V x (R) on H x (ii); (1.0.5) 
D(A) := {u G H 2 (i?) | V ■ u = & (u ■ n = A V" 1 ■ it = 0) ondi?}; 

where V X {R) := {u € C*°°(i?) | V-u = & (u-n = 0AV ± -u = 0)onOi?}, 
are those where we shall consider the evolution of the NSE on. 

In Chapter|21we prove the existence of weak and strong solutions for 1(1.0.1(1 - 
((l.U.4(l , as well as its uniqueness and continuous dependence in the initial data. 
For the case of weak solutions we proceed as in |T5] for the case of no-slip 
boundary conditions. To obtain strong solutions we just have to ask for some 
regularity on the initial data. 

In Chapter |31 for strong solutions, we prove controllability, by means of low 
modes forcing, of Galerkin approximations of the infinite-dimensional system 
associated with (11.0.1(1 - 1(1.0.41) . 

In Chapter^lwe prove the so called controllability in observed component, again 
by means of low modes forcing. In other words we prove that the projection, 
onto any finite dimensional subspace (the space spanned by the modes we want 
to observe), of the attainable set from any point is surjective. 

Finally, we end with Chapter [3] and with the prove of L 2 - Approximate Con- 
trollability the is a straightforward corollary of some tools we have presented in 
Chapter E| 

The author is gratefull to A. Agrachev and A. Sarychev for the inspiring and 
hclpfull discussions on the subject and, for the sugestions in the improvement 
of the text. 

The author would like to thank FCT (Portuguese Foundation for Science 
and Tecnology) for financial support and, SISSA-ISAS (International School for 
Advanced Studies) for hospitality. 



Chapter 2 

Existence, Continuity and 
Uniqueness. 



2.1 The Spaces. 

2.1.1 Recollection of Auxiliary Material on the Spaces 
L 2 (i?) and H 1 ^). 

Recall that since L 2 (R) is a Hilbert space for the scalar product 

(u, v)i :— / uvdx, (2.1.1) 
Jr 

then L 2 (i?) is a Hilbert space for the product topology and the scalar product 
is 

(u, v):= u -vdx. (2-1.2) 
Jr 

We note that 

2 

(u, v) = / u-vdx = / UiVi dx = (ui, t>i)i + (u 2 , v 2 )i- 
Jr l=1 Jr 

The norms associated with the previous scalar products shall be represented by 

\u\\ := (u, u)i (2.1.3) 
M 2 := (u, u) (2.1.4) 



We note that 



M 2 :H Ul | 2 + M 2 



so, the norm | • | is the product norm | - | x x I * 1 1 - Similarly the Sobolev space 
H 1 (R) is a Hilbert space for the scalar product 

((«, w)) i := (u, w)i + J2(§^, = («, v)i + (Vu, Vv). (2.1.5) 
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then H 1 (i?) is a Hilbert space for the product topology and the scalar product 
is 

2 2 

((u, w)) := «j))oi = («, w) + ^(V^, V«i) (2.1.6) 

j=l i=l 

The norms associated with the previous scalar products shall be represented by 

\\u\\ 2 01 := ((«,«))oi (2.1.7) 

Nlo : = ((«»«))o (2-1.8) 

We note that 

IMIo : = IMIoi + IMIoi 
so, the norm |j • || is the product norm || • || i x || • ||oi- 

2.1.2 Some Properties of the Spaces V and H. 

Lemma 2.1.1. H coincides with the closure ofD\(R) in L 2 (i?). 

Proof. It is well known from the study of the NSE with full Dirichlet boundary 
conditions that H is the closure of T>(R) in L 2 (R), where V{R) is the set of 
solenoidal (or divergence free) smooth functions with support in R. Since 

V{R) C 2>i(J2) C H 

we conclude that H coincides with the closure of T>\ (R) in L 2 (R) . □ 

We use the same argument, since T>i C V C H to conclude that 

Corollary 2.1.2. H is the closure of V in L 2 . 1 

In the study of NSE some classical imbedding and compactness theorems 
are frequently used, we start by presenting some of them which we shall need 
(for our particular equation). 

Proposition 2.1.3. For all u <G H 1 (R) and for all q, 1 < q < oo we have 

\\u\\li(r) < Ci(?)|M|oi- 
Corollary 2.1 A. For all u e H 1 (i?) and for all q, 1 < q < oo we have 

\\uh"(R) < Ci(g)||u|| . 
Proposition 2.1.5. The imbedding 

H\R) -» L 2 (R) 
u u 

is compact. 

Corollary 2.1.6. The imbedding 

U\R) -» \?{R) 

u u 

is compact. 

1 In other words V is dense in the closed space H 
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Corollary 2.1.7. The inclusion 



i:V -> H (2.1.9) 

u i — ► 11 



is continuous and compact. 



Proof. By Corollarv l2.1.4l inclusion ((2. 1.9(1 is continuous and, by Corollary 12. 1.61 
it is compact. □ 

2.1.3 Poincare Inequality. Equivalent Norms. 

We have the following Lemma: 2 
Lemma 2.1.8. 

• For any seminorm p in H^i?) satisfying "(p(a) = A a e M 2 ) => a = " 
we have 

\u\ < c(||Vu|| L 2 (fl)2 +p(tt)), Vw S H x (i?) 

and, 

• T/ie seminorm 

p(u) := J |u-n|dT 

satisfies the required condition. 
Remark 1. The product space L 2 (i?) 2 is a Hilbert space for the scalar product 



A consequence of this lemma is a Poincare-like inequality: 

Vti S V[ \u\ < c|| Vw|| L 2 (fl)2 ]. (2.1.10) 
Corollary 2.1.9. The norms \\u\\ := ||Vm||l2(^)2 and \\ ■ ||o are equivalent in 

v 3 ' 



2.2 The Duals of V and H. 

From now we shall consider the space V endowed with the norm || • ||. From 
Corollaries l2 . 1 . 21 12 . 1 . 71 and !2 . 1 . 91 we obtain that the inclusion l|2 . 1 .9(1 (^considering 
V endowed with || • ||) is dense continuous and compact. 

Since H is an Hilbert Space due to the Riesz Representation Theorem we 
can identify H with its dual H' . In this way we arrive to the inclusions 

V dH = H' dV' (2.2.1) 



2 See El subsections II. 1.4 and III.2.2 
3 Note that || ■ ||o = I ■ I + II ■ II- 
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where both inclusions are linear, dense and continuous. Indeed we already know 
that the first inclusion has these properties. In the second inclusion we consider 
the map i', dual to the map i defined in l|2.1.9[) . 

i':H'^ V 

< i(v), h> =<v, i'(h) >, 

that is linear and continuous — we consider V' endowed with the classical norm 
ll/H := sup||„|| =1 | < /, v > |. The density of H' in V is a corollary of the 
following Lemma, which one can find in 0], section (II. 6), 

Lemma 2.2.1. Let G : D{G) C E — > F be a (nonbounded) closed operator 
with D(G)=E. Then 

(i) N(G) = RiG') 1 - 
(ii) N(G') = RIG) 1 - 



(Hi) N(G) ± D R(G') 
(iv) N(G') ± = R(G). 

Where X means closure of X; N(G) stays for the Kernel of G and R(G) 
for the image of G. Applying this Lemma to our inclusions the injectivity of 
the first inclusion implies the density of the second — (i). Moreover we also 
see that, since the first inclusion is dense, it follows the injectivity of the second 

(«)• 

Remark 2. The domain D(i) of i is V. By the continuity of i we conclude its 
closedness. Recall that a operator G : E — > F is said closed if for any sequence 

t \ ■ nfn\ u +u A Un ^ u inE h \ u E D ( G ^ 
(u n ) in D(G) such that < we have < 

Au n -> / inF \Au=f 



2.3 Fourier Series. A Basis in T>i{R) for H and 
V. 

We start by noting that if in the case of the general rectangle R :— [a\, 02] x 
[bi , 62] o-i < 0-2 , b\ < 62 we make the change of variables 

Z\ := xi - oi, z 2 := x 2 - bi, 

since the differential operators V, A are invariant under tranlations, the Navier- 
Stokes Equation does not change in the variables [z%, z-i). To simplify the 
exposition from now we shall deal with the rectanlge 

R := {(xi,x 2 ) G M 2 I < xi < a; < x 2 < b}. 

First we note that under the condition u ■ n = on dR, the condition V ■ u = 
on dR equivals |^ = — on dR, because on the left and right faces we 
have = and, on the top and botton faces |^ = 0. It is well known that 
{sin( I1 J i ) I n £ No} is a basis for the functions in L 2 ([0, L]) vanishing on 
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and L and, {cos( B J £ ) | n e N} is a basis for the functions in L 2 ([0, L}) with 
vanishing velocity on the boundary points and L. 4 So 

r . / nnx\ ( rrnrx\ "i 
jsinl J cos( — — J | ne N , me Nj 

is a basis for L 2 (R) vanishing on the right and left faces of the rectangle and, 

r / nirx\ . ( miTx\ "i 
|cos ( J sin I — — J | n e N, m e N j 

is a basis for L 2 (R) vanishing on the botton and top faces of the rectangle. 

Now consider a function u £ H, this function, being an element of L 2 (i?), 
can be written as 



u ■■ 



'«1 



'£ ne N„Mi fc sin(^) CO s(^)' 



m6M 



u 2 J \ EmeNc u 2k cos(^) sin(^) ' 



u 



Since u is solenoidal — V • u = 0, we conclude that 

*A = /E feeN g^sin(^)cos(^) N 
vW i,E feeN ^2 fe cos(^)sin(^) / 

and that for each k e Nq we have 

fcl7T fc 2 7T 

wife 1- u 2 k—r- = 0. 

a o 

This means that = =: u k and, we arrive to 



^2 u kW k , W k 



*0 

We put 



=M2L S i n (^£i) C Os(^f^ ) 

*^ CO s(*^)sin(^f^) y ' 



W := {W k | fee Nq}. (2.3.1) 



2.3.1 Fourier Characterization of H and V. 

We start by computing the scalar product (Wk, W z ) between two elements of 
W: 

(W k ,W z )= [ W k -W z dx 
Jr 

f k 2 Z 2 1T 2 . /fcl7TXl N ,k 2 lTX 2 , . fZlTTXl-, (Z 2 -KX 2 -. , 

- J R -p- Sln (— ) cos (— ) Sln (— ) cos (— ) d - 

/ fcizi7r 2 /fci7ra;i, . ,k 2 -KX 2 , ,Z\-KX\. . ,z 2 -kx 2 , 

+ / - — cos( sm — - — ) cos( sin — - — ) ax 

J R a z a ' b ' a ' b ' 



>R 

'0 iik^z 

z 



k^ iik = z. 5 



4 Hcre N is the set of natural numbers and No := N \ {0} 
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where k := — n 2 + |§-^ . Hence the square of the L 2 (i?) norm of each element 

in the orthogonal family W equals —k 2 ^-. 

Given u, v € H the scalar product between them is then 

where 

fcGNg fceNg 

The norm of u in _ff results 

feGNg 

and u € H means that 

— fcu 2 , < +00. 

If u, v € y we want to compute the scalar product 

((u, «)) = ((V«i, Vwi)) + ((V« 2) V«a)). 

First we note that 



Vm 2 



' E k enl cos(^) cos(^) N 

,E fceNi « fc (-*P=^) sin(^) sin(^) / 

I J2 k eK M^rT) cos(^) cos(^) 



and we obtain analogous expressions for V«i and V«2- Hence, doing the com- 
putations, we arrive to 

((«> w )) = X! k 2 ^u k v k . 



The norm of it in results 



i 2 = ^£*v fc ( 2 - 3 - 3 ) 



and u G V means that 



fc 2 u^ < +00. 

ken* 



From caracterizations ()2.3.2J) and (|2.3.3JI we can see easily, as referred in Corol- 
lary ^. 1.91 that the norms || • ||o and || • || are equivalent in V. Indeed, — k goes 

5 This follows from the fact that the family of sines is orthogonal in L 2 ([0, 7r]), as well is 
the family of cosines. 
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to +00 when max{fci, fe} does 6 so, there is N £ No and a constant C > 1 such 
that if maxjfci, fe} > N we have — k > 1 and, 

-k<c ~ k2 - ( 2 - 3 - 4 ) 

max \k\ , k-2 }< AT max {fci , &2 }< N 

Therefore \u\ < C||w|| and then || • ||o and || • || are equivalent. 

When u £ V we can not guarantee, in general, neither Au £ H nor Am £ L 2 . 
Anyway if we compute (Au, v) for u, v £ V we obtain 

(Au, v) = Y k 2 u 2 k < +00 = -{{u, v)). (2.3.5) 

fcGNg 

Hence the result of (Aw, i>) is finite and is properly defined even if Au ^ L 2 and 
so the operation (Au, v) seems to have no sense in itself. Below we will give a 
sense to it. 



2.4 The Operators A and B. 
2.4.1 The Operator A. 

First we note that a consequence of the identifications H2.2.1JI , the scalar product 
in H of / £ H and u £ V is the same as the scalar product of / and u in the 
duality between V and V 

V/ £HVueV{< /,u >=(/,«)]. (2.4.1) 

For each u G V, the form 

-> R (2.4.2) 

V H- > ((ti, U)) 

is linear and continuous on V. Therefore there is an element of V we shall 
denote by Au satisfying 

< Au, v >= ((u, v)), Vu £ V (2.4.3) 

and A is clearly linear and continuous; we can also see that A is an isomorphism 
between V and V . Indeed the form a(u, v) :— ((u, v)) is bilinear, continuous 
and coersive on V x V so, by Lax-Milgram Theorem we can conclude that A is 
an isomorphism between V and V . 
By equation (|2.3.5J) we see that 

< Au, v >= (-Au, v) Vu, v £ V. (2.4.4) 



6 For M := max{fci, k2\ and m := rnax{a 2 , b 2 }, we have — k = w 2 (^- + > t 2 -^- that 
goes to 00 when M does. 
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Further Properties of A. 

As we have said above for v e V, Av can be not in L 2 . Now we define a subset 
D{A) of V, we shall call the domain of A, such that 

A : D(A) -» H. 

The domain of A is defined by 

D(A) := {v E H 2 | u e V & V x ■ v = OonT} 

= {weH 2 | V • v = & (v ■ n = A V 1 " • u = OonT)} 
= closure ofI?i(i?) inH 2 . 

Where T := <9i? and H 2 is the Sobolev space H 2 = (W 2 - 2 ) 2 where 

W 2 ' 2 := {u e L 2 \ ^eL 2 , VH<2} 
which is a Hilbert space for the scalar product 

^ v ) 2 ■■= E (-^S-. = 2. (fl£T> a^)- 

ie{l,2} [q|<2 
l«l<2 

The norm corresponding to this scalar product is defined by 



M 2 + IMI 2 + Mf 2] , 



with 



'«" a ™- = E 2= E (ifl?«i a + i#*i a + i^iVi 2 ) 7 



ip]-" 1-^5 

|a|=2 |a|=2 

We can see that the norm | • 1 2 is equivalent to that | • I2+ coming from the 
scalar product defined by 

(it, v) 2+ = (u, v) 2 + (df 2 u, di >2 v). 

Indeed \u\ 2 < \u\ 2 + < 2\u\ 2 . 

Now we look for a Fourier caracterization of the elements on D(A): For that 
we need to write down the expressions for the second order derivatives of it: 





d\u = 

7 d? stays for ^ and, d\ 2 stays for g jf dx2 ■ 



f-^j sin(^i; 

cos (*^) sin(^2 
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Computing (u, v) 2 + for u, v £ D(A) we obtain 



ab v—v — ab 



(u, v) 2 + = — £ ~ ku kVk + -j £ k 
ab 



UkVk 



A fc 1 ■ h 



b 2 a ) \a 3 J U 3 / 



UkVk 



" £ (~fc + P - k 3 )u k v k . 



J2(-k + k 2 -P)ul 



4 

Hence for an element u of D(A); 

I |2 « 6 

M2+ = j 

Using the constants C and TV of the equation (|2.3.4(l and choosing C\ such that 
£ fc 2 <Ci £ -fc 3 , (2.4.5) 

max{fci , /C2 max{fci , /C2 }<-/V 

we conclude that 

£ -P < £ (-jfe + P - P) < ((C + l)Ci + 1) £ -P. 

We define 

[u, u)[2] := (w, u)[2]_ + (9f 2 u, 9 2 j2 w). 

Then all the norms |tt|a+; \ u \"2\ M[2l- an d |xt| [ 2 ] are equivalent. Moreover for 
the last one we have the nice representations: 

(«> w )[2] = 7^ ~k 3 u k v k ; 



u 



- - ab ^ -put 



4 

Now we compute 

V(V-it) +V X (V" L -it), u 
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so 

Au = V(V-u) + V x (V x -u) (2.4.6) 

and, for u solenoidal we have Am = V ± (V 1 * • u). Moreover for each u £ D(A) 
we have 

V- Au = V-V-^V 1 - -it) = 0; 
n • Au = n • V X (V X • u) = OonT; 

AueH. 8 (2.4.7) 

The operator A coincides with —A in -D(A) because, due to (|2.4.1|) and to the 
definition of A we have 

< Au, v >= ((u, v)) = (-Au, v) =< -Au, v >, Vw e V". 

2.4.2 The Operator S. 

We define the form 

b : (H^i?)) 3 -► E 

(«, u, 10) i— > / Ui(diV j)w j dx . (2.4.8) 

Lemma 2.4.1. The form b defined in (|2.4.8|l is trilinear and continuous on the 
product space (H 1 (i?)) 3 . 

Proof. The trilinearity is clear: b(u, v, w) — ((u • V)u, w). For the continuity 
we start to compute | J R Ui(diVj)wj dx\: Since u, v, w G H 1 (i?), by proposition 
12.1.31 we have that 

u h wj € L 4 {R) & diVj € L 2 (R). 
By Holder Inequality we obtain that Ui(diVj)wj S L l (R) and that 

Ui(diVj)u>j dx\ < / \ui(diVj)u>j\dx < \ui\ L 4.\diVj\ L 2\wj\ L A. 



Therefore 

|6(U, V, U))\ < |U| L 4|VW|( L 2)2|W| L 4. 

So, by Corollaries EH and [HOI 

\b(u, v, w)\ < C\\u\\\\v\\\\w\\. 



□ 



The form b, being continuous in (H 1 (i?)) 3 , is continuous in V 3 and, for each 
pair (u, v) <E V 2 we define the operator B(u, v) € V by 

B(it, u):F^M (2.4.9) 

w ^< B(u, v), w >= b(u, v, w) (2.4.10) 

and we set 

B{u) := B(u, u) £ V Vu G V. 



'Since u 6 H 2 , Au e L 2 



2.5. CLASSICAL AND WEAK FORMULATIONS. 
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2.5 Classical and Weak Formulations. 

Classicaly the existence problem for (|1.U.1|) - I|1.U.4[1 . renaming F := F + v 
amounts to finding a vector function 

u : R x [0, T] -> M 2 

and a scalar function 

p : R x [0, T] -> E, 

such that 

tit + (w • V)u + Vp = ^Au + F in i?x]0, T[; (2.5.1) 

V-tt = in Rx]0,T[] (2.5.2) 

u-n = on dRx]0,T[; (2.5.3) 

V ± -?i = on dRx]Q, T[; (2.5.4) 

w(x, 0) = tt (a;) in i?. (2.5.5) 

Where F and wo are given and defined in R x [0T] and i? respectively, uq is 
the position at time of the system. 

Lemma 2.5.1. For u, v € V we have (Au, v) = — ((u, v)). 

Proof. First we note that (Am, v) = — (A^u, A^v) where Ai is the linear map 
defined by 

A^u = A^ u k w k ■= (-k)^u k W k . 



It is clear that A 2 maps V onto H continuously. Indeed for u E V 

ab 

T 



l^l 2 = T E -H-k)ut = Ml 2 < oo 



kern 



and, the continuity follows from the linearity and from \A? (u — v)\ 2 — \\u — v\\ 2 . 
Therefore it is enough to prove the Lemma for u, v £ T>\(R), by continuity it 
will be true for u, v G V. 

2 f 4 2 2 

(Au, v) = / AuiVi dx — — / Vui ■ Wvi dx + / V • (WuiVi) dx 

= -((u, v)) / WiVu-i ■ ndr. 

Now we compute 



wiVMi-ndr= / wi9iuinidr+ / wi9 2 Min2dr. 

The first term vanishes because v\ vanishes where ni does not 9 so we have 
v\S7u\ ■ ndT = + / vid2Ui\i2dT 

JfbUft 



/ 9iu 2 Uin 2 dT = 0. [by V x - it = on T], 



9 ni vanishes on the top and botton faces of the rectangle and vi vanishes on the left and 
right ones. 
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where ft and (below) f r stay for the botton, top, left and right faces of 
R. 

Analogously we conclude that 



v 2 Vu2 ndT = / t^c^wini dT = 0. 

r JflUfr 

Therefore 

(Au, u) = -((it, u)) for u, v £ T>i(R). 

□ 

Lemma 2.5.2. For u, v 6 V we /iawe (Au, u) = — (V • u, V • u). 
Proof. Again it is enough to prove the Lemma for u, v € T>i(R). By (|2.4.(j[) 



/ Au-urfx= / V ± (V ± • u) ■ vdx 
Jr Jr 



=-/V 

JR 


■«)(v x 


• u) cte + 


Jr 


= -/V 

Jr 


•«)(v x 


• i;) eix + 




=- /V 

Jr 


. W )(V^ 


• i>) c&r + 


Li* 


= -/V 

Jr 


■«)(v x 







Where [ Ml := W2 I. □ 



"2/ \— Ml 



Corollary 2.5.3. The norms u i— > ||-u| andu^ |V -u| coincide inDi(R) and 
then in V . 

If u and p are classical solutions of <|2.5.1|) - (|2.5.5[) — u <E C 2 (Rx [0, T]), p€ 
^•(R x [0, T]). Then clearly u <E L 2 (0, T, V) and if we fix v e X>i(iZ), take the 
scalar product and use Lemma 12.5.11 obtain 

u) + v)) + u, u) = (F, u). (2.5.6) 

By continuity the previous expression holds for each v E V. w So is natural to 
define the following weak formulation of problem l|2.5.1|l - H2.5.5(l : 



10 



vie — ► V in || • |[ =>■ V}, — y V in | • |. (u, ■) continuous in | ■ | and ((u, •)) continuous in || ■ | 
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Problem 2.5.1. Given 

F G L 2 (0, T, V') (2.5.7) 

& 

u eH (2.5.8) 

to find 

u G L 2 (0, T, V) (2.5.9) 
satisfying (in the distribution sense) 

— (u, v) + u((u, v)) + b(u, u, v) =< F, v>, Vw G V, (2.5.10) 

and 

u(0)=u . (2.5.11) 

Remark 3. Note that (|2.5.9|l is not sufficient to give sense to ()2. 5.11|l . But we 

will show that if we have in addition ()2.5.10JI then u coincides almost everywhere 
with a continuous function giving a meaning to (12.5.11(1 . 

Now we present a proplem equivalent to Problem 12. 5. II 

Problem 2.5.2. Given 

F S L 2 (Q, T, V'), (2.5.12) 



& 

to find 
satisfying 
and 



u e H, (2.5.13) 
u e L 2 (0, T, V), v! e L\0, T, V') (2.5.14) 
u + vAu + Bu = F on ]0, T[, (2.5.15) 



u(0) = u . (2.5.16) 
To verify this equivalence we will need the following lemmas 

Lemma 2.5.4. Let X be a Banach space with dual X' and let u and g be 

two functions belonging to L 1 ^, b, X). Then the following three conditions are 
equivalent 

1. u is a.e. equal to a primitive function of g, 

«(*)=£+/ g(s)ds, £ G X, a.e.t G [a, b]; 



2. For each test function (f> G T>(]a, b[), 
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3. For each r\ G X' , 



— < u, rj >=< g,i]> 



in the scalar distribution sense, on ]a, b[ 

In particular i/Q£3is satisfied, u is a.e. equal to a continuous function [a, b] — » 
X. 

The proof of Lemma \'2 . 5 . 41 can be found in JHj section 3.1. 
Lemma 2.5.5. If u G L 2 (Q, T, V), the function Bu defined by 

< Bu(t), v >= b(u(t), u(t), v), Vw € V, a.e.t <E [0, T], 
belongs to L\0, T, V). 

Proof. For almost all t, Bu(t) is an element of V. The measurability of the 
function 

[0, T] 9 t ^ Bu{t) e V, 

is easy to check. Indeed it follows from the measurability of 1 1— >• u(t) € V and 
from the continuity of u > B(u) £ V'. By the continuity and trilinearity of b 
on V we have 

\\Bw\\ v > <C\\w\\ 2 , VweV, CeR + . (2.5.17) 

Hence 

\\Bu(t)\\ v > dt < C ||u(<)|| 2 ^ < +oo. 
Jo 

□ 

Now let u satisfy both l|2. 5.9(1 and (|2.5.10() then, due to the identity Q2.4.1[l . 
to definition of A (given in l|2.4.H|l ) and to the previous lemma we can write 
(DJ as 

— < u, v >=< F — vAu — Bu, v >, \fv G V. 
dt 

By the linearity and continuity of A : V — > V we have that 

Au e L 2 (0, T; V') (2.5.18) 
so, F - vAu -Bue L^O, T; V). Therefore, by Lemma |2~CTI 11 and 

f u'ei^O, T; V") 
1 u' = f — vAu — Bu, 

u is a.e. equal to a continuous function [0, T] — > V'. Therefore l|2.5.11|l is 
meaningful and any solution of Problem 12.5.11 is a solution of Problem 12.5.21 
Since any solution of Problem 12.5.21 is clearly a solution of Problem 12.5.11 we 
conclude that these two problems are equivalent. 



'Note that V, being a Hilbert space, is reflexive. 
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2.6 Existence 

2.6.1 Fourier Transform. Fractional Derivatives. 

We want to prove the existence of a solution for Problem 12.5.11 for that we 
will need another compactness theorem envolving fractional derivatives of a 
function. Given a function / from R into a Hilbert space X\ we denote its 
Fourier Transform by 

/+oo 
e- 2 " tT f{t)dt. 
-OO 

The derivative (in t) of order 7 of / is the Fourier Transform of {2%-kt) 1 f or 

DffTr) = (27rr)T/(r). 
Now let Xq, X, X\ be Hilbert spaces such that 

XoClCli 

where, the injections are continuous and the first injection is in addiction com- 
pact. Given 7 > 0, we define the space 

X , Xi) = {ve L 2 (R, X ) I Djv G L 2 (R, X^}. 

H 7 (M., Xo, Xi) in a Hilbert space for the norm || • ||-ht(k,x , x-A defined by 

IMIw-?(K,Xo,x 1 ) = IMIi 2 (R,x ) + IIM 7 £|U 2 (R,Xi)- 
For any set Kctwc define the space 

W£,(R, X Q , Xi) = {ue 7i 7 (R, X , X t ) \ suppu c K}. 

Now we can state the compactness theorem we will need: 

Theorem 2.6.1. Let Xq C X C X\ be Hilbert spaces with both inclusions being 
continuous and the first one being also compact. Then for any bounded set K 
and any 7 > the injection o/W^(R, Xq, X\) into L 2 (R, X) is compact. 

The proof can be found in JS] section 3.2.3. 

2.6.2 The Existence Theorem. 

In the proof of the Existence Theorem 12 . 6 . 41 we will need the following Lemma: 

Lemma 2.6.2. Fixing the first variable in V , the form b defined in H2.4.8(l 
results skew- symmetric in the last two variables, i.e., 

V?i G V\/v, w G H 1 [6(?i, v, w) = —b(u, w, v)]. 

Due to the trilinearity of b the previous Lemma is equivalent to the following 
corollary 

Corollary 2.6.3. Fixing the first variable in V , we have 
Vu G Wd G H^it, v, v) = 0]. 
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Proof. We prove the statement for v S C°°(R). Then by the continuity of b it 
holds for v E H . 

For ueV and v € C°° (R) we have 



fe(u, u, i>) = / UidiVjVj dx — / — ^ — da; 

= E/ v.(u^f) dx -j2[(v. U )^fdx 

3=1 jR 2 J=1 JR 2 

= ±^^-, = ±1^.^ = 0. 

j=i JI j=i J1 



□ 



Theorem 2.6.4. Given F and uo satisfying i|2.5.12|) and (|2.5.13|) . There is at 
least one function u satysfying Ij2.5.14|l - (|2.5.16(l . 

Remark 4. The proof that follows is completely analogous to that for the case 
of full Dirichlet conditions that can be found in ]15^ . 

Proof. We start by defining, for each m S No, an approximate solution u m of 
(|2.5.10|) as follows: 

u m := J2 <® W * (2- 6 - 1 ) 

maxjii , 22 }<m 

{{u m )'{t), Wi) + v((u m (t), Wj)) + b(u m (t), u m {t), W 3 ) =< F(t), Wi >, 

(2.6.2) 

t e [0, T], max{ji, j 2 } < m, 
u m (0) = it™. (2.6.3) 

Where u™ is the orthogonal projection of uq onto span{Wi | maxjii, i 2 } < m}. 
From lj2.6.2|l we obtain the nonlinear system of differential equations in the 
functions it™, max{ii, i 2 } < m: 

E (uTYmwuw^+u Yl «r((Wi,Wi)) 

maxjii , 22 }<m max{ii , 22 }<m 

+ ]T ^r^^^^^.W^, (2.6.4) 

max{ii , 22}<m 
max{/i, ^l^^ri 



2.6. EXISTENCE 



23 



that reduces to the ODE's system 1 

(«?)'(*)(-jt) + ^ 



]T uTu^KWi, W h Wi) =< F(t), Wj >; 



max{ii , *2 

max{ji, j 2 } < m. (2.6.5) 

Now we note that (|2.6.3fl is the same as the m scalar conditions 

u™(0) = the projection of uo onto span{Wj} = uoj. 13 (2.6.6) 

Therefore system (|2.6.5|1 with initial condition l|2.6.6|l has a maximal solution 
defined in [0, t max [ and, we will see that t max — T. Indeed, if t max < T then 
|tt m (i)| should tend to +oo as t goes to t max . Below we show that |w" l (£)| 
remains bounded in [0, T] so we have that the maximal solution is defined in 
[0, T]. 

A priori estimates: We multiply l|2.6.2|) by u™, max{ji, j 2 } < m and add 
the equations obtained. Taking Corollary 12 . 6.31 into account, we arrive to 

{(u m )'{t), u m (t)) + v\\u m {t)f =< F(t), u m (t) > . 

Hence 

d u m (t)\ 2 + 2v\\u m (t)\\ 2 = 2 < F(t), u m (t) > 



dt 



<2||F(t)|| v ,||« m WII<Hl^ l WI| 2 + -|l^)lly' 14 



v 



so. 



±\u"\t)\ 2 + v\\u m m 2 <l\\F{t)\\ 2 v ,. (2.6.7) 
In particular ^\u m {t)\ 2 < ^\\F(t)\\ y, and, integrating on [0, s] we obtain 

\u m (s)\ 2 < W| 2 + - f \\F(t)f v , dt < \u \ 2 + - f \\F(t)f v , dt. (2.6.8) 



v Jq v _ 

Therefore 

the sequence (u m ) remains in a bounded set of L°°(0, T, H). (2.6.9) 
Now we integrate (12.6. 711 over [0, T] and obtain: 

\u m (T)\ 2 + v [ T \\u m (t)\\ 2 dt< |<f + - [ T \\F(t)\\ 2 v , dt 



< | U0 | 2 + - / \\F{t)\\ 2 y, dt. 



In section f2, 31 we have seen that the family W is orthogonal in both H and V. 



12 

13 If we write uo = X] 4eN 2 u 0i Wi. 
14 By expanded (^/ua - -^b) 2 > with a = ||« m (t)|| and b = \\F(t)\\ v , 
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Therefore 

the sequence (u m ) remains in a bounded set of L 2 (Q, T, V). (2.6.10) 
Now we extend u m to the entire real line putting 



u m (t) 



u m (t) iff G [0, T] 
iff 4 [0, T] 



and the Fourier transform of u m will be denoted by u m . 

Application of Theorem 12.6.11 We start by computing the integral 

\ T \ 2 ~<\u m (T)\ 2 dT. (2.6.11) 

-OG 

Equation (|2.6.2(1 with u m replaced by u m results in 

j t (u m (t), Wj) =< f m , Wj > + «\ Wj)6o - (u m (T), Wj)S T ; 

max{ji, j 2 } < rn (2.6.12) 

where 8q, St are the Dirac distributions at and T and, 



f m = F — vAu m — Bu m , 

f m(t) . = !>(*) on[0,T] 

\0 outside[0, T] ' 

Using the Fourier Transform (|2.6.12|) becomes 

2TTtT(u m (T), Wj) =< f m (T), Wj > 

+ (u™, Wj) - (u m (T), Wj) exp(-27riTr). (2.6.13) 

Multiplying (|2.6.13|) by g™"(r) and adding the obtained m equations we arrive 
to 

2iriT\u m {T)\ 2 =< f m {r), u m (T) > 

+ (tig 1 , u m (T)) - (u m (T), u m (r)) exp(-27r i Tr). (2.6.14) 

By (j2.5.17|) we have 

||/ m Wlk' < \\F(t)\\ v ,+v\\u m (t)\\+C\\u m (t)\\ 2 -, so 

\\f m (t)\\ v > dt < f \\F(t)\\ v , +u\\u m (t)\\+C\\u m (t)\\ 2 dt (2.6.15) 
Jo 

by (|2.5.12|l and (|2.6.10|l the integral (|2. 6.15(1 remains bounded. Hence 15 
sup ||/ m (r)|| y / < const, Vm G N . 



5 It is known that |/(t)| < C j R \f(t)\ dt. See for example [T2]. 
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By (|2~Q| both |u m (0)| and \u m {T)\ are finite. Then by (j2A14)l 

\r\\u m (r)\ 2 < C^uTWW+WW < D\\u m (r)\\. 16 (2.6.16) 

Fix 7 < j and define the real function Q(x) := x ^_^ , a; € [0, +oo[. Q is 
continuous and bounded, 17 then we can find a constant D\ £ R + such that for 
all r e M: 

Q(H)<£»i 

l + H 



|r| 27 < Di : 



l + | r |i-2 7 

Therefore the integral (|2.6.11|l is bounded by 

by m 

/+oo \\f.m( \|| r+oo 
i!n$UW- ■*"<*>■**■■ 

The integral L> 3 ||u m (r)|| 2 dr, by the Parseval Equality and l|2.6.10|) . there 
is a constant -D4 such that 

/+00 
\\u m (T)\\ 2 dT < L> 4 , Vra (2.6.17) 
-OO 

For the integral D2 i^ r P^ ^ T we a PP^y Schwartz Inequality and Parseval 
Equality to obtain 

1 1 

ft £° ttf^ * * ft (/T irr*W*j ' (jf "•-»"'*) * 

and, this product is finite and bounded as to — > +00, i.e., there is a constant 
D 5 such that 

r+oc l| fi m( T )|| 



° 2 J_ ]TW^k dT - ° 5 > Vm - (2 - 6 - 18) 

By (12.6. 17|) and (|2.6.18|) we conclude that the integral l|2.6.11|l is finite: 

/+00 
\T\^\u m (r)\ 2 dT<D i + D 5 =:E. 
-OO 

The finiteness of l|2. 6.11)1 with (|2.6.1()|l implies that 

the sequence (u m ) remains in a bounded set of 7i 7 (R, V, H). (2.6.19) 

The limit. We now need the following lemmas that can be found in 0] section 
III.6: 



16 Ci, C2 and D are constants, 
"lim^+oo Q(x) = 1, Q(0) = 



Smce 7 < §, 1+ | T fi-», 6 i 2 (K). Indeed J (1+|T| 1 i- a - T)a = 2/+°° (1+t1 1 -^ )2 - 
Putting x = 1 + r 1— 27 we see that the last integral equals 2 J^ 00 ~2~ 1 _ 1 2 7 ( x — l) 1-27 < 

_Ji_ — 1 anc j ^ ne j as ^ i n t e g ra l converges if 2 — > 1 «-> 7 < ;j- 

~ ~ 1-27 
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Lemma 2.6.5. Let E be a reflexive Banach space and let (x n ) a bounded se- 
quence in E. Then there is a subsequence (aV(?i)) of (x n ) such that x a / n -\ — 1 x, 
for some x G E. 19 

Lemma 2.6.6. Let E be a separable Banach space and let (/„) a bounded 
sequence in E' . Then there is a subsequence (/ CT ( n )) of (/„) such that f a i n ) — S / ; 
for some f G E' . 

Lemma [2.6.61 and (|2.6.9() implies the existence of a subsequence (u aljn ^) of 
u m and u G L°°(0, T, H) such that 

Analogously, Lemma [2.6.51 and H2.6.1U|) implies the existence of a subsequence 

( u a(<r(m))) of u <T( m ) and v g L 2(Q ; Tj guch th&t 

The sequence (u m ) is in the space 7iJ[, T ](M, if) which injection in L 2 (R, if) 
is compact due to Theorem I2.6.1I 21 Then i|2.6.19[l implies the existence of a 
subsequence f3(a(a(m))) of a(a(m)) and w £ L 2 (0, T, if) satisfying 

u /5(a( CT (m))) ^ W; ini 2 (0, T, if). 
We put rj :— (3 o a o a and we obtain 

u = v = w G i 2 (0, T, if) n L 2 (0, T, V) n L°°(0, T, if) 

and 

u v( m ) Wj ini°°(0, T, if); (2.6.20) 

^ Uj ini 2 (0, T, V); (2.6.21) 

u v(m) _^ Uj i nJ L 2 (0, T, if). (2.6.22) 

Indeed i°°(0, T, if) = (L x (0, T, if))' and i 2 (0, T, if) c L l {Q, T, if). So for 
each / € i 2 (0, T, if): 

\ U "M(/) -►«(/). 

Since both w, v £ L 2 (0, T, if) = (i 2 (0, T, if))', we can rewrite the previous 
expressions as 

vf^r 2 m r mvf f( - uV{m)) ~* f{u) 

Vfe(L (0,T,H)) | /(u „ (m)) 
Then u = u G i 2 (0, T, if). 

On the other side (i 2 (0, T, if))' C (i 2 (0, T, V))' because a converging se- 
quence in L 2 (0, T, V) converges in L 2 (0, T, if) too. 
So for each / G (i 2 (0, T, if))' we have 

f/(^M) -►/(„) 
1/(u"M) -»■ f(w). 



19 We shall use the symbols: — »• for weak convergence; — for weak-star convergence and, 
•» for strong convergence. 
20 If wanted, without lack of generality we may assume u m — ^* u, inL°°(0, T, H). 
21 With V = X , & H = X = X l . 
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But since both w,v & L 2 (0, T, H) we have that v = w. 

Ending. To complete the proof we will need the following Lemma: 



Lemma 2.6.7. If a sequence (u n ) satisfies u n — u in L 2 (0, T, V) and u n — > u 
in L 2 (0, T, H), then for any vector function w with components in C (R x 

[o,m 

b(u n {t), u n {t), w(t)) dt~> [ b(u(t), u(t), w(t)) dt. 



Proof. From 



b(u n , u n , w) - b(u, u, w) = b{u n - u, u n , w) + b(u, u n - u, w) 



we obtain 



b(u n (t), u n (t), w(t)) -b(u(t), u(t), w(t))dt 



< \b(u n (t)-u{t),u n (t),w(t))\dt+ \b(u(t),u n (t)-u(t),w(t))\dt 



(using l|2.7.1|l below, and Lemma \'2. 6. 21 ) 

i-T r 



<c 



\u n (t) - u(t)\-*\\u n (t) - u(t)p\u n (t)\?\\u n (t)p\M\ 



+ \u(t)\-^u(t)p\u n (t)-u(t)\^\\u n (t)-u(t)p\\w\\ 



dt 



<C X U \u n (t) - u(t)\\\u n (t) - u(t)\\ dty\\u n {t)\\ LW>v) 

+ Ci||«(*)Hi»(o,r J To(y \u n (t)-u(t)\\\u n (t)-u(t)\\dty ; 
but u n is bounded in L 2 (0, T, V), because u n u in L 2 (0, T, V), then 



b(u n (t), u n (t), w(t)) - b(u(t), u(t), w{t))dt 



<C 2 \\u n {t)~u{t)\\l 2 ^ 



T, H)' 



Hence \ J b(u n (t), u n (t), w(t)) — b(u(t), u(t), w(t))dt\ goes to as n goes to 
+oo. □ 



22 Note that since both L 2 (0, T, V) and L 2 (0, T, H) are both Hilbert spaces then each of 
them coincide with its dual. But we can not do these identifications here if we want to compare 
their duals, because the spaces H and V the norms are different. If we do the identifications 
we will obtain (L 2 (0, T, H))' D (L 2 (0, T, V))' and that is not true. 
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Multiplying (|2.fi.2Jl by a function <f> £ C' 1 ([0, T]) and, integrating we obtain: 

((u m )'{t), Wjtftt)) dt+ [ v((u rn {t), Wj(j)(t))) dt 
Jo 

+ I b(u m (t) 7 u rn (t) 7 Wj4>{t))dt = I < F(t), W s <l>{t) > dt 



(u rn (t), Wj(f)'(t)) dt + v ((u m (t), Wj4>{t))) dt 



+ f T b(u m (t), u m (t), W j( p(t))dt = «, W,)#t)(0) - ( W m (T), W,)^(t)(T) 
Jo 

+ / < F(t), Wj<f>(t) > dt. 
Jo 

Due to Lemma r2.6.7l we can take the limite in the nonlinear term and, by Ij2.6.21|l 
we can take the limite in the linear terms. 23 Hence 



'j 

rT 



(u(t), Wj(j>'{t)) dt + v ((u(t), Wj(j){t))) dt 
Jo 

b(u(t), u(t), Wj(f>(t))dt 
{u , Wj)<P{0) - (w(T), WjWT) + [ < F(t), Wj<P{t) > dt (2.6.23) 



o 



Equation l|2.6.23[) being true for all Wj by linearity will be true for any finite 
combination of functions in W and by continuity will be true for all v E V. 
Then taking a test function <fi E 2?(]0, T[) in (|2.6.23(l we conclude that 

{u{t),v<j>'(t))dt + v [ ((u(t), v(j)(t)))dt+ ( b(u(t), u{t), v<j>(t))dt 



=(«o, «)0(O) - («(T), v)(j)(T)+ < F(t), v<f>(t) > dt; (2.6.24) 
then 



(u(t), v(j)'(t)) dt + v {(u(t), v(j)(t))) dt+ b(u(t), u(t), v</>(t)) dt 
Jo Jo 

= < F(t), v<t>(t) > dt (2.6.25) 
what means that 

Equation H2.5.1()(l is satisfied in the distribution sense. 
Multiplying (Ii>.5.l0j> by <p E C l ([0, T]) such that </)(0) = 1, cf>(T) = we obtain 

(u(t),v(j)'(t))dt + v [ ((u{t),v<f)(t)))dt+ [ b(u(t), u(t), vcj)(t))dt 



--(u(0), v)+ j < F(t), Wj<f>{t) > dt. (2.6.26) 



J Note that implies that u^" 1 ) -v u in L 2 (0, T, H). 
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The equation resulting from l(2.fi.24(l with the same <f> as in l(2.fi.2fi(l is 

rT pT pT 

(u(t), vcf>'(t)) dt + is ((u(t), v(f>(t))) dt + b(u(t), u(t), V(j){t)) dt 
'o Jo Jo 

=(u , v)+ < F(t), vcj>(t) > dt (2.6.27) 

From Ij2.fi. 26(1 and 1(2. fi. 27(1 we conclude that 1(2.5.11(1 is satisfied (and, we have 
finished the proof of theorem l2.fi. 4(1 . □ 

2.7 Uniqueness. 

Lemma 2.7.1. 

\b(u, v, w)\ < c\u\i\\u\\i\\v\\\w\?\\w\\%, Vu, v, w £ H^i?). (2.7.1) 
Ifu belongs to L 2 (Q, T, V)nL co (Q, T, H), then Bu belongs to L 2 (0, T, V) and 

\\Bu\\ L ^ TjYt) < c\u\ L ~ (0 , T , H) \\u\\L 2 (0, T, V). (2.7.2) 

See [HI section III 3.2. for l(2~7~T|) . Then lt2~7^|) is a corollary of l|27rTTJl . 

Indeed 

\b(u, u, w)\ = \b(u, w, u)\ < C|tt|4||u||i|u|*||«||i||«;|| 

= C|tt|H|H, w, ueV. (2.7.3) 

In the case u G L 2 (0, T, V) n L°°(0, T, H) we obtain that 

\\Bu{t)\\ 2 v , dt < C [ \u\ 2 \\u\\ 2 dt 
Jo 

<C\\u\\lao, 0:T H) / \\u(t)\\ 2 dt = C\\u\\l x , TtH) \\u\\l 2(0 T V) < +00 
Jo 

Theorem 2.7.2. The solution of Problems \2.5.1l2.53 given by Theorem \2.b'.4\ 

is unique. Moreover it is a.e. equal to a continuous function from [0, T] into H 
and, 

u(t) -> u(h), inH ast^h t\ G [0, T]. {2.7 A) 

In particular 

u(t) — > uq, in H as t — > 0; 
u(t) — > u(T), inH as t — > T. 

For the proof we need the following lemma from |15( : 

Lemma 2.7.3. Let V, H, V' three Hilbert spaces satisfying V C H C V with 
dense and continuous inclusions. If 

u S L 2 (0, T, V) & u' e L 2 (0, T, V) 

then u is a.e. equal to a continuous function from [0, T] into H and 

^\u\ 2 = 2 < it', u > (2.7.5) 
holds in the distribution sense on (0, T). 
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Remark 5. Note that (|2.7.5[1 is meaningful because both 

1 1 y \u{t)\ 2 & t ^< u'{t), u(t) > 
are integrable on [0, T\. 

Proof of Theorem \m% By ()2.5.12j) . 12.5.15)1 . (|2.5.18|> and J^2J wc have that 
u' G L 2 (0, T, V')- Then by (|2.6.20j) . (|2.6.21|) . and Lemma l2~7~31 we have that 

ueC([0, T], if). 

Therefore (I2.7.4|l is satisfied. 

To prove the uniqueness we consider two solutions u, v of the problems and put 
w := u — v. Then 

w =u — v = F — vAu — Bu — {F — vAv — Bv) 
= — vAw — Bu + Bv 

and 
w(0) =0. 

Now we take the scalar product with w in the duality between V and V' and 
obtain 

< w' ' , w > +v < Aw, w >=< Bv, w > — < Bu, w > . 

We note that since u, v G L 2 (0, T, V) & u', v' € L 2 (0, T, V) then w G 
L 2 (0, T, V) & w' G L 2 (0, T, V')- Hence we can apply (|2~7"5f and obtain 

^-Hi)| 2 + 2^||w(t)|| 2 = 2b(v(t), v(t) w(t)) - 2b(u(t), u(t), w(t)) 

= -2b(w(t), v(t), w(t)) = 2b(w(t), w(t), v(t)). 24 

Hence using ()2.7.3|l 

IK*)! 2 + 2v\\w(t)\\ 2 < 2(CH*)IIK*)||IK*)II) 

<2v\\w(t)\\ 2 + ^\w(t)\ 2 \\v(t)\\ 2 25 

Threfore 

^Ht)\ 2 <^\ w (t)\ 2 \\v(t)\\ 2 . 

at 2v 

Since t i-> ^|w(i)| 2 = 2 < tt/(t), io(t) >, £ i-> Ht)| 2 =< w(t), w(t) > and 
t h- > ||w(i)|| 2 =< > are integrable, by Gronwall Inequality (see l|2. 7.6(1 

below) 




24 b(w, v, w) = —b(v, v, w) + b(u, v, w) = —b(v, v, w) + b(u, w, w) + b(u, u, w) = 
—b(v, v, w) + b(u, u, w). 
25 By expanded {yj2va j=bf > with a = \\w(t)\\ and 6 = |to(t)|||v(t)[|. 
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(using the fact that w(0) = 0). Hence 

u = v 

concluding the uniqueness of the solution. □ 

Lemma 2.7.4 (Gronwall Inequality). Let <?, h, y, be locally integrable 
functions sattisfying 

^<gy + h fort>t . (2.7.6) 

Then 

y(t) < y(t )exp( / g(T)dr\+ h(s)exp(- g^drjds, t>t . 
[The proof can be found in 

2.8 Continuity on Initial Data. 

Theorem 2.8.1. The map 

§ : H x L 2 {0, T, V')x]0, +oo[ -> C([0, T], ff) 
(u , F, v) 

is continuous. Here u £ C([0, T], If) is the unique solution of Problems \2.5Al 
WT^ (see Theorem \F71fy . 

Proof. Fix a triple (u , F, is) E H x L 2 (Q, T, V')x]0, +oo[ and consider the 
solution of Theorem 12.7.21 induced by this triple. Then u satisfies 

v! + vAu + Bu = F, u(0) = u . 

Now consider another triple (vq, G, 77) € H x L 2 (0, T, V)x]0, +oo[. The solu- 
tion associated with this triple satisfies 

v' + r]Av + Bv = G, v(0)=v . 

We put 

w := v — u 
and see that w satisfies the equation 

w' = G - F - r/Aw + (v- r))Au - Bv + Bu. 
Taking the scalar product with w we obtain 

< w' , w > —< G — F , w > — ?7||w|| 2 + {y — rf)((u, w)) + b(w, w, u) 

j t \w\ 2 < 2\\G - F\\ v ,\\w\\ - 2 V \\w\\ 2 + 2\y r,\h\\\W\ + 2C\w\\\w\\\\u\\. 
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Expanding 



(yfiHi-cy^NiHi) 2 >o; 

(Jl\\M\-W~v\J^\H\) 2 >0; 



we obtain 



^M a + filMl a < ^HGf-#||L + |i/-n| a ^||«|| a + o a ^M a IMI a . (2.8.1) 

(JiL If II II 

By (p.7.6|) we have 

K*)l 2 

,3 " 



< |w(0)| 2 exp(c 2 ^y ||u(r)|| 2 dr) 

C l / 3 3 \ / 3f 2 ( s \ 

+ (-IIG^-F^II^ + l^-^-IK^fjexpf / \\u(r)fdr)ds 

Jo y V V ' v V Jt ' 

< exp(c 2 ^ T ||u(r)|| 2 dr) (jw(0)| 2 
+ f T (-\\G(s) F(s)f v , + \u- v\ 2 -\\u(s)\\ 2 ) ds^ 



h G (s)-F(s)\\ 2 v , + \v-r,\ 2 -\ 
o Kr ) V 

Now if \v - ?/| < |, i.e., rj G]|, we have 



Ki)| 2 <exp(c 2 ^ jf || U (r)|| 2 dr)(^K0)| 2 

+ f T (l\\G( S ) - F(s)\\ 2 v , +\u- v\ 2 rM s )^) ds 



Fix e > 0. 

Now if we put Eq :— exp^C 2 ^ J Q T ||u(t)|| 2 drj and choose the triple (vq, G, rj) 
such that 

£ -I 

|wo - «o I = k(o)| < -7=^0 2 = :a «; 

1 

||G-F|| £a(0 ,T,v')= f^HGW-J'WII^dsJ < ^l^* =:a fl ; 

77 1 < min{^, 5} =:a v , 



(/ \\u(s)\\ 2 d s y~ 2 E- 



where 5 = ^ 
3^ 

we obtain 

|77j(i)| 2 < e 2 \w(t)\ < e, t e [0, T] 



2(i 



26 Since u £ L 2 (0, T, V) we have that ot v and a v are finite. We can also see that ot v , 
depend only on the (fixed) triple (uq, F, v). 
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and, we have the continuity of §. □ 
Theorem 2.8.2. The map 

§ 2 ■ H x L 2 (0, T, V')x}0, +oo[ -> L 2 {0, T, V) 
(uq, F, v) h- > u 

is continuous. Here u € L 2 (0, T, V) is the unique solution of Problems \2. 5 Al 
WT^ fsee Theorem \£7J$ . 

Proof. As in the proof of Theorem 12 . 8 . II we fix e > and a triple (uq, F, v) S 
H x L 2 (0, T, V')x}0, +oo[ and consider the solution of Theorem 12 . 7.21 induced 
by this triple. Now consider another triple (vq , G, if) e HxL 2 (0, T, V')x]0, +oo[ 
and the solution associated with this triple. Put 

w := v — u 

and taking the scalar product with w we arrive to (|2.8.1|l . Integrating over [0, T] 
we obtain 

-\w(T)\ 2 - -\w(0)\ 2 + f \\w{t)\\ 2 dt <\\\G- F||| 2(0 , Tj v/) 
V V Jo V 2 ( ' ' 

3 3C 2 
+ W-v\ 2 ^\\u\\ 2 L 2 {0tTi V) + -^a-||Hlo([o J r],H)Hli«(o I r,v) 

and, if \v — r]\ < \ we have 

\\w(t)\\ 2 dt< ^\\G-F\\l 2{0tTiV , ) + \u-r ) \ 2 ^\\u\\l H0)T<v) 



+ [ v 2 \\ U Wl 2 (0.T,V) + - )\\ w \\c(0,T,H)- 

Using Theorem l2.8.1l there is 5 > such that if (vq, G, if) satisfies 
\v - u \ < S, \\G - F\\ L 2( T VI) < S, \i] -v\<6, 

then \\w\\ c( p,T,H) < (4)*(^Hl£ a(0 ,T,v) + 
Hence for 

we have ||u>||£2(o,t, v) < £ - Therefore the map §2 is continuous. □ 

2.9 Some More Estimates for the Form b. 

We present here some estimates we shall need later. We have 27 

\b(u, v, w)\ < CK. 

27 These estimates can be found in |13l section III. 3. 2.. They can be obtained by interpola- 
tion (El), generalized Sobolev inequalities ( IWl ) and by a theorem by S.Agmon ([Q)- See 1141 
for indications how to obtain them. 
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where C is a constant and K is one of the following products 



Ml IMI II w \\ 

u|*H*IMI*IMltM 

"I'HIalMIM 
"INIM'IMIa 

u|*IMI^IMIM^IMI 5 



u, v, w G H 1 ^), 

u G H\R), v G H 2 (i?), 10 G L 2 (i?), 
u G H 2 (i?), v G H^i?), 10 G L 2 (i?), 
u G L 2 (i?), i! G H x (iJ), w G H 2 (i?), 
u, u, u> G H x (i2). 



(2.9.1) 
(2.9.2) 
(2.9.3) 
(2.9.4) 
(2.9.5) 



2.10 Strong Formulation. 



Sometimes we want more regularity for the solutions of problem (|l.().4|l . 
Instead of Problems (|2.5.1[l - H2.5.2Jl where we ask for weak solutions we consider 
the following equivalent problems where we look for solutions more regular than 
weak. 



Problem 2.10.1. Given 

F G L 2 (0, T, H) 



k 

to find 



u a eV 



u G L 2 {0, T, D(A)) n L°°(0, T, V) 
satisfying (in the distribution sense) 



— (u, v) + v{{u, v)) + b(u, u, v) — (F, v), Vu G V, 



and 



it(0) = Uq. 
Problem 2.10.2. Given 

F G L 2 (0, T, H), 



k 

to find 

satisfying 
and 



u a G V, 

u G i 2 (0, T, D(A)) n L°°(0, T, V), and 
v! G L 2 (0, T, if) 

it' + z/Au + Bu = F on ]0,T[, 



u(0) = u . 

The equivalence of these problems follows from 

1. A solution of Problem 12. 10.21 is a solution of Problem l2.1U.ll 



2.10.1 
2.10.2 
2.10.3 

2.10.4 
2.10.5 

2.10.6 
2.10.7 

2.10.8 
2.10.9 
(2.10.10 
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2. For a solution u of Problem 12 . 1 . II we have 
\Bu(t)\ < C\\u(t)\\?\\u(t)\\$ a.e. t by (223, so 

|iM*)| 4 < C / T ||«(t)|| 6 ||«(t)||l < Ci f \\u{t)\\l < oo. 
Hence 

(a) G i 4 (0, T, H). 

Since it G L 2 (0, T, and J T |Au(<)| 2 = J Q T < oo we 

have 

(b) Am G L 2 (0, T, if). 

From Il2al and Ij2bj) we have that 

/ - vAu - Bu G L 2 (0, T, H) C L 2 (0, T, F')- 

Since u € L 2 (0, T, D{A)) C i 2 (0, T, V), by Lemma (l23Hl and <2.1().4t . 
we have that 

u = f - vAu - Bu a.e. and u € (7(0, T, l/')"a.e.". 

Hence 

3. m is a solution of Problem l2~TO 
Lemma 2.10.1. The inclusions 

D(A) QV <ZH 
are both dense continuous and compact. 

Proof. We have already seen in Corollarv l2 . 1 . 91 and beginning of section l2~2l that 
the inclusion has the required properties. For the first inclusion we have 

that V, respectively D(A), are the closure of T>i(R) in H 1 , respectively in H 2 . 
From the density, continuity and compactness of the inclusion H 2 C H 1 come 
the same properties for the inclusion D(A) CV. □ 

2.11 Existence. 

Theorem 2.11.1. Given F and u satisfying (|2.10.6|l and 12.10.7j l. There is 
at least one function u satisfying ()2.10.8jl - lj2.10.10jl . 

Remark 6. The proof that follows is completely analogous to that for the case 
of weak solutions and for no-slip boundary conditions that can be found in }1<^ . 

Proof. Outlines: We define an approximate solution u m , for each m G No like 
in H2.6.1|) - (|2.ti.3|l . and arrive to the estimate l|2.6.7jl and conclusions l|2.6.9|) 
and (|2.6.10j) exactly in the same way. For Theorem we need only some more 
estimates: If we multiply 12 . 6 . 21 bv juj 1 and add the obtained equations we arrive 
to 

({u m (t))', Au m (t))+v({u m (t), Au m {t))) + Bu m (t)(Au m (t)) = (F(t), Au m (t)), 



36 CHAPTER 2. EXISTENCE, CONTINUITY AND UNIQUENESS. 



{{{u m {t))', u m (t))) + v(u m (t), u m (t)) [2] +Bu m (t)(Au m (t)) = (F(t), Au m (t)), 
>r 

i|lk m (t)|| 2 + ^|ii m Wlf 2] <|^ m W(^-(i))| + |F(*)||A U -(t)| 

(by EH3) < q^(i)|i|K(i)||| U ™(t)|f 2] + \F(t)\\u m (t)\ m 
(by Young Inequalities) 28 

<^(l« fn (*)l| ] )*+%|(C|« Tn (t)|*||u' n (t)||) 4 

+ ^K(t)|f 2] +C f , 2 |F(i)| 2 

< Jd^Wll) 1 + ^(CKWI'll« m WII) 4 

+ ^ m (*)l[2] + ^(*)! 2 - 



||K"|| 2 + H« m |f 2] < +|^C 2 | U m (i)| 2 ||u™(t)|| 4 + ^|F(t)| 2 . (2.11.1) 

From equation l|2.11.1[l and from Gronwall Inequality we can derive for s G [0, T] 
(using lf2~Ojl and l|2.6.10[) ) 

\\u m (s)f 

<exp(^ T |^C 2 | M m (0| 2 || U m (i)|| 2 *)(||<|| 2 + ^ T ^|F(t)| 2 ^ 

for some constant K\ (independent of m). Hence 

the sequence (u m ) remains in a bounded set of L°°(0, T, V). (2.11.2) 
Now we integrate (|2.11.1|) over [0, T] and obtain: 

|K l (T)|| 2 -|K(0)|| 2 + i , [ T \u m (t)\f 2] dt 

Jo 

-[ £ c V i wi 2 ik"wii 4 + l\nt)\ 2 dt < k 3 

with K3 being a constant (independent of m, using G32a, una and em). 

Therefore 

the sequence (u m ) remains in a bounded set of L 2 (0, T, D{A)). (2.11.3) 

28 Young Inequality: ab < eaP + CY tP b p ' . Where a, b > 0, 1 < p < oo, £ > 0, A + i = 1 and 
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By 12.11.21 12.11.31 and Lemmas 12.6.51 and 12.6.61 we conclude the existence of a 
subsequence such that 

(™) _l „ in T? 



u y ' — u 



in L 2 (0, T, D(A)) 



u *(m) ^ U in L°°(0, T, V) 
This limit is a solution for problems l2~l 0.112. 10. 21 □ 



2.12 Uniqueness. 

Lemma 2.12.1. The inclusions 

D(A) C V C £>(A)' 
are oof/i dense and continuous. 

Proof. We have just seen in Lemma fe.lO.ll that the inclusion -D(A) C V has the 
required properties. For the second inclusion we proceed as in the beginning of 
section l2~2l We identify V with V' 29 The continuity of the second inclusion 
follows from the continuity of the first one. The density of the second inclusion 
follows from Lemma T2.2. II and the injectivity of the first one. □ 

Theorem 2.12.2. The solution of vroblems \2. l(Tll2. 10.21 is unique. Moreover 
it is a.e. equal to a continuous function from [0, T] into V . 

Proof. A solution of problems 12 . 1 . 112 . 1 . 21 is a solution of problems 12.5. 112.5. 21 

as well, which is unique by Theorem l2~7~2l By u G L 2 (0, T, D(A)), 

u G L 2 (0, T, H) C L 2 (0, T, V') C L 2 (0, T, D{A)') 
and by lemmas [2~T2~T1 and l2~7~3l we conclude that u £ C([0, T],V). □ 



2.13 Continuity on Initial Data. 

Theorem 2.13.1. The map 

§ s : V x L 2 (0, T, i?)x]0, +oo[ -> C([0, T], 7) 
(u , F, v)^ u 

is continuous. Where u G C([0, T], V) is i/ie unique solution of Problems \2. 10.H 
\2.10.2\ (see The.orem, \2.12.2\) . 

Proof. Fix e > and a triple (w , F, u) G V" x i 2 (0, T, i?)x]0, +oo[ and, 
consider the solution of Theorem 12 . 1 2 . 21 induced by this triple. Then u satisfies 

u' + vAu + Bu = F, u(0)=u . 

Now consider another triple (u , G, rf) € V x L 2 (0, T, H)x]0, +oo[. The solu- 
tion associated with this triple satisfies 

v' + rjAv + Bv = G, v(0)=v . 

29 Note that here H is not present. We are doing only one identification. 
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We put 

w := v — u 
and see that w satisfies the equation 

w' = G — F - r/Aw + (v - ij)Au - Bv + Bu. 

Taking the scalar product with Aw we obtain 

-||H| 2 + 2r?Mf 2] < 2{G - F, Aw) + 2\v - r)\(u, w) [2] + 2(Bu - Bv)(Aw). 
Choosing 77 such that \v — ?/| < ^: 

j t \\w\\ 2 + v\w\ 2 [2] < 2(G-F,Aw)+2\v-r]\[u,w) m +2(Bu-Bv)(Aw). (2.13.1) 

The last term is the more complicated so we work it a bit: 

(Bu — Bv)(Aw) = b(u, u, Aw) — b(v, v 7 Aw) 

= b(u — v, u, Aw) + b(v, u, Aw) — b(v, v, Aw) 

= —b(w, u, Aw) — b(v, w, Aw) 

= —b(w, u, Aw) — b(w, w, Aw) — b(u, u>, Aw) 

By f2~<Ei|) we have 

\b(w, u, Aw)\ < t7|«;| * ||w|| * ||u|| * [a] 

< Ci||u)|||m|[ 2 ]|w|[2] 

and, by lTHH|) 

1 1 1 1 

\b(w, w, Aw) + b(u, w, Aw)\ < C\w\ 2 |w| ^ |MI \ w \[2] + C\u\^\u\ 2 

< CH'HtflHHl] + C|u|| ] |tt|i||«;||Hp] 

Hence from <|2.13.1[) we have 

^-\\w\\ 2 + v\w\f 2] 
<2\G - F\\w\ [2] + 2\v - 77||u| [2 ]|w|[2] + 2Ci|MN[2]Hp] 
+2C|io|i||«;|||ii;| [ | ] + 2C|u|| ] |u|*|MIM[2]. 
Applying Young inequalities with suitable exponents and constants we arrive to 

^IHI 2 + HHf 2] 

+^\™\h + ™ci\\wr\u\f 2] 

+ (2C) 4 M 2 IHI 4 + ^Mf 2] + ^c 2 mm [2] ihi 2 - 
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Then 

<H| G _ p\ 2 + % - V \*\u\h + -c!\\wf\u\h 
+D 1 \w\ 2 \\w\\ 4 + -C 2 \u\\u\ m \\wf, (2.13.2) 

where D 1 = (f ) 4 (v) 3 - B y Gronwall Inequality fFUfy we have for t € [0, T] 
IK*)I| 2 

< exp ( f -Cl\u{s)\\ 2] + Dxms)! 2 !^^)!! 2 + ^-C 2 \u(s)\\u(s)\ [2] ds) 

• (ll^(0)|| 2 + jf ™\G(s) F(s)\ 2 + % - rilM^] ds ) ■ ( 2 - 13 - 3 ) 

By theorems E33 and EES and from u £ L 2 (0, T, D(A)) D L°°(0, T, V) the 
argument of the exponential is bounded, say less than a constant E, if we choose 
the triple (vo, G, rj) such that both |uo — uq\, \\G — -F||l 2 (o ,t, v) an d W ~ 7 l\ are 
less than S for some S > sufficiently small. 
Now put 

Si :=min{^, J, ^(exp E)~i , ^=(^7) '(exp^TT^ 

1 

' (exp£)~5 j. 



-LfHllnll? 
V3^ 



u lli 2 (0,T, D(A) 



It is clear from (|2.13.3JI that if both |«q — uq\, \\G — F\\ L 2m T ,v') an d |f — rj\ are 
less than Si, we have ||u>(i)|| less than e and, we have the continuity of S s . □ 

Theorem 2.13.2. The map 

S 2s : V x L 2 (0, T, i?)x]0, +oo[^ L 2 (0, T, L»(A)) 
(u , F, v)^> u 

is continuous. Where u G L 2 (0, T, D(A)) is the unique solution of Problems 
\£THJSKIM (see Theorem\nEM). 

Proof. As in the proof of Theorem 12 . 1 3 . II we fix e > and a triple (uq, F, v) S 
H x L 2 (0, T, V) x]0, +oo[ and consider the solution of Theorem l2 . 1 2 . 21 induced 
by this triple. Now consider another triple (vo , G, -q) e HxL 2 (0, T, V')x]0, +oo[ 
and the solution associated with this triple. Put 



w := v — u 



and taking the scalar product with Aw we arrive to H2.13.2(l . Integrating over 
[0, T] we obtain 
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T 

^IMI[2] 



MT)|| 2 -M0)|| 2 + 

< — - HI L 2 (0, T, H) + —\\ u \\ 2 L 2 (0,T, D(A))W ~ V\" 

+ ~^l\\ u \\ 2 L 2 (0,T,D(A))\\ w \\c([Q,T],V) + TD i\\ w \\c([0,T], H)\\ w \\c([0,T],V) 

+ ~^ 2 \\ u \\c([0,T], H)\\ U \\1 2 (0,T, D(A))\\M\c([0,T],V)- (2.13.4) 

By Theorem 12.13.11 for a small enough S > we have that if both \vq — uo\, 
\\F ~ ^IIl 2 (o, t, h) > an d \t) — v\ are less than S, then ||Hlc([o,T], V) < 1- Thus 
from (|2.13.4[) we have 



k" 2 



i) 



[2] 



20 - 2 20 2 2 

— ^2 11^ _ f Wl 2 (o,t,h) + ^2"!MIl 2 (0,T, DfA))! 1 ' ~~ v\ 

v ~ + ^2^lll M lli 2 (0,T, £>(A)) + ~ T D l\\ w \\c(lO,T], H) 

+ ^2 ( - ; ' 2 |l u llc'([0,T],H)ll u lli 2 (0,T,_D(yl))) II HI C([0, T] , V) 

20 

<^II G " HI^(0,T,ff) + ^ " ^ + A»||Hlc([0,T|, V) 

for some constants D2, D3 and all triples (vq, G, rj) satisfying 
\v Q - u Q \ < 5, \\F - G|| L 2( 0)T)J? ) < S, \r) -v\<S. 

Then for some Si smaller than 5 we have that if both \vo — uq\, \\F — G'||l 2 (o, t, h) > 
and \rj — v\ are less than S\, we obtain 



(/„ »-»f 



< e. 



I[2] y 

Therefore the map S2 S is continuous □ 



Chapter 3 

Controllability of Galerkin 
Approximations. 

3.1 The FCE Procedure. 



In this section we present the FCE Procedure, i.e., a procedure of Factoriza- 
tion^ Convexification+ Extraction: 

3.1.1 A Lemma from Linear Algebra. 

A result from Linear Algebra we will need is the following: 

Lemma 3.1.1. Fix a linear space X of dimension N > n + m. Given two 
families V :— {vi \ i = 1, . . . n} and W := {wj \ j = 1, . . . m} satisfying: The 
vectors in V are linearly independent and the vectors in ttW = {ftWj \ Wj G W} 
are linearly independent as well, where tt is the projection onto some space V c 
transversal to span(V), (X = span(V) ®V C ). Then the family V^JW is linearly 
independent. 

3.1.2 Factorization. 

Consider a control-afHne system 



where /, gi are smooth vector fields and [<?,, gj] = i = 1, . .. , r. 

In |2j it is proven that if we decompose the flow of system as 



r 




(3.1.1) 



i=l 
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where g := (gi, g%, ■ ■ ■ , g r ), v '■= i v i: V2, • • ■ , v r ) T and G* v m stays for the flow 
exp f Q gv(r)dT = e gw ^\ w(t) — J v(r)dT, then 

A qo {f + gv){t) = A qo ((e-9V)*f)(t) o {Gl {t) I v{t) i W] 1 (3.1.4) 

Here A x (y) (t) stays for the attainable set at time t from x following the vector 
fields y. 

Similarly if we rewrite system (|3.1.1|l as 

r 

Q = /(«) + E( w ' (*) + v i(t))9i(q) Q G R n , G K 

i=l 

we arrive to 

Ao(/ + ^)W = Ao((e-» v2 ),/i)W o {G* 2(t) I « 2 (t) i K r } (3.1.5) 

where /i((j) := /(<7)+X}I=i v l(^)9i(q) an d where u , i> 2 and F 2 are independent. 
The system q = (e~ gV )*fi(q) is called factorized system. 

Lemma 3.1.2. With (e~9 y )„f 1 and G* 2 as in equation (|3.1.5j) we have 

Ao(M F2 )./i)W°{^ (() h 2 (t)er} 

3A.(Mn,/i)W°fc (t) I v 2 (t) e K r }. 

Proof. Let x G „4 gn (( e -9 y2 )*/i)(<) ° {G* 2(t) | w 2 (i) e R r }. Then there are 
y £ ^ ((e-s y2 )*/i)(t) and a control u(t) € R r f e [0, T] such that a: = yoG^uy 
Let y„ — ► y, y n £ A qo ((e~ gV2 )*fx){t). Hence x„ = y„o G^ (t) is a sequence 
on ^^((e -9 ^ )*/i)(t) ° G^m t nat converges to x. □ 

So system l(3.1.1|) is approximately controllable in time t if 

A q {{e-s v2 )*h){t) o {G* 2(t) I v 2 (t) g R r } = R" Vo e R". (3.1.6) 

If are constant vector fields — 9i{q) — 9i i = 1, ■ ■ ■ r q £ M. n — they 
commute and the systems q = (e~ 9X ^)*fi(q) and q — fi(q + gX(t)) coincide. 
A corollary of this is 

Corollary 3.1.3. System l(3.1.1|l (with g constant) is approximately controllable 
in time t if 

A q (f lx )(t)o{e3 v2 ^} = W n Vgel". 2 

Where 

Ax(q) := fl(q + gX(tj) = f(q + gX(t)) + g(q + gX^v 1 
= fx(q)+gv 1 X(t)£R r . 

In particular it is approximately controllable in time t if 

A q (.hx)(t) = R n Vq£R n . 

1 Here V is independent of v. 

2 Here, for more precision, we should write { e s v (*) | V 2 (t) El r ). 
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3.1.3 Convexification. 

If for some constant vector 7 € K", f(q) +7 belong to the convex set Conv{fx \ 
X e R r }, then for every u 1 e R r 

/(g) + 7 + 5U 1 € Cont;{/ x | X € E r } + span(g) 

C Conv{f x + gw 1 I.d'g R r }. 3 

This means that we can follow any of the vector fields f(q) + 7 + gv without 
changing the closure of the attainable set at time t (recall that convexification 
does not change the closure of attainable set at time t see 8 ). In particular 
system (|3.1.1ll is approximately controllable at time t if 

A q (f(q)+ 1 + gv^)(t) =R" VqeW 1 . 

3.1.4 Extraction. 

Let G be a cone and suppose that 

f(q) + CC Conv{f x I X G M r }. 

Then putting G := span(g) 

f(q) + C + G C /(g) + CW(C) + G C Conw{/ x (g) + G | Iel r } 
= Com;{/ 1Jf ( g ) I I6f}. 

Now from Conv(C) + G we extract the linear space 

G 1 := (G + Gonw(G)) D (G - Conv(C)). 

We shall call the directions on G 1 "extracted" directions. Since clearly G C G 1 
because € G, those directions in G will be called "old" directions and, those 
in G 1 \ G "new" directions. 

Adding new directions does not change the closure of attainable sets so, we can 
say that system (13.1. It is appoximately controllable in time t if the "bigger" 
system q = f(q) + giv 1 is, where v 1 € W 1 , r\ (> r) is the dimension of G 1 and 
<7i is a matrix whose r\ columns are vectors spanning G 1 . 

3.1.5 Iterating FCEs. 

Iterating FCE Procedures we obtain an increasing sequence 

G =: G° C G 1 C ... C G j C ... 

of subspaces of controlled directions without changing the closure of the attain- 
able set at time t. Obviously if for some p E N we have G p — R™, then the 
controllability in time t is an immediate consequence of Corollary 13.1.31 (note 
that in such a case we can set for V(t) any vector from W 1 ). 



3 Here span(g) means the span of the columns of g. Conv(A) stays for convexification of 
the set A. 
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3.2 The Projection onto H. 

Let L 2 (i?)t the subspace of L 2 (i?) defined by 

L 2 (R) t := {v e L 2 {R) \ v • n = OonT}. 
Let v € L 2 (i?) t . We can write v as 



i'2 



EfcgNxNn W 2fc COS 




Now we want to write w as a sum of an element u £ H (a solenoidal element) 
and a gradient of a function g e L 2 {R) — v = u + Vg. 
We put 



'E^NoxN^sm^) 

v EfceNxN v 2k COs(^)si 



u := ^ u k W k ; q := ^2 9feC0s( 

fceNg feeN 2 \{(o,o)} 




(3.2.1) 

k 2 TTX2 ' 



Then 

Vq 



v E feeN2 -¥?^os(^)si 




For € N \ Nq we find 



-qk^r=v lk , for fceN o x{0} 

b 



-q k h t=v 2k , for fce{0}xN 



and, for k g Nq we obtain 



wife, fc G N x {0} 
b v 2k , k e {0} x No 



ki it 
b 

k 2 TT 



fcl7T k 2 TT 

qk =vik + —r- u k (3.2.2) 

a b 

k2K fclTT 

;— 9fc =W2fc u k . (3.2.3 

o a 

Multiplying (|3~2~2f by and (|3~2~5|l by ^ and then adding the products 

we obtain 

k 2 n k\ir T 

—r-vu v 2k = ku k , 

b a 

obtaining in this way u from v. 
Similarly we can obtain 

/k2^\2 fkiTT.2 fcl7T k 2 7T 

qk - ( — ) qk = — vik + ~rV2k 

v b ' a ' a b 
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or 

T kilT k 2 TT 

kq k = v lk + —r-v 2k - 

a b 



G L 2 (i?) t 



Therefore, given 

\ = /E feeNo xN«i fe sin(^)coste 
^ VE feeN x No ^cos(^)sin(^ 

we can write it, in a unique way, as a sum v = u + Vg, where u e H and 
q € L 2 (R). So the projection map onto the solenoidal (divergence free) space H 

P w : L 2 (i?) t -» ff 

is well defined and we have 

u := = V U^v lk - ^v 2k ) (3.2.4) 

ketq 

q = q\+qi 

Ea /kiirxi\ -r-^ b /k 2 irx 2 \ 

Vik j^ co i—)- E W2 ^ cos (— J 

feeN x{o} fce{o}xN 

(3.2.5) 

Elf fciTT fc 2 7r\ /fci7rai\ (k 2 -nx 2 \ 

k[ Vlk — +V2k -b-) C0S {—) C0S {—)- (3 ' 2 - 6) 

It remains to show that u e H and qi, q 2 & L 2 (R). That follows from the fact 
that both vi and w 2 arc in L 2 : 

feeNg 

< + T < (|| 



fcGN x{0} fce{0}xN o 

a 2 + b 2 , y-^ , 2 a ^ 2 

fcGN o x{0} fee{0}xN o 

<^J^(IKII| 2 + IKII! 2 ); 



El / fci7r fc 2 7r \ ab 



, fc 2 \ a 

fcSN§ \ 



^ E + M)*t < -=L=(|M| L 2 + ||«2 fc || L2 ) s 
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3.3 Spectral Method. 

We want to use Spectral Algorithm, i.e., we want to study NS equation in coor- 
dinates corresponding to the basis of eigenfunctions W defined in 12.3.f|) . 

We consider the equation 

u t + (u ■ V)it + Wp = Ait + F + v; 



and write 



u k W k - 



Then we compute (u ■ V)u = 



u ■ Vmi 

u ■ Vli 2 



Vw 2 = 



£ fce M|M-¥)^COs(^)cOs(^) 

'E feeN ^^(-^)sin(^)sin(to) 



Then 



7i ■ VUl 

n^ n (-— )(-—)— sm^-^-J oos( 



n 2 7rx 2 \ 
6 / 



/mi7ra;i \ 



/m 2 7ra; 2 
C °\—b— 



Eni-n . m 2 7r m 2 7r /?ii7ra;i\ . /n 2 7rx 2 \ 



'rni7rxi\ /m 2 7rx 2 
x sin I sin I 



E7T" 7T 



V 6 
• (ni + mi)irx 



sm( j+sm( 



{ni — mi)TTXi 



cos 



(n 2 + m 2 )7rx 2 



+ cos ^ 



(n 2 ~ m 2 )TTX2 



b J \ b 

(rii + m^TTXi \ . /(n\ — m\)'KX\ 



-)+sin(< 



y~] u m u n -j—-Tnin2~T~n2 sin( 

t-^ _ 4ab b L V a / v a 

/(ra 2 + to 2 )7tx 2 \ /(n 2 - m 2 )7ra 2 
I b J +COS l b 



cos 
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We then arrive to 
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u ■ Vui 

Ett 2 7T. , . f (n\ + m\)-KX\\ f(n 2 +m 2 )irx 2 
u m u„— min2-(m 2 - rc 2 )sin^ J cos(^ 

71-2 i"/ \ /(ni +mi)7rxi\ /(n 2 -m 2 )7rx 2 
+ 2^ u TO u„— mm 2 -(m 2 +n 2 ) sin j cos^ 



tt 7T / (m - mi nxx \ /{n 2 +m 2 ttx 2 \ 
+ 2^ u mW„— min 2 -(m 2 - n 2 )sm^ j cos^ J 

7T 2 7T /(m - mi)7rxi\ / {n 2 - m 2 )nx 2 \ 
+ 2^ u mW„— min 2 -(m 2 + n 2 )sm^ j cos^ J. 

(3.3.1) 





a 


(ni 


+ mi)7TXl 




a 


(ni 


— mi)irxi 




a 


(m 


— mi)irxi 



Similarly: 



u ■ Vw 2 



E, n 2 ir mi7r mi7r . fn\iTXi\ (n 2 Ttx 2 \ 
u m u n (- — )_(-_ _) Bm ^__jco B (— j 



TOl7TXl\ /TO 2 7TX 2 
S1 H— ~ 



+ 



Erii7T TOi7r m 2 7r /ni7rxi\ . /"n 2 irx 2 \ 
u m u n cos sin — - — 

a a o V a / V o / 



/mi7TXi 
X COS I 



7TXl\ / 
— J C0S l 



u m u n -— m\n 2 -m\ 
Aao a 

m,nGNg 



— COS 



m 2 irx 2 \ 
b ) 
(ni + 7711)77x1 



( a ) +C °i 



(m — mi)7TXl ' 



'(n 2 - m 2 )iTX 2 \ 



si H — 1 — ) sm ( — r — ; 



E7T 7T [ / 

«mMn7T m l n 2-1l COS 
4ao a L V 



(ni + mi)7TXi 



) + cos( 



(ni — mi)7TXi 



. /(n 2 +TO 2 )7TX 2 \ / 

sin( j - cos( 



(n 2 — TO 2 )7TX 2 
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and then, 



U ■ V« 2 

4ao a 

E7T 2 7T 
u m u n —— min,2—{mi - m) cos 

+ ^ u m u n ^—rmin 2 — {mi + ni) cosf 
4ao a V 



= ^ u m u n — T min 2 — (ni - mi) cos 



(ni 


+ mi)irxi 




a 


(ni 


+ mi)nxi 




a 


(ni 


— mi)7TXi 



6 

E7T" 7r / (n x - m^KXi \ . f(n 2 - m 2 )irx 2 \ 
u m u„— min 2 -(-mi - nijcos^ J sm^ - J. 

(3.3.2) 



Grouping the terms in sum 1(3. 3.1(1 and 1)3.3.2(1 envolving the product u m u n we 
obtain 



u ■ Vui 

E7T 2 7T, , . /( 

u m u n — rn A n-(m 2 ~ n 2 ) sm^- 



m,nGF 
m<n 



71-2 v/ , \ • f( n i +m 1 )TTX 1 \ ( 
+ u m u n -^mV n-(m 2 + n 2 )smy J cosy 

m<n 

7:2 m ^ / \ ■ ({ni-mi)izxx\ / 
+ " m u„— mv n—(m 2 — n 2 )sm^ J cos^ 

m<n 

71-2 A 77 / \ • /("I _ TOl)7TXl\ / 

+ > u m u n — — -to A n—(m 2 + n 2 ) sm I cos 

^— ' 4ao o V a / V 



("-2 


+ m 2 )nx 2 




b 




- TO 2 )7T22 




b 


("2 


+ m 2 )nx 2 




b 


("2 


- TO 2 )7ra;2 



b 



. 2 7r 1"/ x • /2ni7TXiN 4 

+ U "M nifi2 6™ 1 a J ^ ' 



3.3. SPECTRAL METHOD. 
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and, 

U ■ Vli2 

E7r 2 7r / (ni + mi)«i \ . / 
u m u n to A n— (ni — mi J cos sin 
4ao a V a / V 



m,nSr 
m<n 



7T 2 7T /(m +TOi)7TXi\ . / 

+ > u rn u n m V n— (mi — n\ ) cos I sin 

z — ' „ ' 4ao a V a / V 



m.ngNJ 



m,nGP 
m<n 



("2 


+ m 2 ^ 


)nx 2 




b 




(n 2 


- m 2 


)irx 2 




b 




("2 


+ m 2 


)nx 2 



77 . \ - mi)7rxiN . {n 2 -m 2 )'KX2\ 

+ } u m u n — — -m A n— (— mi — ni) cos sin 

4ab a \ a ) \ o J 



Aab a 

Where the order under the sum sign, in (|3.3.3[1 and l|3.3.4|) . is the lexycographical 
one and, by definition: 

m V n := m\n 2 + ri\m 2 \ m A n := m\n 2 — n\m 2 . 

Put 

7T 2 7T 2 

Cm n = "r~r TO V n & C„ „ = ——rm A n. 

m,n ^ m,n 4flfo 

We rewrite 5 



-(V ■_ 



m<n 

+ UmUnC ra,»l( m2 _ n 2 )sin( 

m.riGNQ 

+ 51 u ro tt n C4 )n r(TO2 +n 2 )sin( 





+ mi \irxi 




a 




+ mi \ttxi 




a 


ni 


— mi \irxi 




a 


|ni 


— mi \irxi 


a 



\n 2 


+ m 2 \nx 2 




b 


\n 2 


- m 2 \nx 2 




b 


\n 2 


+ m 2 \nx 2 




b 


\n 2 


- m 2 \nx 2 



E 2 7T 2 TT /2niTTXl\ 
u nJ-7 n i n 2T(2n 2 )sm[ I; (3.3.5) 



l 4ab b 



4 This last sum cooresponds to the sum over the diagonal {(m, n) S (Nq) 2 | m = n}. 
5 m < n implies mi < m. 
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and 



U ■ Vli 2 

\n 2 + m 2 \irx 2 

[ill — mijcosi I sui) 

' a 

m<n 



= 2^ u mU n C m n -(ni - toi)cos^ jsm^ 

^ i \ (\ni + m x \-KXi\ , (\n2-n 
+ 2^ w m u„C m „-(mi -mjcos^ J sm^ - 

ra2<n 2 

11 , \ f\ n i +m 1 \irx 1 \ . ~ m 2 \irx 2 \ 
+ 2^ Um«nC mi „-(-mi +ni)cos^ jsm^ - J 

m<n 

m2>H2 

/-.V \ fl n l ~ rn 1 \nxi\ . (\n 2 +m 2 \irx 2 
+ 2^ u m u n G m n -(mi+ni)cosy J smy - 



m,nGF 
m<n 



n/\ . \ (\ni-mi\nx\\ . ( \n 2 - m 2 \irx 2 \ 
+ 2^ MmUnC m „-(-mi -ni) cos Jsm^ - J 

m<n 

m2<H2 

/ s (\ n i - mx\-KXi\ . (\n 2 - m 2 \irx 2 \ 
+ 2^ u m,u n C^ n -{m 1 +n 1 )cosy J sm^ - J 



m<n 
m2>n 2 



p 2 f2 T/„ > . {2n 2 nx 2 \ . 



Now we project (u • V)w in -ff : 6 First we put 



n( )m := (rii — mi, 112 — m 2 ); n{ — h)rn := (ni — mi, n 2 + m 2 ); 

n(H — )m := (rii + m 1} n 2 — m 2 ); n(++)m := (ni + mi,n 2 + m 2 ); 

(a,6)+:=(|o|, |6|), «,kZ 7 



and then the projection can be written as 



6 See section IO 

7 Z represents the set of integer numbers. 



3.3. SPECTRAL METHOD. 



51 



P v [{u- V)u] 

m<n 

\ ~> u m u nC^ n n , 

+ 7 / I \ x+ ^+-^(n( + -)m) + 

2 (n(H — )m)+ 

m<n 
m 2 <ii2 

\ ^ u m u nCm,n „ 2 ti/ 

+ T = T === TT +~ H"(n(+-)m)+ 

o (n(H — )m)+ 



m 2 >n 2 



+ / / ,v x+ £> -+ ty (n(-+)m) + 

„,2 (" — him + 

m<n 

\ ~> u m u n C m _ , 

+ E 7 7 ^ = £) --^(n(— )m)+ 

m<n 
m 2 <n 2 

+ E (3-3.7) 

, (n( — )m)+ 

m<n 
m 2 >n 2 



where 



_ (n(++)m) 2 f 7r 7r (n(++)m)^7r 7r 

D ++ = T(m 2 - n 2 ) - — (ni - mi) 

(X (X 

1 (n(H — )TO) 2 f 7T7r (n(H — )m)^7r 7r 
T o a a 

2 (n(H — )m)tinr, . (n(H — )m)firir 

D\_ = ^ / ;2 T m 2 + n 2 - — (ni - mi) 

oo a a 

_ (n( — h)m) 2 f 7T7r (n( — h)m)i~7r7T. 

D -+ = — / T m 2 - "2 - — mi + m 

a a 

, fn( )m)j7r7r, . fn( )m)firir, . 

D — = — / T (m 2 + n 2 - ZJ >—!! _ mi _ ni 

a a 

^? M ) m )^ 7r7r / \ M )m)f IT IT . , 

D 2 _ = ; n i + n2 _ >_>i mi + y 

b b a a 



52 CHAPTER 3. CONTROLLABILITY OF GALERKIN APPROXIMATIONS. 
Equation (|3.3.7|1 can be reduced to 



Hence 



P v [(u- V)u] 

Eu m u n C m n _ 
===={n - m)W (n(++)m)+ 

, (n(++)m)+ 

\ ^ u mUnC m n _ 

, m(H — m)+ 

rri2<n2 

\ ^ u m u n C m n _ 

' 1^ 7=7T \ a+ (n ~ m )^("(+-)»0 + 

, (n(H — )m)+ 

m<n 
m2>"2 

v— \ u m u n C mn _ 

m<n 

v— \ u m u n Cm n _ _ NTI . 

2 (n( )m)+ 

m 2 <ri2 

v— \ u m u n C m n _ 

- 1^ , , x v+ (n-m)Wi;„ ( __ )ro)+ . 



m,n 
m<n 
m 2 >ri2 



P v [(w- V)u] 

======r(n - m)W (n(++)m) + 

, (n(++)m)+ 

-===== (n - m)sign(n 2 - m 2 )iy ( „ (+ _ )m) + 
, mM — )m) + 

m<n 

^— \ u m,U n Cm n /_ _ NTI , 

, (n( — h)m) + 

+ Z_v ~ 7 ? \ t=(" ~ ra)sign(n 2 - rn 2 )W {n( — )m)+ 
— , (n )m + 

m<n 
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53 



or, what is the same 
P v [(u- V)tt] 

E u mUnCm,n /_ _ -, TJT 
===={n - m)W (n(++)m)+ 
n (n(++)m) + 

m<n 

^—v u mU n Cm. n ._ _. . . . 

+ 2^ ~ 1 7 \ =( Tt - TO ) S1 g n ("-2-m 2 )Ty (ra(+ , )m) + 

— „ (n(H — m)+ 

m<n 

\ ■* u mU n C^ n _ 

+ ■====(n-m)sign(n 1 -m 1 )W {n( _ +)m) + 

— , (n( — + )m) + 

- = = =(n ~ ^)sign(ni - mi)sign(n 2 - m 2 )W {n( — )m)H 
(n( )m)+ 



m < n 



(3.3.8) 



3.4 The Infinite ODE System. 



We are now able to write an infinite system of ODEs related with the N-S 
Equation (|1.0.1|) 

u t + (u ■ V)w + Wp = vAu + F + v. 

As we can seen in the strong formulation of the N-S Problem fProblem l2.10.lfl . 
at each time r the last equality is an equality between elements of V' . However 
if we project this equality onto the solenoidal space H we obtain 

P v (ut + (u ■ V)u + Vp) = P v (vAu + F + 

and, by i|2.4.7|) and the fact that r € D(A) a.e., we obtain 

u t (r) = -P v Bu(t) + vAu{t) +F + v(t). (3.4.1) 

We suppose F is solenoidal, otherwise we just take its solenoidal part. 
From 

u = ^ u kW k 



we derive 



Ut 



^2 (u k )tWk 



Au = kUk 



m<n 
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Therefore we may write l(3.4.1|) in the form of the infinite system of ODEs: 



U k = 



= — -(n-m) 



ra,nEN; 
(n(++)m)+=fc 

= — (n — m)sign(ni — mi)sign(n2 — 7712) 

k 

m,nEff 

)m)+=fe 

\ ~* H m M n O m „ _ 

+ y, f (n — mjsignni — mi 

m<n 

(ri( — h)m) + = fc 

\ ~< u m u nC m n _ , 

+ 2^ S — (n — TO)sign(n 2 - m 2 ) 

m,n6N* 
m<n 

(n(H — )m) + =fc 

+ vku k + F k + V k . (3.4.2) 



3.4.1 Galerkin Approximations and G-Saturating Sets. 

Trying to make the writting simpler we introduce some notation referring to the 
coeficients appearing in system <|3.4.2|) : 

Definition 3.4.1. 



C++ . = n " t/wt (n-fh) 

Aab ( m ( ++ ) n )+ 

C;;„ := —====(« - m)sign(ni - mi)sign(n 2 - m 2 ) 
Aab ( TO ( ) n )+ 

4a& (tti( — h)n)+ 

Cm,™ : = 7T==== « - m Sign n 2 - to 2 . 
4a6 ( TO (-| — ) n )+ 

Let 5 be a set of (coordinates of) r linearly independent vectors in H. Write 
also g for the matrix whose columns are those r vectors. 

Rewriting system (|3.4.2|) with the directions in g as the controlled ones, we 
obtain 

u = f(u)+gv v e W 

or, equivalcntly 

u k = (f{u)) k + {gv) k v e R r , keN 2 
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where 

(/(«))* : = J! U mUnC^ + n 



(n(++)m)+=fe 



u m u n C n 



(n( )m)^ — k 



(n( — h)m) +=fc 



m.nENJ 

(n(H — )m)+=fe 



B k (u, u) + vku k + Fk, 



defining B k (u, u) := (f(u))k — iykuk + Fk)- Let us apply a FCE Procedure 
to this infinite system. Applying factorization to this system we obtain, 8 the 
factorized system 

{hx{u)) k = (fx(u))k + gv 1 = (f(u))k+B k (u,gX) + Bk{gX,u) 

+ iy~k(gX) k + B k ((gX) 1 (gX))+gv 1 (3.4.3) 

Now we put 

(V {gX)) k := (fiu^k+gv 1 ; 

{Vi{gX)) k := B k {u, g X) + B k {gX,u) + isk{gX) k ; (3.4.4) 
(V 2 (gX)) k := B k ((gX), (gX)); 

and we note that Vo, Vi, and V 2 are respectively, independent, linear and bilin- 
ear on gX vector fields. 

Now, ginen X £ W we have 

V (g(-X)) = Vo(gX); Vi(g(-X)) = -Vi(gX) and V 2 (g(-X)) = V 2 {gX). 

So, 

/(«) + B(gX, gX) = ± (f x (u) + /_*(«)) £ Conv{f x (u) | X £ W}. 
The set 

{B(gX, gX) \ X £ W} = {B(Y, Y) \ Y £ span(g)} 

is a cone. 

8 See scction [ij.ll 
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Next, from 

Conv({B(gX, gX) \ X e R r } + G = {B(Y, Y) \ Y € span(g)}) + G, 
where G := span(g), we extract the subspace 

G 1 := (<3 + Conv{B(Y, Y) | Y e G}) n (g ~ Conv{B(Y, Y) \ Y e G} 

Now we present two definitions: 

Definition 3.4.2. A finite set of vectors g C D(A) 9 is said Saturating for 

system (|3.4.2f) if the sequence (G J )„ e N of subspaces of H defined recursively by 

G° := G = span(g); 
G j+1 := (G j + Conv{B(Y, Y) \ Y S G J '}) 

n ((P - Conv{B(Y, Y) \ Y € G 3 }) n D(A); 10 

satisfies 



[Jg* = h. 

We say that g is G-Saturating for system (|3.4.2fl if span(g) = span{Wk \ 
k £ JC 1 } for some finite set of modes K. 1 € J-V(Nq), and there exists a sequence 
of subspaces £F such that H° := span(g) and 

H j+1 C (w + Conv{B{Y, Y) | Y e H j }) n (w - Conv{B(Y, Y) \ Y G H j }) ; 



|J = H; 



and, besides there exists a finite subset K,j +2 € J-T^Nq) suc/i that 
span{Wfc | fc e /C J + 2 }. 

Remark 7. in i/ie previuos definition the condition 



(A) \jDi=H 

iGN 

is equivalent to 

(B) \/x E H \j —+ oo only if \x — Hjx\ — > 0] 

where D l stays for either G % or H l . Indeed if (B) is satisfied and x £ H, we 
have HiX € D l and so, (11^2;)^^ is a sequence in lJ ieN -D l converging to x. Thus 

Conversely if (A) is satisfied and x S H , we can find jo G N such that 
\x° — G J0 | < f and, for each n G No we can find j n € N and x n £ D ln such that 
in > jn-x and \x-x n \ < then also \x-Uj n x\ < 11 So (\x-Uj n x\) nS ® 
(that is a subsequence of (\x — Tljx\) converges to zero. Since the sequence 
(D l )nz?t is increasing (in the inclusion sense), the sequence (\x — TljxDj^ is 
decreasing so, it converges to zero (as does its subsequence (\x — Tlj n x\) n( zfi). 

9 Recall that D(A) C H has been denned in II. I). 51 . 

10 We intersect with D(A) because, it is helpful for the study of controllability issues of the 
infinite system. 

11 Because \x — Iij n x\ = min{\x — y\ \ y G G 3 ™}. 
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Definition 3.4.3. A Galerkin approximation of system (I3.4.2|l is the same 
system with the additional condition 

k,n,m eG e TV(Nl). 

From now we shall look for the existence of a G-Saturating set K} for system 
(|3.4.2[1 . The existence of such a set means that for any given finite set of modes 
O C Ng, we want to observe, there exists p S N such that O C H p . 

For j < p the construction of fP' +1 depends only on vectors of iP C H p . 
Therefore from the /C p+1 -Galerkin approximation with K) (1 < j < p) as set 
of controlled modes, i.e., with fP _1 as set of controlled directions, we can arrive 
by a FCE Procedure to the same Galerkin approximation with K? +1 as set of 
controlled modes. 

Iterating FCE procedures, after p steps (starting with K. 1 as set of controlled 
modes), we arrive to the /C p+1 -Galcrkin approximation with JC P+1 as set of 
controlled modes and, the approximate controllability at time t of such system 
is an immediate consequence of Corollarv l3.1.3l The approximate controllability 
at time t of the K. p+1 -Galerkin approximation with K. 1 as set of controlled modes 
follows from the fact that a FCE Procedure does not change closure of attainable 
set at time t. 

Hence the controlling of the modes in K} in the /C p+1 -Galerkin approximation 
we "observe" approximate controllability in the state space H p — span{JC p+1 ) 
and so also in the "O-space" . 



3.5 Looking for Saturating Sets. 

We still working in the rectangle [0, a] x [0, b]. We set two elements n, m E Nq 
such that m < n. We have that 



n — m = 



l> 2 



and, in this equation mi = n± rri2 = So since m < n we have that 

n — fh = implies mi ni & 1112 7^ ni. Then 

a 2 m\-n\ (m - m{]{n x + mi) , , 

n-m = <^ — = j j = ■ -. 3.5.1) 

b z mt,-nf, (m 2 - n 2 )(m 2 + n 2 ) 

Now since mi < n\ which, from H3.5.1|l . implies m 2 > n 2 we obtain 

n — fh — & m < n => n\ > mi & 712 < m,2. (3.5.2) 

By (l3~5~2l we obtain 

Corollary 3.5.1. Under the condition m < n A (mi = ni V ri2 > m2) we have 

C++ = ^ m A n = ^ C~ = 

and 

C — *" C'^ — 
m,n / / m,n" 
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Proposition 3.5.2. The set JC 1 := {(m, n 2 ) G N„ | m, n 2 < 3} \ {(3,3)} is 
G-Saturating for system (I3.4.2[l . 12 

Proof. For any j G No put 

/C J := {(m, 7i 2 ) G Nq | rn, n 2 < j + 2} \ {(j +2J + 2)}. 

We shall prove that the directions in H^ 1 := spa^/O 7 ) can be obtained by a 
suitable FCE Procedure from the directions in H^ 2 (j > 2). 

For every pair (n, to) G JC 1 , n ^ to and for every A G M + we define the 
vectors v^ m , w m,n G ^ r defined by 

( v m,n)n = A; 
( v m,n)m = lj 

. ( w m,Jfc = 0, fe G /C 1 \ {n, m}; 

(>m : Jn = A; 
( W m,n)"i = — lj 

( w m,n)fc = °, & G /C 1 \ {n, m}. 



and 



The vector fields 



and, 



1 2 1 = ' W + 7m '" 

: n ! = /W _ A 7m,« 



where 



7m, n — C mn 9(„(__) m )+ + C m J,9(„( — h)m)+ 
+ ^n,ti^(n(+-)m)+ + ^m,n^K++]m)+i 

belong to Conv{f± v > k n , f± w ^ n I A € K + , to, n E JC 1 , to}. 

Now we shall extract from {/x(t) I G ^ r } a family as follows: 

First we define for each pair (m, ra) G /C 1 , n ^ m the vector <5 min in R N by: 

$m,n — C"m,n e (n( )m)+ "F Crai e (n(-+)ra) + 

+ Cm,n e (nH — )m) + + ^mli e (ii(++)ra)+ (3.5.3) 

Now we select the subfamily Fs 1 '■= {5 m ,n | (to, n) € S\ C (if 1 ) 2 } of this 
vectors, where 

S X = {((1,2), (2,1)); ((1,1), (2, 3)); ((1,2), (2,2)); 

((1, 1), (3, 2)); ((2, 1), (2, 2)); ((1, 1), (1, 3)); ((1, 1), (3, 1))}. (3.5.4) 



12 We take (3, 3) off only because a writting reason. We shall apply a induction procedure 
in the proof and without these "corners" it becomes simpler. 
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The vectors of this family are precisely 

_9n 2 (b 2 -a 2 ) 15^ 2 (a 2 -6 2 ) 
6 M)W ~ 4a6(a 2 + 6 2 ) ^ + 4ab(9a* + & 2 ) 6(1 ' 3) 
157T 2 (a 2 ~b 2 ) ; 7r 2 (& 2 -a 2 ) 

+ 4a6(a 2 + 96 2 ) e(3,1) + 4a6(a 2 + &2) e (3,3) 
_ 7r 2 (36 2 + 8a 2 ) 57r 2 (8a 2 + 36 2 ) 

d(1 ' 1} ' (2 ' 3) " 4a6(4a 2 + 6 2 ) g(1 ' 2) " 4afe(16a 2 + 6 2 ) 6(1 < 4) 
57r 2 (8a 2 +36 2 ) 7r 2 (36 2 + 8a 2 ) 

+ 4a6(4a 2 + 9b 2 ) C(3 ' 2) 4a6(16a 2 + 96 2 ) 6(3 ' 4) 
96tt 2 36tt 2 
d(1 ' 2) ' (2 ' 2) " ~2a(16a 2 + 6 2 ) e(M) + 2a(16a 2 + 96 2 ) 6(3 < 4) 

^ 2 (86 2 + 3a 2 ) 5^ 2 (3a 2 + 86 2 ) 

d(1 ' 1) ' (3 ' 2) " ~4a6( a 2 + 46 2 ) e(2 ' 1) " 4a&(9a 2 + 46 2 ) e(2 ' 3) 
57r 2 (3a 2 +86 2 ) 7r 2 (86 2 + 3a 2 ) 

+ 4a6(a 2 + 166 2 ) S(4 ' 1} + 4afe(9a 2 + 16fe 2 ) 6(4 ' 3) 
9a?r 2 3a?r 2 
d(2 ' 1) ' (2 ' 2) " "2&(166 2 + a 2 ) e(4a) " 26(166 2 + 9a 2 ) 6(4 ' 3) 
2a7r 2 a7r 2 

"(1,1), (1,3) = 5(62 + a2) e (2, 2 ) - &(fe2 + 4fl2) e ( 2 ,4) 

2&7T 2 6tt 2 
^^■C^ - ~^T^) 6(2 ' 2) + a(a 2 +46 2 ) e ^' 

Projecting the vectors in this subfamily on the space span{ek | k € JC 2 \ K. 1 } 
we obtain 

7T 2 (b 2 -a 2 ) 

Ul5 MM> = Aab(a 2 + b 2 ) e ^ (3 - 5 - 5) 

_ . _ 5^ 2 (8a 2 + 36 2 ) ^ 2 (36 2 + 8a 2 ) 

n^(i,i) )(2 ,3) - - 4a6(16fl2 + 62) e (M) - 4a6(16a 2 + 96 2 ) e(M) 

_ . 9&7T 2 36tt 2 

niS {1>2Um - 2fl(16fl2 + 62) e ( i,4) + 2a(16a2 + %2) e (3,4) 

_ 57r 2 (3a 2 +86 2 ) 7r 2 (86 2 + 3a 2 ) 

nid (bD,(3, 2 ) - 4a6(a 2 + 16& 2 ) e(4 ' 1} + 4 a fe(9a 2 + 166 2 ) 6(4 ' 3) 

9a?r 2 3a?r 2 
n^( 2 ,i),( 2l2 ) - - 26(16&2 + a2) e (4,D - 2 fe(16fe 2 + 9a 2 ) e(4 ' 3) 



(l,i),(i,3)-- 6(62 + 4a2) e ( 2,4) 
6tt 2 

IIl ^ 1 )'( 3 - 1 ) = a (a 2 + 46 2 ) e ^- 

Now we consider the case a / d, i. e., if our rectangle R is not a square 
(the case of the square will be considered below in Remark EJ. If a ^ b these 
projections are linearly independent so, by Lemma |3. 1 . II the 15 vectors of the 
family {efc | k G K. 1 }UFs 1 are linearly independent in span{K, 2 ) — H 1 and then 
they span span{ek | k € /C 2 } C H . 
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Now we extract the linear space 

H 1 := (h° + Conv{XS m , n | (to, n) 6 Si, A 6 K}) 

n (h° - Conv{\6 m , n | (to, n) £ Si, A e R}j 

from the space 

'h° + Conv{B(Y, Y) \ Y G H }) n (h° - Conv{B{Y, Y) | Y G ff }' 



Since Conu{A(5 m . n | (to, n) G Si, A G M} ) coincides with span{A(5 mj „ | (to, n) G 
Si, A e 1}], we have that 



H 1 = span(fC U {5 m>n | (to, n) G Si}) = span(JC ) = H . 
The following proposition completes the proof. □ 

Proposition 3.5.3. From the directions in W = span(/C : ' +1 ), (j > l) 13 we 
can obtain the directions in iF + = span{K.^ +2 ) by a FCE Procedure. 

Proof. We consider two cases "j even" and "j odd" . 
• j even: In this case 

W +1 := {(m, n 2 ) G N 2 , I m, n 2 < j + 3} \ {(j + 3, j + 3)}. 

can be written as 

W +1 := {(m, n 2 ) G N 2 \ n u n 2 <2 P +l}\ {(2p +l,2p+ 1)}, 

setting p = Then p > 2. 

As we did before in the case "j = — > j = 1", we extract a subfamily -Fs J+1 := 
{<5 m , n | (m, n) G S J+ i C (/P+ 1 ) 2 } of {5(Y, Y) | Y G fP} where now the 
"selection" is 

S j+1 = {((l,2),(2p,2p-l))} 

U{((l,l),(2z,2p+1)) | z = l, ...,p}U{((l,p+l)(2,p+l))} 
U{((l,l),(2p+l,2z)) | s=l,...,p}U{((p+l,l),(p+l,2))} 
U{((s,l),(s,2p+1)) | s = l,...,p} 
U{((l,*),(2j>+1,«)) | s = l,...,p}. 

If we write explicitely the vectors of Fs j+1 we obtain quite long expressions, for 
example we have that , 2p 2p _i) equals 

(a 2 (-3 + 2p) + 6 2 (-l + 2p)){ix + 2 P 7r) 2 sign(-3 + 2p)sign(-l + 2p) _ 
4a6(6 2 |l-2p| 2 + a 2 |3-2p| 2 ) ' 



J The case j = was already seen in the proof of proposition 13.5,21 
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so, here we will not write those vectores explicitely. Anyway, those vectors are 

tf(l,2),(2p,2p-l) = =C '(l,2),(2p,2p-l) e (2p-l,2p-3) + C(l^),(2p,2p-l) e (2p-l,2p+l) 
+ C '(l,2),(2p,2p-l) e (2p+l,2p-3) + C '( t l^),(2p,2p-l) e (2p+l,2p+l); 

<5(ia),(2z,2p+l) =C, (ia),(2 Z: 2p+l) e (22:-l,2p) + ,(2z,2p+l) e (2z- l,2(p+l)) 

+ C '(ia),(2z,2p+l) e (22:+l,2p) + ^l),(2z,2p+l) e (2z+l,2(p+l)) 

z = l, ...,p; 

<5(l,p+l),(2,p+l) =C(l^+l),(2,p+l) e (l,2(p+l)) + C '(l^+l),(2,p+l) e (3,2(p+l)); 
<5(l,l),(2p+l,2z) =C, (ia),(2p+l,2 Z ) e (2p,2^-l) + C(M) : (2p+l,2 Z ) e (2p,2* + l) 

+ C '(l,l),(2p+l,2 Z ) e (2(p+l),2^-l) + C(^),(2p+l,2 Z ) e (2(p+l):2^ + l) 

z = 1, . .. , p; 

<5(p+l,l),(p+l,2) =C, (p+l,l),(p+l,2) e (2(p+l),l) + C '(ptl,l),(P+1.2) e (2(P+l).3); 
<5(s,l),(s,2p+l) =C(! 5 ,l),( s ,2p+l) e (2s,2p) + C(!s4),( S: 2p+l) e (2s,2(p+l)) 

s = 1, . .. , p; 

5(l,s),(2p+l,s) =C(M),(2p+l, s ) e (2p,2 S ) + C( l l^),(2p+l :S ) e (2(p+l),2 S ) 

s = 1, . . . , p; 

And, projecting them onto the space span{ek | k G /C J+2 \ /C J+1 } we arrive to 
the family Hj + iFs j+1 whose elements are 



n.f+l<5(l,2),(2p,2p-l) -C(t^2) : (2p,2p-l) e (2p+l,2p+l); 

nj+l5(l,l) 1 (2 Z) 2p+l) =C, (M),(2 Z ,2p+l) e (2z-l,2(p+l)) + C '(14),(2 Z ,2p+l) e (2z + l,2(p+l)) 

Z = l, ...,p\ (3.5.6) 

n j + l (5 (l,P+l),(2,p+l) =C, (l^ + l).(2,p+l) e (l,2(p+l)) + C, (ti+l).(2,p+l) e (3,2(p+l)); 

(3.5.7) 

nj+l^(l,l) I (2p+l,2z) =Cci,l),(2p+l,2 2 ) e (2(p+l),2 Z -l) + G (l^l),(2p+l,2z) e (2(p+l),2z+l) 

z = l, ...,p; (3.5.8) 

n j+l^(p+l,l),(p+l,2) = C '(p+l,l),(p + l,2) e (2(p+l),l) + C '(ptl,l),(p+l,2) e (2(p+l),3); 

(3.5.9) 

nj+i5( a) i),( S) 2p+i) =C(^),(^2p+i) e (2s,2(p+i)), s = 1, . . . , p; (3.5.10) 

nj+l5(l, s ),(2p+l, s ) =C(^),(2p+l, S ) e (2(p+l),2 S ) J S = 1, . . . , p. (3.5.11) 



No one of the cocficients appearing in these expressions vanishes because all 
pairs (m,n) satisfy m < n A (mi — m V n 2 > m 2 ) and because no one of the 
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following expressions vanish 

(l,2)A(2p,2p-l) = -l-2p; 
(1, 1) A (2z, 2p + 1) = 1 + 2p - 2z, z = 1, . . . , p; 
(l,p+l)A(2,p+l) = -(l+p); 

(1, 1) A (2p + 1, 2z) = -1 - 2p + 2z, z = 1, . . . , p; 
(p + l,l)A(p + l,2) = l+p; 

(s,l)A(s,2p+l) = 2ps, s = l, ...,p; 
(l,s)A(2p+l,s) = -2ps, s=l, ...,p. 

Hence we can see that these vectors are linearly independent. Indeed it suffices 
to prove that: 

• The vectors n i+ i5 (1)1 ) > (2 i 2 P +i), (z = Ihi (|3.5.6[l) andn 3 - + i5( liP+1 ) i ( 2j p + x) 
are linearly independent; and 



• The vectors n j ' + i5( 1)1 ) > (2 P+ i i 2), (z = lin l|3.5.8|l) andn j - +1 i( p+1)1 ) i ( p+li2 ) 
are linearly independent; 

But that comes from 

C (M),(2,2 P +i) = (3 + 2p)(36 2 + 4o 2 p(l+p)) 
C iC P+ m^ P+ i) % 2 (l+p) 
. (-l + 2p)(36 2 + 4a 2 p(l+p)) = Cl ("iq),(2,2p+i) 

>(2,P+1) 

and 

C ( + iri),(2 P +i.2 Z ) = (3 + 2p)(3a 2 + 46 2 p(l+p)) 
^MMp+M) ~ 9a 2 (l+p) 

(-l + 2p)(3a 2 + 46 2 p(l+p)) _ g(M),(2,2 P+ i) 
3a 2 (l+p) C '(i^+i),(2,p+i) 

By Lemma ED"]] the (2(p + l)) 2 - 1 of F Sj+1 U {e fe | fee /C J+1 } are linearly 
independent and then they span 

span{e k \ 1 < fc x , fc 2 < 2(p + 1)} \ {(2(p + 1), 2(p + 1))} 
= span{ek | A; € /C J+2 }. 

Since with i5 mjn also A5 m .„ belongs to {-B(Y, F) | 7 6 iP} for every A £ K, wc 
can select the linear space W +l defined by 



H j+1 := (h j + Conv{XS mj7l | (to, ra) e A € I 

n (iF - Con<u{A5 ro ,„ | (m, n) e S j+ i, A e M}) 

from the space 

(W + Conv{B(Y, Y) | Y G fP'}) n (iP - C*onu{B(r, F) | Y 6 7P}Y 
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Since Conv{\5 m . n | (to, n) G Sj+i, A G M}^) coincides with span{A<5 mi „ 
(to, n) G Sj+i, A G M}j , we have that 

W +1 = span{1C 3+1 U {<5 m ,„ | (to, n) G S^+i}) = span{1C 3+2 ) = H j+1 . 



Now we study the other case 
• j odd: In this case 



/C7+ 1 := {(n 1; n 2 ) G N 2 , | m, n 2 < j + 3} \ {(j +3,j + 3)}. 



can be written as 



V +1 := {(m, n 2 ) G N§ | m, n 2 < 2p} \ {(2p,2p)}, 
setting p = i±2. Then p > 2. 

We extract the subfamily i*s,, +1 := {<5 m ,„ | (to, n) G Sj+i C (K' J+l ) 2 } of 
{B(F, y) | Y eW] where now the "selection" is 



S j+1 = {((l,2p-l),(2p-l,l))} 

U{((l,l),(2z-l,2p)) | z = 2, ...,p}U{((l,p)(3,p+l))} 
U {((1, 1), (2p, 2z - 1)) | z = 2, . . . , p} U {((p, 1), (p + 1, 3))} 
U{((l,l),(2«,2p)) | a = l, ...,p-l}U{((l,p),(2,p+l))} 
U {((1, 1), (2p, 2s)) | s = 1, . . . , p - 1} U {((p, 1), (p + 1, 2))}. 

Those vectors are 

£(l,2p-l),(2p-l,l) =C(i, 2 p-i),(2p-i,i) e ( 2 (p- 1 ). 2 (p- 1 )) 

+ C '(1^2p-l),(2p-l,l) e ( 2 (P- 1 )' 2 P) 

+ C '(l,2p-l),(2p-l,l) e ( 2 J'' 2 (P- 1 )) + C '(l^p-l),(2p-l,l) e ( 2 P^ 2 P) 
<5(l,l),(2z-l,2p) =C, (l7l),(2 Z -l,2p) e (2( 2 -l),2p-l) + C '(M),(2 Z -I,2p) e (2( 2 -I),2p+1) 
+ C '(l,l),(2 Z -l,2p) e (2z,2p-l) + C '( t ia),(2z-l,2p) e (22;,2p+l) 

z = 2, . .. , p 

<5(l,p)(3,p+l) =C, (l,p)(3,p+l) e (2,l) + C(i^,)(3,p+l) e (2,2p+l) 
+ C, (l,p)(3,p+l) e (4,l) + C '( , l^)(3,p+l) e (4,2p+l) 
<5(l,l),(2p,2z-l) =C( M ) : (2p,2 Z -l) e (2p-l,2(*-l)) + C (1 + } (2p i2z _ 1) e(2p_l,2z) 
+ C, (l,l):(2p,2 2 -l) e (2p+l,2(2:-l)) + G (M),(2p,2z-l) e (2p+l,2z) 

z = 2, . . . , p 
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■(p,l),(p+l,3) = C (p~l),(p+l,3) e (l,2) + C '(p,l),(p+l,3) e (1.4) 

+C (p,l),(p+l,3) e (2p+l,2) + C '(pa),(p+l,3) e (2p+l,4) 

5(l,l),(2a,2p) =C(l,l),( 2iS ,2p) e (2 S -l,2p-l) + C(14),(2 s ,2p) e (2s-l,2p+l) 
+C, (l,l),(2s,2p) e (2s+l,2p-l) + C(ta),(2^2p) e (2«+l,2p+l) 
S = 1, . .. , p- 1 

; (l,p),(2,p+l) =C (i^),(2,p+l) e (l,l) + C '(li),(2,p+l) e (l,2p+l) 
+C, (l,p),(2,p+l) e (3,l) + C '(li),(2,p+l) e (3,2p+l) 

5(l,l),(2p,2s) =C? (M),(2p,2 5 ) e (2p-l,2s-l) + C '( : L4),(2p :2s ) e (2p-l,2 S +l) 

+ C '(l,l),(2p,2 S ) e (2p+l,2 ;S -l) + C (l^l),(2p,2 S ) e (2p+l, 2 H-l) 
S = 1, . .. , p- 1 

(p,l),(p+l,2) =C (p4),(p + l,2) e (l>l) +C, (pa),(p+l,2) e (l>3) 

+C, (p,l),(p+l,2) e (2p+l,l) + C '(p4),(p+l,2) e (2p+l,3); 



And, projecting them onto the space span{ek \ k € JC' J+2 \ /C J+1 } we arrive to 
the family TLj + iFs j+1 which elements are 



n,-+l<5(i,2p-l),(2p-l,l) - C '(^p-l),(2p-l,l) e (2p I 2p) 

nj + l<S(l,l),(2z-l,2p) =C(l^L) : ( 2z -l,2p) e (2(^-l),2p+l) + C(ia),(2^-l,2p) e (2^,2p+l) 

z = 2,...,p (3.5.12) 

n j+l<W)(3,p+l) =C, (l^)(3,p+l) e (2,2p+l) + C '(l^)(3 : p+l) e (4,2p+l) 
nj+l5(l 1 l),(2p,2z-l) =C(l,l),( 2 p,2 2 -l) e (2p+l,2(2:-l)) + C(ia),(2p,22-l) e (2p+l,22) 

z = 2,...,p (3.5.13) 

n j+l<W),(p+l,3) =C, (pa),(p+l,3) e (2p+l,2) + C '(p4),(p+l,3) e (2p+l:4) 
11^+1(5(1,1), ( 2s ,2p) =C(l^L),(2 s ,2p) e (2s-l,2p+l) + C(t4),(2 s ,2p) e (2s+l,2p+l) 

s = l,...,p-l (3.5.14) 

nj + l<5(i, p) ,(2,p + l) =C(i^),(2,p + l) e (l,2p+l) + C '(l^),(2,p+l) e (3,2p+l) 
nf+l(5(i,i),(2p,2 5 ) = C '(iri),( 2 p,2 iS ) e (2p+l,2 S -l) + C '(^4),(2p,2 S ) e (2p+l,2 S +l) 

s = l, ...,p-l (3.5.15) 

nj+l^(p,l),(p+l,2) =C(pri),(p+l,2) e (2p+l,l) + C '(p4),(p+l,2) e (2p+l,3); (3.5.16) 



No one of the coeficients appearing in these expressions vanishes because all 
pairs (m, n) satisfy m < n A (mi = m V n 2 > m 2 ) and because no one of the 
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following expressions vanish 



(l,2p-l)A(2p-l,l) = l-(2p-l) 2 ; 

(1, 1) A (2z - 1, 2p) = 1 + 2(p - z), 2 = 2,..., p; 

(l,p)A(3,p+l) = l-2p; 
(1, 1) A (2p, 2z - 1) = -1 - 2(p - «), z = 2, . . . , p; 
(p,l)A(p+l,3) = 2p-l; 

(l,l)A(2 S ,2p) = 2(p- S ), s=l,...,p-l; 
(l,p)A(2,p+l) = l-p; 

(1, 1) A (2p, 2s) = 2(s -p), s = 1, . . . , p - 1 
(p,l) A (p+1,2) = p- 1. 



Hence we can see that these vectors are linearly independent. To see this is 
enough to see that: 



• The vectors U J+ iS^^^ 3 ^ p) , (z = 2 in i(3.5.12(ll and ilj+ifyi^^p+i) are 
linearly independent; 



• The vectors IIj + x^(i,i),(2p,3) ! (z = 2 in l|3.5.13|l) and Ilf+i J^x),^^^) 
are linearly independent; 



• The vectors Hj+iS(i,i),{2,2 P )i (s = 1 in 13.5.14(1 1 and Hj+iS^^^^+i) are 
linearly independent; 



• The vectors Ilj+i^i,!)^^), (s = 1 in 1(3.5.15(1^1 and II,- + i£(p,i),(p+i,2) are 
linearly independent; 



But that comes from 

C (Ci).(3,2 P ) _ (3 + 2p)(8fe 2 + a 2 (4p 2 - 1)) 
^)(3,p + i) (l + 4p)(86 2 + a2(2p+l)) 



(2p - 3)(86 2 + a 2 (4p 2 - 1)) _ C t%X^ P ) 



(2p-l)(86 2 + a 2 (2p+l)) C++ )( 3 !p+1) 



C (M),(2 P ,3) _ (3 + 2p)(8a 2 + 6 2 (4p 2 - 1)) 



^(ti).(P + i,3) (l + 4 P )(8a 2 + 6 2 (2p+l)) 
(2p - 3)(8a 2 + 6 2 (4p 2 - 1)) _ ^( + i + i),(2 P ,3) 



(2p-l)(8 a 2 + & 2 (2p+l)) C++ (p+li3) 
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C (i\),(2.2 P ) _ 2(1 + P)(36 2 + a 2 (4p 2 - 1)) 



2(36 2 + a 2 (V-l)) _ C++ ) (22p) 



3& 2 + a 2 (2 P+ l) C++ )(2;p+1) ' 
and 

g (M),(2 P ,2) = 2(l+p)(3a 2 + 6 2 (V-l)) 
C (p:i),(p + i,2) (3p+l)(3a 2 + 6 2 (2p+l)) 
2(3 Q 2 + b 2 (4p 2 -l)) _ g( + ia),( 2p ,2) 



3& 2 + a 2 (2 P+ l) C ( ; i)(p+1 , 2) 

By Lemma ETU the (2p + l) 2 - 1 of F Sj+1 U {e fe | fc G /C 3+1 } are linearly 
independent and then they span 

span{e k | 1 < ki, k 2 < 2p + 1} \ {{2p + l,2p+ 1)} = span{e k \ k G K, 3+2 }. 

Again, since with S m , n also A<5 mi „ belongs to {B(Y, Y) \ Y E H J } for every 
A £ R, we can select the linear space H 3+1 defined by 

H j+1 := fiF + Conv{\5 m , n | (m, n) G Sj+i, A € 

n (iF - Conv{\5 m , n | (to, ra) G Sj+i, A € I 

from the space 

'iF + Coni;{B(y, F) | F G iP}) n f/f j - Conv{B(Y, Y) \ Y G i/ 3 } 



Since Conw{A<5 mjn | (to, n) G Sj+i, A € R} ) coincides with span{A^ 



(to, n)6S;+i,A6R}), 



we have that 



W +1 = span{JC J+1 U {c5 m ,„ | (to, n) G S J+X }) = span(K. j+2 ) = W +1 . 

□ 

Remark 8. Starting with span(IC 2 ), by Proposition \3. 5.7ft we can obtain the 
directions on span{K, 3 ), (j > 2) iterating FCE procedures. Since that Propo- 
sition is valid for any rectangle (including the square) we can say that K- 2 is 
G-Saturating for all rectangles. 

Remark 9. We can show that K, 1 is saturating for the square too and so com- 
plete the proof of Provosition TB. 5. 'ft The only step that is not clear in the square, 
is how to arrive to span{JC 2 ) (or to a superspace of its) from spaniJC 1 ). In the 
case of the square the extracted family in the end of proof of Proposition Vd.d.'lA 
( Case j even: j = in our present case) does not lead to a family of linearly in- 
dependent vectors when projected in span{ek | k (fc IC 1 } because 115(1^)2,1) = 0. 
Then we proceed as follows: First we select the family Fg 1 \ {<5(i.2),(2,i)}- Pro- 
jecting this family onto span{K 2 \ ({(3, 3)} U/C 1 )) we obtain a family of linearly 
independent vectors. So "adding" the vectors in {e^ | k € IC 1 } to those Fs 1 we 
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obtain a family of 4 2 — 1 — 1 linearly independent vectors spanning the space 
span{ek | k G JC 2 \ {(3,3)}}. So proceeding as before we reach the directions 
span(K 2 \ {(3, 3)}). Then we select the family 

. _3ir 2 (5a 2 +b 2 ) 9n 2 (5a 2 + b 2 ) 

"(1,1), (2, 4) - 2a5(9a 2 + 52) e (i.3) - 2ab(25a 2 -b 2 f^ 

ir 2 (5a 2 + b 2 ) 3ir 2 (5a 2 + b 2 ) 



' 2ab(a 2 + b 2 ) e(3 ' 3) ~ 2ab(25a 2 + 9b 2 ) 6(3 ' 5) 
. 7r 2 (5a 2 + 36 2 ) 77r 2 (5a 2 +36 2 ) 

~ - 4ab(a 2 + b 2 ) 6 ^ 4a6(25a 2 + 6 2 ) e(1 ' 5) 
77r 2 (5a 2 + 36 2 ) 7r 2 (5a 2 



4ao(a 2 + 9b 2 )' e{3 ' 1) 4a6(25a 2 + 9b 2 f {3 > 5) 
_2l7r 2 (6 2 -5a 2 ) 277r 2 (5a 2 - b 2 ) 

~ Aab(9a 2 + b 2 ) 6(1 < 3) + 4a6(25a 2 + b 2 ) C(1 < 5) 

37r 2 (5a 2 -b 2 ) 2l7r 2 (b 2 -5a 2 ) 
+ 4a&(a 2 + b 2 ) e(3 ^ 3) + 4a6(25a 2 + 9o 2 ) 6(3 ' 5) ' 

Projecting in span{ek | k G JC 3 \ (JC 2 \ {(3, 3)}) } we obtain 

9n 2 (5a 2 + b 2 ) 
n<5 (M),(2,4)=- 2ab{25a 2_ b 2f(^) 

w 2 (5a 2 +b 2 ) 37r 2 (5a 2 +6 2 ) 
+ 2a6(a 2 + b 2 ) 6(3 ' 3) 2ao(25a 2 + 9b 2 ) 6(3 ' 5) 
7^ 2 (5a 2 + 3o 2 ) ^ 2 (5a 2 + 3o 2 ) 

U6 ^^ - - 4a&(25a 2 + 6 2 ) e(1 ' 5) + 4 fl 6(25a 2 + 96 2 ) e(3 ' 5) 
277r 2 (5a 2 - b 2 " 
' 4a6(25a 2 + b 2 )' 
37r 2 (5a 2 -& 2 ) 21^ 2 (6 2 -5a 2 ) 
+ 4a&(a 2 + b 2 ) 6(3 ' 3) + 4a6(25a 2 + 9& 2 ) 6(3 ' 5) ' 

Since we are working in the square a = b and no one of the coeficients appearing 
in the last expressions vanish. We compute 

(f~i — h si-\ — 

°(1,1),(2,4) °(1,1),(2,4) °(1,1),(2,4) 
^(1,2), (2,3) ^ (1,2), (2,3) 

°(1,4),(2,1) °(1,4),(2,1) °(1,4),(2,1), 

157r 2 (125a 6 + 75a 4 o 2 - 5a 2 o 4 - 326 6 ) 
16a 3 o 3 (a 2 + o 2 )(25a 2 + o 2 )(25a 2 + 96 2 ) 

(since a = b) 

2880a 6 45 



n5 (i.4),(2,i) =7717^2— T2T e (i,5) 



28288a 12 442a 6 



Hence the vectors are linearly independent. If we joint to the three 8 m . n vectors 
the vectors in {e^ \ k G /C 2 \ {(3,3)}} we obtain a family of 4 2 + 1 linearly 
independent vectors spanning the space {ek \ k G K. 2 U {(1, 5), (3, 3), (3, 5)} }. 
We arrive in this way to span(K. 2 U {(1, 5), (3, 3), (3, 5)}) D span(JC 2 ). Since 
JC 2 is G-Saturating for the square so is JC 2 U {(1, 5), (3, 3), (3, 5)}. Hence K. 1 is 
saturating for the square. 
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3.6 Exact Controllability of Galerkin Approxi- 
mations. 

Write the /C^-Galerkin approximation of the NSE system l|3. 4.2(1 . with JC 1 as 
set of excited modes, in the concise form 

ilk = B k {u) + vAuk + F k + v k k e K, 1 
V :\u k =B k {u) + vAu k +F k k e K N \ K> (3.6.1) 

u e R KN . 

In ^Hj E. Weinam and J. Mattingly proved the Full Lie Rank Property for the 
2D NSE with periodic conditions and for some class of few low modes controls. 
Similarly we prove that our equation also is "full Lie rank", i.e., Lie brackets at 
each point span the ambient space R KJV . 

Before we have proved that for all N £ No and all t > the system [ 1)3. 6.10 . N] 
is time-i approximately controllable: 



Vu£R KN A u (T N ){t) =R KN 

where Tm is the family of vector fields of system [ 1)3.6. lfl .N]. i.e., 

T N = {£(•) + i/A(-) + F KN + v | v£ R K1 }. 

Next we prove the (exact) controllability of system [ i)3.6.1)) .N]. For that we need 
to compute some Lie brackets. 

3.6.1 Lie Brackets. Full Lie Rank Property. 

From Tni we set the vector fields 

V° :=B + vA + F KN , X l :=V° + ^- 

oui 

where F KN is the projection of F onto R 1 ^ and km — #IC N and, compute 



so 



fc^(n++i) + k=(n-\ — i) + 

k=(n — H) k— (n i) + 

i<n i<n 



where S^ k '^ is the Kroncckcr delta function: S^ k '^ — 



1 if k = i 
if k ^ i ' 



8V i 

V 1 ' 1 := [X 1 , V 1 } = [X j , \X\ V }} = — 
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so. 

V 3 ' l =^,j, j>i. 
At each given point u G K KJV the family of brackets 14 

H 1 := {V° ± i^-, a 7iJ | n G R"\ (i, j) G Si} 

span a superspace of WL K2 . Indeed, = (V° + - (V° - . 15 

Now we have a technical difficulty, the vector fields have coeficients with 
not so nice expressions and, if we compute Lie brackets envolving them we 
will obtain even more complicated expressions. To avoid these expressions we 
prove by (finite) induction on i G {1, . . . , kn} that each constant vector fied 
■J^—, n G K- 1 is a linear combination of brackets: 



• For i — 1 we take the family H 1 . 

• Inductive Step: The induction hypothesis is: 

"There is a family of brackets W~ x = {W^ \ j = 1, . . . , M p _i} such 
that every constant vector field i G KP~ l can be written as a linear 

combination of its elements: 



d 

dui 



A/p-i 

i=i 



Then for all i G K. p ~ l 
d 



V 1 := 



du,. 



V° 



G span{\Wj, V°] | ] = 1,..,,M H } 



and, for each i, n G K. p 



V" 



r d 

du r , 



, V 



G span{[W k , [Wj, V }} \k,j = l,..., M p _i}. 



Since the vectors in {gf-, 7i, n | G /C p_1 , (i, n) G Sp-i} span R Kp and 
can be written as a linear combination of brackets, then also the vector 
fields -Tp-, i G K. P \K. P ~ 1 are linear combinations of brackets. The wanted 

aui 5 \ 

family is W~ x U {[W fc , [W,-, V ]) \ k, j = 1, . . . , M p _ x }. 



Therefore, for all iV G No, system [ (13.6. It . N] is a full-rank bracket generating 
system. From that and from its approximate controllability 16 we conclude its 
controllability. Unfortunately for fixed time the bracket generating property is 
not suficient to conclude controllability from approximate controllability. To 
achieve controllability at time t we shall need some lemmas which proofs can be 
found in 181. 



14 Including the elements of J 7 ^ we consider brackets of "length" 0. 

15 Recall that Si has been chosen so that {e n , Si 4 | n S K. 1 , (i, j) G Si} span l" 2 . 

16 Approximate controllability at time t implies, trivialy approximate controllability. 
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3.6.2 Zero Orbits and Zero Ideal. 

Definition 3.6.1. A zero-time orbit Nq u throughu of a family of vector fields 
T is the set 

v 

N 0u ■= {uoe tlVl o . . . o e*^ | pe No, Vi G T, t i ^, s £ j t i = 0} 

i=l 

Definition 3.6.2. The derived algebra of T , denoted T> er (T), is the set of 
all linear combinations of iterated brackets 17 

The zero-time ideal, denoted T(T), is the span of elements in T> er (!F) and 
differences of the form X — Y with X and Y in T . 

Lemma 3.6.1. Let T be any family of analytic vector fields on an analytic 
manifold M. Let N be an orbit of T and, Nq be a zero orbit of T contained in 
N . Then we have the following: 

• Each connected component of Nq is an orbit ofT(J-); 

• For each u S Nq, the tangent space of Nq at u is equal to the evaluation 
of T(J-) at u; 

• The dimension of T U (T) is constant as u varies on N. It is equal either 
to dim(Lie u (!F)) — 1 or to dim{Lie u {!F)) ; 

• c?im(Lie li (^ r )) = dim(T u (J-)) iff X(u) S T U (T) for some X G T . 

Lemma 3.6.2. Suppose that T is a family of vector fields on M such that 
both T and its zero-time ideal T (J 7 ) are Lie-determined (the evaluation of Lie 
brackets at each point span the tangent space to the orbit). In addiction, assume 
that T contains a complete vector field. Then 

• A u (T)(t) is a connected subset of some zero orbit Nq z through some ele- 
ment z € M. 

• A u (^F){t) has a nonempty interior in the manifold topology of the zero- 
orbit where it is contained. Moreover, the set of interior points is dense 
in A u {T){t). 

Coming back to our system [ l|3.6.1() .N]. by Lemma \3 . 6 . 1 1 and due to the fact 
that V°(0) = e we have 

dim(Ljeo(^ r v)) = dim(Xo(-? r A r )) = 

which means that the zero-time orbit Nq through has dimension kn and, since 
that dimension is constant in all points in the unique orbit R KJV of the system, 
we conclude that Nq is a union of connected components of dimension kn- Since 
the dimension of that components is kn their topology coincide with that of 
R KJV and, from the fact that the zero-time orbits form a partition of R KJV we 
conclude that M rejv is a union of connected open sets. Therefore there is only 
one zero-orbit, it is the whole state space M. KN . 

17 Brackets of "length" > 1, considering the elements of T brackets of length 0. 
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By Lemma f3.fi. 21 and by the fact that V° is a complete vector field which 
follows from the estimate \u(s)\ < \u(0)\ + ^||F||y-/ (see estimate (|2.6.8[) with 
F = F), the interior intA u (^N)(t) of the attainable set from u at time t is 
dense in A u (!F)(t), where the interior and density are relative to the topology 
of M KJV because that is the topology of the zero-orbit. Hence we arrive to the 
equality 

intA u (T N ){t) = A u {T N ){t) = R KN 

for all t > 0. 

Now we can prove the controllability at time t of system [ (|3.6.1|I .N]: Let u, z 
be two elements in R KJV . Since the intersection of two open dense sets stills open 
and dense, we may take a point 

w € intA u {T N ){t/2) n intA z (-F N )(t/2). 

[Note that the family — Tn '■= { — V | V E J-} satisfies the requirements of 
lemmas I3.fi. II and I3.fi. 21 ] 
Then we can write 



w = u o e tlVl o • • • o e tnVn 



w = zoe~ SlWl o ••■ oe - SmW " 



" t 

V i €F N ,t i >0,J2 t i = p 



So, z is reachable from u in time t: 

z = u o e* 1 Vl o • • • o e 4 "^ o e s ™ Wm o ■■■ o e SlWl . 
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Chapter 4 



Controllability in Observed 
Component. 

4.1 Controlled N-S Problem. Existence, Unique- 
ness and Continuity. 

In chapter we have presented the classical, weak and strong formulations of 
(2XTJI- $H>M- So for the controlled version with F + v in the place of F we 
have existence, unicity and continuity in the data (uq, F + v, v), because our 
control is an essencially bounded function, so that F+v will belong where F did: 
(L 2 (0, T, V) or L 2 (0, T, H)). If we consider the initial data as (u , F, v, v), 
by Theorems 12.6.41 12.8.11 and 12.8.21 we easily conclude that (considering again 
the external force depending on time): 

Corollary 4.1.1. Given 

(u , F, v, v) £ H x L 2 (0, T, V') x L°°(0, T, V')x]0, +oo[, 

there is at least one weak solution u G C([0, T], H) for Problem \2. 5. 1\ with F + v 
in the place of F. 

and, 

Corollary 4.1.2. The maps 

§ : H x L 2 (0, T, V) x L°°(0, T, V')x}0, +oo[ C([0, T], if) 

(u , F, v, v) h-> w 

and 

S 2 : H x L 2 (0, T, V') x L°°(0, T, y') x ]0, +oo[ -> L 2 (0, T, 7) 

(u , F, w, i/) u 

are continuous. 

Similarly, by Theorems 12. 11. II 12. 13. II and 12. 13.21 we can easily see that 
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Corollary 4.1.3. Given 

(n 0) F, v, v) e V x L 2 (0, T, £T) x L°°(0, T, ff)x]0, +oc[, 

there is at least one weak solution u <E C([0, J 1 ], V) /or Problem \2.10.1\ with 
F + v in the place of F. 

and, 

Corollary 4.1.4. TVie maps 

§ s : V" x L 2 (0, T, V) x L°°(0, T, i?)x]0, +oo[ C([0, T], V) 

(u , F, v, v) t-> u 

and 

S 2s : V x L 2 (0, T, V) x L°°(0, T, H)x]0, +oo[^ L 2 (0, T, 
are continuous. 

4.2 Change of Variables. 

If we make the change of variables 

u = y + Iv 

where I is the primitive operator — pfo](i) = f v(t) dr, from 

v! = -vAu - Bu + F + v 

we arrive to the equation 

y = -vA(y + Iv) - B{y + Iv) + F. 

Note that the function v appears only implicitly in the last equation. Now we 
forget that Iv is a primitive of an essentially bounded function and replace it 
by P in the equation. Since v is a low modes forcing it takes value in a finite- 
dimensional space and, It; being a primitive we have Iv <G C([0, T], D(A)). But 
we take P in the larger space L 4 (0, T, D{A)). 



4.3 Weak Case. 

We want to study the following equivalent problems 14.3.11 and 14.3.21 



4.4. EXISTENCE. 
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Problem 4.3.1. Given 

F S L 2 (0, T, V'), P e L 4 (0, T, D(A)) (4.3.1 
& 

y eH (4.3.2 
to find 

yeL 2 (0,T,V) (4.3.3 
satisfying (in the distribution sense) 
Vv e V : 

^(y, «) + + P, «)) + b(y + P, y + P, v) =< F, v >, (4.3.4 
and 

y(0) = yo- (4.3.5 

Problem 4.3.2. Given 

F e L 2 (0, T, V'), , P G L 4 (0, T, (4.3.6 
& 

Vo G ff, (4.3.7 

to /md 

y G i 2 (0, T, V), j/ G L x {0, T, V') (4.3.8 
satisfying 

y' + vA(y + P) + B(y + P)=F on ]Q,T[, (4.3.9 
and 

2/(0) =2/o- (4.3.10 

Remark 10. The equivalence of these problems can be shown the way we proved 

the equivalence of •problems Iff. .5.71 and t^.5.S\ (see chapter\^). This equivalence 
follows from 

• Fe L 2 (0, T, V')\ 

• A(y + P)eL 2 (0, T, V). 
Indeed \\A(y + P)\\ v , < \\y + P\\ and y + P £ L 2 (0, T, V); 

• B(y + P) 6 L x (0, T, y')- 
Indeed ||P(y + P)||y, < C||y + P|| 2 and y + P e L 2 (0, T, 7). 

So tfiai F - + P) - P(y + P) e L x (0, T, V')- 

4.4 Existence. 

We have the Theorem 

Theorem 4.4.1. Given F, P and yo satisfying ()4.3.6J) and i[4.H.7[l . There is 
at least one function y satisfying (|4.3.8|I - (|4.3.1UI) . 
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The proof is analogous to that of Theorem 12.6.41 Basically it differs only in 
some estimates we compute now: Following the proof of Theorem l2 . 6 . 4l presented 
in chapter[21 with the suitable adaptations to problems 14.3. 114.3.21 we define an 
approximate solution 

y m = £ v™W w * 

max{«i, i2}<m 

for each m € No and arrive to the equation 



< (y m )\ y m >+»< A(y m + P m ), y m > 

+ < B(y m + P m ), y m >=< F, y m > . 1 (4.4.1) 

We note that 

. < A(y m + P m ), y m >= \\y m \\ 2 + ((P m , y m )) 

• b(y m + P"\y m + P m ,y m ) 

= b(y m + P m , y m + P m , y m + P m ) - b(y m + P m , y m + P m , P m ) 

= -b(y rn + P n \ y m + P m , P m ) 

= -b(y m + P m , y m , P m ) - b{y m + P m , P m , P m ) 

= -b(y m + P n \ y m , P m ) = ~b(y m , y m , P m ) - b(P m , y m , P m ). 

Hence from 114.4. ID we have 



±\ y m \ 2 + My m \\ 2 = -M(P m ,y m )) 

- b(y m , P m , y m ) + b(P m , y m , P m )+ < F, y m >; 

hence 

±\y^ +2 u\\y m \\ 2 < 2H|P||||y ro ||+C|y ro |||^||||P||+C||P|| 2 ||t/ m || + ||F||v'||y ro ||; 
thus 

4lv m l 2 + "HiHI 2 < M\p\\ 2 + ^-\y m \ 2 \\P\\ 2 + ^Vll 4 + ^ll^ll 2 - (4.4.2) 

From equation (|4.4.2|l and from Gronwall Inequality we can derive the estimates 
(|4.4.3j) and (|4.4.5j) below: 

\y m (s)\ 2 <e X p£^-\\P(t)\\ 2 dt[ 

4(H|P(t)|| 2 + ^-||P(i)|| 4 + \\\F{t)\?) dt) . (4.4.3) 

Since P S L 4 (0, T, D(A)) C L 4 (0, T, V) C L 2 (0, T, V) and F S L 2 (0, T, V), 
(|4.4.3[) shows that 

y m remains in a bounded set of L°°(0, T, H). (4.4.4) 



1 Where P m is the projection of P onto span{Wi | max{ii, 12} < m}. 
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From 

\y m {T)\ 2 -\y m {Q)\ 2 + v f \\y m {t)\\ 2 dt 
Jo 

< [ T Uv\\P(t)\\ 2 + —\\P(t)\\ i + kF(t)\\ 2 )dt+ [ ^-\y m (t)\*\\P(t)\\ 2 dt. 
Jo v v ' Jo v 

(4.4.5) 

By (|4.4.4|l the last integral is finite. Then (|4.4.5|l shows that 

y m remains in a bounded set of L 2 {0, T, V). (4.4.6) 
The rest of the proof is completely analogous. 

Remark 11. As we can see in the estimates (I4.4.3|) and (|4.4.5J) . considering 
P G i 4 (0, T, V) is sufficient to guarantee existence of weak solutions. 

4.5 Uniqueness. 

Theorem 4.5.1. The solution of Problems \4- 3 '.1\4- 3 '.H^ given by Theorem \4-4-l\ 

is unique. Moreover it is a.e. equal to a continuous function from [0, T] into 
H. 

The continuity follows from 

. A(y + P) G i 2 (0, T, V') 

• F G L 2 (0, T, V') 

. B(y + P) G L 2 (0, T, V'). 

Indeed these expressions imply 

y G L 2 (0, T, V) k y' G L 2 (0, T, V). 

The first two expressions we already know to be true. The last one follows from 

\\B(y + P)\\<C\y\\\P\\+C\\P\\ 2 ; 

P G L 4 (0, T, D(A)) c L 4 (0, T, V); y G L°°{0, T, H). 

To conclude the uniqueness we consider two solutions y and z of problems 
(I4.3.1|l - (|4.3.2() . The difference w := y - z satisfies 

w' = ~vAw - B(y + P) + B{z + P) 

from which we derive 

^-\w\ 2 + 2v\\w\\ 2 < 2C|t«|||«;||||a: + P|| 

^M 2 <^M 2 lk + ^ll 2 

dt 2v 

H*)| 2 < \w(0)\ 2 exp f ^-\\z(t) + P(t)\\ 2 dt = 0. 
Jo Lv 
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4.6 Continuity. 

Theorem 4.6.1. The map 

Y : H x L 2 (0, T, V') x L 4 (0, T, D(A))x]0, +oo[ -> C([0, T], i?) 

is continuous. Where y is the unique solution of problems H4.3.1|) - (|4.3.2J) corre- 
sponding to the data (jjq, F, P, v). 

Proof. The proof is analogous to that of Theorem 14.1.21 We fix a quadruple 

(y , F, P, v) e H x L 2 {0, T, V) x i 4 (0, T, D(A))x]0, +oo[ 
and e > 0. Then consider another quadruple 

(z Q , G, Q, 7])eHx L 2 {0, T, V) x L 4 (0, T, D(A))x]0, +oo[. 

Put 

V := ¥(2/0, F, P, !/) & z := Y(z„, G, Q, rj) 

so that 

?/ + vA{y + P) + B(y + P) = F k z' + f]A(z + Q) + B{z + Q) = G. 
Putting w := z — y we obtain 

w' = G-F - rjAw - r)A(Q - P) + (v - n)A(y + P) - B(z + Q) + B{y + P). 
Taking the scalar product with w we obtain 

< w', w > = < G - F, w> -?7||w|| 2 - r]((Q - P, w)) 

+ - rf){{y + P, wj) + b{y + P, y + P, w) -b(z + Q,z + Q, w). 

(4.6.1) 

Now we estimate the last term of the last equality: 

b(y + P,y + P, w) -b(z + Q, z + Q, w) 

= b{y, y, w) + b{y, P, w) + b{P, y, w) + b{P, P, w) 
- b(z, z, w) - b{z, Q, w) - b(Q, z, w) - b(Q, Q, w) 

= b(y, Vi w ) ~ b (z, z i w ) + b (y, p , w ) ~ b (z, Q> w ) 
+ b(P, y, w) - b{Q, y, w) + b(P, P, w) - b{Q, Q, w). 



We put 



z 2 

Zz 
Za 



= Hy, y, w) - b{z, z, w) 
= b(y, p > w ) ~ b(z, Q, w ) 
= b(P, y, w) - b(Q, z, w) 
= b(P, P, w)-b(Q, Q, w). 



Then 



\b(y + P,y + P,w)-b(z + Q, z + Q, w)\ <\Z X \ + \Z 2 \ + \Z 3 \ + \Z 4 \ (4.6.2) 
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Since 



\ z i\ = My, y> w ) - b (z, z , w )\ = I - b(w, y, w )\ 

< CHIMIN, by(H3J (4.6.3) 

\ z 2\ = \b(y, P, w) - b(z, Q, w)\ = \b(y, P, w) - b(y, Q, w) - b(w, Q, w)\ 
= \b(y, P-Q,w)- b(w, Q,w)\< \b(y, w,P-Q)\ + \b(w, Q, w)\ 

< C\y\\\w\\\P-Q\i\P-Q\f 2] +C\w\\\w\\\\Q\\ by tHOt and (ITOl . 
<Ci\y\\\w\\\P-Q\ [2] + C\w\\\w\\\\Q\\, by|-|<||-|| 2 (4.6.4) 

\ z 3\ = \b(P, y, w) - b(Q, z, w)\ = \b(P, y, w) - b(Q, y, w) - b(Q, w, w)\ 

= \-b(P-Q,w,y)\< C\P - Q|*|P - Q\f 2] \\w\\ \y\ by jE3 

<Cl|P-Q| [2 ]HIM: (4-6.5) 

\Z 4 \ = \b(P, P, w) - b(Q, Q, w)\ = \b(P - Q, P, w) + b(Q, P, w) (4.6.6) 
- HQ, Q, to) | < \b(P - Q, P, w)\ + \b(Q, P-Q, w)\ 

<C||P-Q||||P||HI + C||Q||||P-Q||H| bv E1H1) . (4.6.7) 

From EO, and glHT) we obtain 



-||uf < 2\\G - F\\ v >\\w\\ - 2r,\\w\\ 2 + 2 V \\Q - P\\\\w\\ + 2\v - r)\\\y + P\\\\w\ 
ill 

+ 2C\w\\\w\\\\y\\ + 2C 1 |y|||to|||P - Q\ [2] + 2CMIMIIIQII 
+ 2C' 1 |P-Q| [2] |H|| 2 /|+2C||P-Q||||P|||| U ;||+2C||Q||||P-Q|||| U ;||; 



Hence 

^-H 2 + v\\w\\ 2 < -\\G - f\\ 2 v + 8r?||Q - P|| 2 + -W- v\ 2 h + p\\ 2 

at rj r] 



Now if rj satisfies 



we have r\ €]-|, 4r[. Hence, from the last equation we obtain 



^H 2 + ^Hl 2 < -||G - Pf v , +\Q- P|f 2] (12„C 2 + ^4C?|y| a ) 



+ ^l^-^l 2 b+P|| 2 +^c 2 H 2 (bll 2 +IIQII 2 )+^C' 2 ||p-Q|| 2 (||p|| 2 +||Q|| 2 ). 

(4.6.* 
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By the Gronwall Inequality: 

Hs)\ 2 < expQf ^C 2 (\\y(t)f + \\Q(t)\\ 2 ) dt) (\w(0)\ 2 

+ f \^\\G(t) - F(t)\\l, + \Q(t) P(t)\f 2] (12uC 2 + ^Cf\y\ 2 ) 
Jo L v v 



10 

+ ^\u- V \ 2 \\y(t) + P(t)\\ 2 + ^C 2 \\P(t) Q(i)|| 2 (||P(*)|| 2 + \\Q(t)\\ 



Now we consider two cases P = & P^O. 
If P ^ and if Q satisfies 

\\P - Q|U 4 (0,T,£>(A)) < 



dtj. 
(4.6.9) 



IPI 



L 4 (0,T,n(A)) 



and since ||j/||l°°(o,t, h) =. D < +oo we have that 



|w(s)| 2 < exp 



+ C 4 



C 3 (/ o \\y(t)\\ 2 dt+(J o \\P(t)\\ 4 dty 
jT \\G(t) - F(t)f v , dt + (jf \Q(t) - P(t)\f 2] dt) ' 



+ / \\y(t) + P(t)\\ 2 dt+ / ||P -OH 4 * 



IPII 4 * 



Thus 



\w(t)\* < C 5 {\w(0r + \\G- F\\l^ TtVI) + \\Q- P\\i H0 ^ D(A)) 

+ \V - V\ 2 + \\Q - P\\ 2 LH0,T.D(A)^j ■ 

We rewrite the last equation as 



Ht)| 2 < C 2 6 (J W (0)| 2 + \\G- F||i 2(0 , T , vl) + WQ- P\\i HQ ,T.D(A)) + W- V[ 
Then (in the case P ^ 0), if the quadruple (zo, G, Q, 77) satisfies 
\v-T)\< min{-, — }; 

ME, ^ ■ f ll P IU 4 (0,T,_D(A)) e \ 

\\P-Q\\l*(o,t,d(A)) < mm| , — J; 

\zo - 2/0 1 = |w(0)| < ^r; ||G- P||i,2(o,t,v') < ^r; 

we have that |u>(s)| < e. Therefore we have the continuity of Y in all quadruples 
such that P^O. Note that Cq depends only in the fixed quadruple (y , F, P, v) 
where we are studying the continuity. 2 

2 C(; depends in R and T too, but T and R are fixed in the statement of the Theorem. 
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In the case P = we can start from a bound for Q, say that Q satisfies 
II ( 5IIl 4 (o,t, d(A)) < CV (so that the argument of the exponential is bounded by 
a constant depending only on C<j and (yo, F, P, v), but not on Q). Then from 
(|4.6.9[1 we can conclude that 

\w(t)\ 2 <C a (\w(0)\ 2 + \\G-F\\ 2 L2{0 ^ vl) 

+ IIQII L 4 (0, T, D(A)) "T" L 4 (0,T, D(A)) I ■ 

If IIQI|l 4 (o.t,.d(y1)) < 1 we can rewrite the last equation obtaining 

\w(t)\ 2 < C 9 2 (j W (0)| 2 + ||C — F|| 2 L2(0 , T ,v. + \\Q\\h {0 ,T,D(A)) + W~V\ 2 ^) 
So in the case P = we choose, for the quadruple (zo, G, Q, rf), the bounds 

V £ £ 

\v-7]\ <min{-, — }; ||Q||l 4 (o : t,d(A)) < min{l, — }; 
\z -y \ = \w(Q)\ < \\G - F\\ L 2 {0)T)VI) < 

and Y is continuous in this case too. □ 
Theorem 4.6.2. The map 

Y 2 : H x L 2 (0, T, V) x L 4 (0, T, D(A))x]0, +oo[ -> L 2 ([0, T], V) 

(yo, -F, P,v)^y 

is continuous. Where y is the unique solution of problems (14.3. lfl - (|4.3.2I> corre- 
sponding to the data (yo, F, P, v). 

Proof. We Fixe > 0, (y , F, P, v) G HxL 2 {0, T, l/') xi4 (0, P D(A))x}0, +oo[ 
and put y := Y 2 (y , F, P, v). 

Let (zq, G, Q, rf) be another quadruple in the product H x L 2 (0, T, V') x 
i 4 (0, T, L>(A))x]0, +00 [. Put z := Y 2 (z , G, Q, rf) and w:= z - y and like in 
the proof of Theorem 14.6.11 we arrive to l|4. 6.8(1 : 



|H 2 + ^IIHI 2 < ^\\G- Ff v , + ^\\Q- P\\l(l + AC!\y\ 2 ) 

+ ^|^-r / | 2 || y +p|| 2 +^c 2 |H 2 (lly|| 2 +lig|| 2 )+-C' 2 ||P-Q|| 2 (||P|| 2 +||Q|| 2 ). 

V V V 

Integrating over [0, T] 



,(T)\ 2 -\w(0)\ 2 + -\\w(t)rdt< 



' '-'|2 



2 j 



^||G(t)-F(i)|| 2 y , 



16 16 

+ -||Q(t) - P(t)|| 2 (l + 4C 2 |y(t)| 2 ) + -\ v v \ 2 \\y(t) + P(t)\\ 2 

^C7 2 |^( £ )| 2 (||^(t)|| 2 + ||Q(t)|| 2 )+^C7 2 || J-( £ )-Q(t)|| 2 (||^(^)|| 2 + ||C5(t)|| 2 ) 



dt. 
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Looking at the last two terms we see that 

\\y(t)\\ 2 + \\Q(t)\\ 2 dt = [ T \\y(t)\\ 2 dt + J T \\Q(t)\\ 2 dt 

2 

li 4 (0,T, D(A)) 



<D 1 + D 2 \\Q\\ 2 



and 



T r r 

2iid/ + mi2j + , / \\r>/4.-\ nc + Mi2|i / n/ + Mi2 . 



\\P(t)-Q(t)\\'\\P(tWdt+ / \\P(t) - Q(t)\nQ(tW dt 

Jo 

<D 3 \\P - Q\\l*(0,T,D(A)) (\\P\\l*{0,T,D(A)) + IIQ|ll 4 (0,T,£>(A)) 

so, we arrive to 

\\w(t)fdt <d\\G- F\\h {0 ,T, V+C4Q- P\\l H 0,T,D(A) 

+ C 3 \v- ri\ 2 +C 4 \w\ 2 c([0 T] H) {1 + \\Q\\h {0 ,T:D(A)) 

+ C 5 \\P — Q\\ L 4 (o,T, i3(A)(ll-f II L 4 (0,T, D{A) + WQW 

where the constants do not depend on the quadruple (zq, G, Q, rf). 
If we choose Q close to P, say \\P — Q\\l 4 (o,t. d(A) < 1 we have 

T |kWII 2 ^<Ci||G-#||i 2(0 ^ !y ,+C 2 ||Q-P||| 4(0!T ^ (A) 
+ C 3 \v- n\ 2 + D 4 \w\ 2 c([Q T] H) (l + ||P||L 4 (o,T,r>(A)^) 

+ D S \\P- Q\\l*(0,T,D(A) ( 1 + \\ p \\lho,t,d(A) 



or 

rT 







IHtWdt^CiWG-FWt^y,) 



+ D 6 \\Q-P\\ 

L 4 (0, T, D(A) 

+ C 3 \v-n\ +D 7 \w\ cmTl 



H) 



(where the constants do not depend on the quadruple (zq, G, Q, rf)). 

By Theorem l4.6.1l there is a 8 > such that if both \yo — Zq\, \\G— -F||z,2(o,t, v)> 

\\P-Q\\lHo,t,d(A) and \q-v\ are less than <5, then Hc([o, t], h) < ^Jrb;- Hence 

for some <5i smaller that 8 we have that if both |yo — zq\, \\G — F\\l 2 {o,t,v')i 

i 

\\P ~ QI|l 4 (o,t, d(A) and |?7 — v\ are less than 8\, then ( J Q T \\w(t)\\ 2 dtj < e. 
Therefore the map Y 2 is continuous. □ 

Remark 12. As said in R.emark 1 1 1\ considering P 6 £ 2 (0, T, V) is enough 
to guarantee the existence of weak solutions but, to have the continuity of Y 
or Y 2 we have to consider P varying in L 2 (0, T, D(A))-norm. Variation in 
L 2 (0, T, V)-norm seems to be not stronger enough. 
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4.7 Strong Case. 

Consider the equivalent problems 14 . 7 . 1 1 and 14 . 7 . 21 
Problem 4.7.1. Given 

F e L 2 (0, T, H), P G £ 4 (0, T, D(A)) (4.7.1) 
& 

yo € V (4.7.2) 

to /ind 

?y G L 2 (0, T, D(A)) n L°°(0, T, V) (4.7.3) 
satisfying (in the distribution sense) 

yvev : 

j t {y, v) + v({y + P, «)) + % + P, y + P, v) = (F, v), (4.7.4) 
and 

y(0) = yo- (4.7.5) 

Problem 4.7.2. Given 

F G L 2 (0, T, if), , P G L 4 (0, T, £>(A)) (4.7.6) 
& 

?yo G V, (4.7.7) 
to find 

y G L 2 (0, T, £>(A)) n i°°(0, T, V), y' G L 2 (0, T, if) (4.7.8) 
satisfying 

y' + vA(y + P) + B(y + P) = F on ]0,T[, (4.7.9) 
and 

1/(0) = wo. (4.7.10) 
The equivalence of these problems follows from 

• A solution of Problcm l4.7.2l is a solution of Problem l4.7.1l 

• Given a solution of Problcm l4.7.1l 

- - F G L 2 (0, T, if); 

A(y + P)eL 2 (0, T, ff). 

Indeed ||A(y + P)|| 2 = \y + P\f 2] and y + P € L 2 (0, T, D(A)); 

-% + P)eL 2 (0,r,4 

See Remark H~3l below. 

So that F - vA{y + P) - B(y + P) G L 2 (0, T, if) C L 2 (0, T, V). ByEHand 
y G L 2 (0, T, V) we have y' = F-vA(y+P)~B(y+P) a.e., andy G C([0, T], V) 
"a.e.". 

Remark 13. Py definition and trilinearity we have 

B(y + P)(h) = h{y + P,y + P, h) 

= b(y, y, h) + b(y, P, h) + b(P, y, h) + b(P, P, h) 
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or 



B(y + P)= B(y, y) + B(y, P) + B(P, y) + B(P, P). 
From <|2.9.2|) we obtain 



(4.7.11) 



\B(y, v)\<C\y\*\\v\\\v\f 2]] 

\B(y, P)\<C\y\?\\y\\i\\P\\i\P\f 2] ; 

\B(y, P)|<C|P|i||P||i|M|5|j/|| ]; and 

\B(P,P)\<C\P\i\\P\\\P±. 



Thus for some constant C\ 



y)l 4 < Cilly|| 6 |y|f 2] e L^O, T, R); 
\B(y, P)| 4 < ai||y|| 4 |P|f 2] 6^(0, T.R); 
\B(y, P)\ 2 < Ci|P|f 2] ||y|||y| [2 ] € L\0, T, K); and 
\B(P, P)| 2 <Ci|P|f 2l G i x (0, T, M). 



5o, P(y, y) andB(y, P) are elements of L (0, T, P) and, P(P y) and B(P, P) 
are elements o/L 2 (0, T, P) f/uts, 6y (|4. 7. 1 l|l . B(y + P) is a sum of elements of 
L 2 (0, T, P). 

4.7.1 Existence. 

Theorem 4.7.1. Given F, P and uo satisfying 14.7.6J1 and (|4.7.7|l . There is 
at least one function y satisfying (|4.7.8|I - I|4.7.1U|I . 

Proof. As we have done for Thcorem l4. 4. II starting from approximate solutions 
we arrive to the conclusions 14.4.4|) and 14.4.6|) . 

An approximate solution y m := Y. m ^{i l . A2 }< m vT( t ) w i satisfy 



((y m Y, v) + v{A{y m + P m ), v) + (B(y m + P m ), v) = (F, v), v e V, 



((y™)', Ay m ) + v{A{y m + P m ), Ay m ) + (B(y m + P m ), Ay m ) = (F, Ay m ). 



(B(y m + P m ), Ay m ) 
=b(y m + P m , y m + P m , Ay m ) 

=b(y m , P m , Ay m ) + b(y m , y m , Ay m ) + b(P m , P m , Ay m ) + b(P m , y m , Ay m ). 



setting y™iWi for v and summing up we arrive to 



(4.7.12) 



We have that 



(A(y m + P m ), Ay m ) 




and for the term with B we have 
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By (gUSH we have 

\b(y m , P m , Ay m )\ < C\y m \*\y m \f 2] \\P\\\y m \ [2] = C\y m \^\\P\\\y 

\b(P"\P m ,Ay m )\<C\P\^\P\f 2] \\P\\\y m \ [2] 
\b(P m , y m , Ay m )\ < C\P\i\P\f 2] \\y m \\\y m \ [2 

From l|4.7.12jl we obtain 

d 
dt 



m\ 2 
[2] 



3 

m| 2 

21 



-||2/ m || 2 + 2u\y m \f 2] < 2v\P\ [2] \y m \ [2] + 2\F\\y m \ [2] ; 



+ 2q y m | 5 ll^llly m lf 2] + 2C| 2 ;"^|| 2/ m ||| y ™| [ | ] 

+ 2qP|5|P|| ] ||F|||y m | [2] +2C|P|^|P|| ] || 2/ m |||y™| [2] ; 

by Young inequalities 
d 



Jt i\y m \\ 2 + Hy m l[ 2] <c 1 |p|f 2] + c 1 |p| 2 ; 



2| |y m || 4 



+ C 2 \y m \ 2 \\P\\ 4 + C 2 \y r > 
+ C3|P||P| [2] ||P|| 2 + C3|P||P| [2] ||2/ m || 2 . (4.7.13) 

By Gronwall Inequality 



||y m WH 2 < exp(jf C 2 \y m (s)\ 2 \\y m (s)\\ 2 + C 4 \P(s)\f 2] ds) (j|y(0)|| 2 

+ f Ci|P( S )|f 2] + CxlP^I 2 + C 2 | y '»( S )| 2 ||P( S )|| 4 + C 5 \P(s)\f 2] ds 



By 1033)1, [jlXHJl . P S L 4 (0, T, P(A)) and P G L 2 (0, T, P) we have 
||t/ m (i)|| 2 < Ki for some constant K\. 

Hence 

y m remains in a bounded set of L°°(0, T, V). (4.7.14) 
Integrating 14.7.1311 over [0, T] we obtain 



|| ym(T) ||2_ ||ym(0)|| 2 + / ly m {t)l 2 2]dt < / Cl \P{t)\\ 2] + C 1 \F(t)fdt 

Jo Jo 

+ f C 2 \y m \ 2 \\P\\ 4 + C 2 \y m \ 2 \\y m \\ 4 + C 5 \P\f 2] + C 6 |P| 2 2] ||y m || 2 dt. 

Again, by (l4"4~4l) . BXTHl . (I4.7.14L P S L 4 (0, T, and P e L 2 (0, T, H) 

we arrive to 



T 

"' ' ' ' J for some constant K? 



\y m (t)\{ 2] < K 2 

Hence 

- 2 



y m remains in a bounded set of L z (0, T, D(A)). (4.7.15) 
The rest of the proof is analogous to that of Corollary 14.1.11 □ 
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Remark 14. Contrary to the case of weak solutios (see Remark \ll\) . we can not 

guarantee from the previous sketch that considering P £ L 2 (0, T, V) is enough 
to guarantee the existence of a strong solution. 

4.7.2 Uniqueness. 

Theorem 4.7.2. The solution of Problems [4-7.1^4-7-^ gitien by Theorem \4-7.1\ 

is unique. Moreover it is a.e. equal to a continuous function from [0, T] into 
V. 

Proof. The uniqueness follows from Theorem 14.5.11 and from the fact that a 
solution of problems 14.7. 1I4.7~51 is a solution of problems 14.3. 114. 
By y' £ L 2 (0, T, H) we have y' £ L 2 (0, T, D{A)'). Since y £ L 2 (0, T, D{A)) 
and the inclusions 

D(A) eye D(A)' 

are dense and continuous we have that y £ C([0, T], V). [See Lemma 1.2 in [T5] 
section III. 1] . □ 

4.7.3 Continuity. 
Theorem 4.7.3. The map 

Y s : V x L 2 (0, T, H) x L 4 (0, T, D(A))x}0, +oo[ -> C([0, T], V) 

(2/0, F, P, f)^y 

is continuous. Where y is the unique solution of problems l|4.7.1|l - l|4.7.2|l corre- 
sponding to the data (jjq, F, P, v). 

Proof. The proof is analogous to that of Theorem 14.6.11 We fix a quadruple 

(y , F, P, v) £ V x L 2 (0, T, H) x L 4 (0, T, D(A))x}0, +co[ 
and e > 0. Then consider another quadruple 

{z , G,Q,n)£Vx L 2 (0, T, if) x i 4 (0, T, D(A))x]0, +oo[. 

Put 

y := Y s (y , F, P, 1/) & z := Y s (z , G, Q, rj) 
Putting w := z — y we obtain 

w' = G-F- rjAw - r]A{Q -P) + {v- rj)A{y + P) - B{z + Q) + B(y + P). 
Taking the scalar product with Aw we obtain 

(to', Aw) =<G- F, Aw > -f/Mfa] ~ ViQ - p - w h] 

+ (v - t])(y + P, io)p] +b(y + P,y + P, w) -b{z + Q, z + Q, w). 

(4.7.16) 

Now we estimate the last difference: 

b(y + P,y + P,w)-b(z + Q,z + Q, w) 

= b(y, y, Aw) + b(y, P, Aw) + b(P, y, Aw) + b(P, P, Aw) 
- b(z, z, Aw) - b(z, Q, Aw) - b{Q, z, Aw) - b(Q, Q, Aw). 
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We put 

Z\ := b(y, y, Aw) — b(z, z, Aw) 

= Hv, y, Aw) - (b(w, z, Aw) + b(y, z, Awn 

— b(y, y, Aw) — (b(w, w, Aw) + b(w, y, Aw) + b(y, w, Aw) + b(y, y, Aw)) 

= —b(w, w, Aw) — b(w, y, Aw) — b(y, w, Aw); 
Z 2 := b(y, P, Aw) - b(z, Q, Aw) 

= b{y, P, Aw) - (b(w, Q, Aw) + b(y, Q, Aw)) 

= b(y, P, Aw) - (b(w, Q, Aw) + b(y, Q - P, Aw) + b(y, P, Aw)) 

= -b(w, Q, Aw) - b(y, Q-P, Aw); 
Z 3 := b(P, y, Aw) - b(Q, z, Aw) 

= b(P, y, Aw) - (b{Q, w, Aw) + b(Q, y, Aw)) 

= 6(P, y, Aw) - (b(Q, w, Aw) + b(Q - P, y, Aw) + b(P, y, Aw)) 

= -b(Q, w, Aw) - b(Q - P, y, Aw); 
Z± := b(P, P, Aw) - b(Q, Q, Aw) 

= b{P - Q, P, Aw) + (b(Q, P, Aw) ~ b{Q, Q, Aw)) 

= b{P - Q, P, Aw) + b(Q, P-Q, Aw); 

We have 

\b(y + P, y + P, w)-b(z + Q, z + Q, w)\ < \Z X \ + \Z 2 \ + \Z 3 \ + \Z±\ (4.7.17) 
and, by l(!HH3|) 

l^al < Cl^lil^lljll^Hl^ip] +^1^15^1^11^111^1^ +C| 2 /|ib|| ] || W |||^| [2] , 

(4.7.18) 

I^I^CHiHljHQiiH^ + aiyi^MljiiP-oHHp], (4.7.19) 
13,1 <c|Q|*|Q|| ] HIH[2] + c|i'-Ql i |J > -QlJ ] NIH[2] ) (4-7.20) 

|Z 4 |<C'|P-Q|i|P-Q|| ] ||P||Hp ] +C'|Q|2|Q|| ] ||P-Q||H [2] . (4.7.21) 
From l|4.7.17|) and H4.7.16|l we obtain 

^|M| 2 < 2|G - P|M[ 2] - 2»?Mf 2] + 2 V \Q - P| [2] M[2] 

4 

+ 2| I ,-7 ? ||y + P| [2] | W | [2] +2^|Z i |. 
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If 77 satisfies \rj — v\ < §, i.e., | < 77 < we obtain 

^||HI 2 + HH?a] < 2|G - F|H[ 2] + MQ - P\[2]\w\ [2] 

4 

+ 2|i/-» 7 |||, + P|p ] H[a] +2^|^| 

i=l 

so, applying Young inequalities, for appropriate constants (independent of the 
quadruple (z , G, Q, 77)) we have 

|lMI 2 + ~M? 2] < Ci|G - P| 2 + C 2 |Q - P|f 2] + C z \v - nf\y + P|f 2] 

. 1 2 1 1 „ ,. 1 1 4 1 |„,,|2||„.||4 , L.IL.I ||„..||2 , U,.|2|| / -i||4 i L.IL.I II D /OI|2 







+c 5 r 

Jo 



+ c, (MiHI 4 + \wW + \yWy\mWwW 2 + MlQII 4 + \v\\v\m\\P Oil 

Hg||g| M i| w || 2 +|p-oi|p-oi [2] ||y|| 2 +|p-oi|p-oi [2] ||p|| 2 +ioiiQ| [2] ||p-oii 2 

(4.7.22) 

By the Gronwall Inequality: 

M*)ll 2 

<ex P f f T c 4 (Hf)| 2 |Kt)|| 2 + \y(t)\\y(t)\ m + dtj (\\w(0)\\ 2 

IOW - P(t)\f 2] (1 + \y(t)\\y(t)\m + h(t)f (4.7.23) 

+||P(t)|| 2 + |Q(i)||Q(%]) + \G(t) - P(i)| 2 + \ v - v \ 2 \y(t) + P(t)| a p] 
+K*)| 2 ||y(i)|| 4 + Kt)| 2 ||Q(t)|| 4 dtj. (4.7.24) 

Given a constant E > 0, by Theorems 14.6.11 14.6.21 and by the existence of a 
constant C such that 3 

|a|<C|M| & \\f\\ L * { o,T,v>)<C\\f\\ L > {0 , T , H) 

for alia € V", / € i 2 (0, T, P). 

there is <5 > such that if both 1 1 3/0 — ^0 1 1 , \\G- F\\ L 2 {0 ^ TtH) , \\P-Q\\l*{o,t,d{a)) 
and \rj — v\ are less than S we have 

\H\c([o.T].H) <E and f \\w(t)\\ 2 dt < E. 

Jo 

Therefore, if we choose in addition S < £ , we have 



3 The existence of such a constant comes from the continuity of the inclusions V C H C V' . 
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\\w( s )\\ < D i 
Jo 



and, because 

l-T 



\\wmV + \\G-F\\h{o,T, H) 

\Q(t) - P(t)\[2] (i + lv(*)l[2] + |P(*)|f a] + \Q(t)\f 2] 

+ \u- r/\ 2 + \\w\\ 2 c{[Q T] H) {l + ||Q||t*(o,r,r»(A))) 



dt 



J \Q(t) - P(t)\f 2] (l + |y(t)| [2] + |P(*)|f a] + \Q(t)\f 2] ) dt 

<\\Q - P\\%*(0,T, D(A)) M|1|U 2 (0,T,K) + l|y|U 2 (0,T,D(A)) 
+ 11^11 L 4 (0,T, D(A)) + WQW L 4 (0, T, D(A)) I 

<Kx\\Q - P\\ 2 LH ^ TtD(A)) + (5+ \\P\\l*{o,T:D(A))) 2 ^ 
<K 2 \\Q - P\\l*(o,t,d(,A))- 

we have 



\w(s)\\ 



< K s 



(4.7.25) 



(with K3 independent of the quadruple (zo, G, Q, rj)). 

Then for some 5\ smaller than 8, and using Theorem l4.6.1l we have ||w(s)|| < e 
if both ||y -z |j, \\F - G\\ L 2 {0 T H) , \\P - Q\\ l *(o,t, d(A)) and |r? - v\ are less 
than 5\. Thus the map Y s is continuous. □ 

Theorem 4.7.4. The map 

Y 2s : V x L 2 (0, T, H) x i 4 (0, T, D(A))x]0, +oo[ L 2 ([0, T], D(A)) 

(yo, p, P,v)^y 

is continuous. Where y is the unique solution of problems (I4.7.1[l - H4.7.2|) corre- 
sponding to the data (yo, F, P, v). 

Proof We fix (y , F, P, v) e VxL 2 (0, T, H)xL A (Q, T, D(A))x]0, +00 [, £ > 
and, put y := Y 2s (y , P, -P, f)- 

Let (zo, G, Q, ri) be another quadruple in the product space V x L 2 (0, T, H) x 
i 4 (0, T, L>(A))x]0, +00 [. Put z := Y 2s (z , G, Q, 77) and to := z-y and like in 
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the proof of Theorem 14. 7.31 we arrive to H4.7.22[l : 



^IM! 2 + 2 M[2] <Ci\G- F\ 2 + C 2 \Q- P|f 2] + C 3 \v - V \ 2 \y + P|f 2] 

+ c* 4 ^M 2 IHI 4 + \w\ 2 \\y\\ 4 + \y\\y\[2]\\w\\ 2 + M 2 IIQII 4 
+ \y\\y\[2]\\P - Oil 2 + IQIIQI[ 2 ]IHI 2 + |JP-Q||JP- QlmWvW 2 

+ |P-Q||P-Q| [2] ||P|| 2 + |Q||Q| [2] ||P-Q|| 2 



^IM| 2 + 2 kl[ 2] <Ci\G- P| 2 + C 2 |Q - P| 2 2] + C 3 \p - V \ 2 \y + P| 2 2] 
+ C b (^wf + \\w\\ 2 + \y\f 2] \\w\\ 2 + WwfWQf + \Q\f 2] \\w\\ 2 

+ l^-QI[ 2 ](i + |y|[2] + |P| 2 2] + IQI 2 2] 



(where the constants are independent of the quadruple (zq, G, Q, r/)). 
Integrating over [0, T] 



J (t)\[2] dt < D \\w\\c([q,t\,v) + \\F ~ ^\\l 2 (o,t, h) + \\Q ~ -P|Il 4 (o,t, D{A)) 

+ W-V\ 2 + \\w\\ 2 C([0,T],V)(\\ w \\c([0,T],V) + 1 + IMll^O.T, 23(A)) 

+ \\Q\\1 4 (0,T,D(A)) + II ( 3IIl 4 (0,T,D(A))) 

+ ll-P - Q\\h(0,T,D(A)) (i + \\y\\L2(0,T,D(A) 

+ \\P\\ L 4 (0,T, D{A)) + IIQII L 4 (0, T, D(A)) 



By Theorem 14.7.31 given a constant E > 0, there is a (5 > such that, if both 
|] 2/0 - zo||, ll-P - QIIl 4 (o,t,d(A)), IIP - G\\ L 2 {0 T H) , and |r/ - i/| are less than S 
we have 



IMIc([o,t], V) < E. 
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Thus, 



Ht)|f 2] dt < D 1 \\F- G||i a(0jTiH) + HQ- P\\l, {0 , T ,D(A)) + W- V\ 2 

+ \\ w \\c([0,T],V){\\ w \\c([0,T],V) + 1 + l|2/|lz,2(0,T,Z)(A)) 
+ \ S+ \\P\\L 4 (0 : T,D(A)ij + (S + ||-P||l 4 (0, T, D(A)) 

+ \\ F - ( 3IIl'1(o,t,_d(A))( 1 + \\v\\l 2 (p,t,d{A) 

+ \\ P \\l*(Q,T,D(A)) + ( S + \\ P \\l\0,T,D(A))) ) 



or 



(t)\f 2] dt<D 2 



\\ F ~ G\\l*(0,T,H) + \\Q ~ ^'lli 4 (0,T,_D(A)) 



W - V\ 2 + IHIc([0,T],V) 



(with Z?2 independent of the quadruple (zq, G, Q, rf)). 

Hence for some Si smaller than S we have \w\%, < e 2 if both ||yo — £o||j 

||F - G|| L 2( 0i t, jj), II -f Qll l 4 (o, t, d(A)) an d I 7 ? - H are less than Thus the 

map Y 2s is continuous. □ 



4.8 Continuity in Relaxation Metric. 

We begin with a definition: 

Definition 4.8.1. The relaxation metric in ^ 1 ([0, T], R d ) is defined by the 
norm 



\q Wrr ■— max 



to 



9(r) dr 



(4.8.1) 



where \\ ■ || R <j is i/ie Zi norm in 



If nothing in contrary is stated we consider the spaces M d (d € No) endowed 
with Zi-norm — ||x|| R a = ||a;||i 1 := X^f=i Nil- 



Remark 15. It is easy to check that (|4.8.1|l is a semi-norm and, since functions 
in L 1 ([0, T], WL d ) coinciding on a set of measure T are identified we can conclude 
that (|4.8.1|l is a norm. 

Consider, also, the w- relaxation metric on L 1 ([0, T], M. d ) defined by the 
norm 



\g\\wrx ■= max 
te[o,T] 



g(r) dr 



(4.8.2) 
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Remark 16. The identity map 

(V([0, T], R d ), || • |U) ^ (V([0, T], R d ), || • \\ wrx ) 
is continuous. 

Now we fix a finite subset F C Ng, put T := #F and define the space 
Sw '■= span{Wk \ k € F}. From the fact that the dimension of Sw is finite, 
we conclude the equivalence of all norms defined in it. In particular the norms 
|| • IIrj- and | • | [2] are equivalent on From this equivalence we have the 
following: 

Proposition 4.8.1. The map 

l:££x([0, T], M*) - C([0, T],S ¥ ) 

is continuous. Where lufe(i) := f Q Ufc(r) dr and t/ie subscript "wrx" means that 
we are considering relaxation metric on the set L°°([0,T], R^). 

Proof. The continuity of the map 

I:L£x([0, T],R^)^C([0, T],R^) 
(■Ufe(t))feeF ^ (Iffc(i))fceF 

is trivial. □ 

Recall that by definition, the map § of Corollary 14.1.21 gives us the weak 
solution, belonging to C([0, T], H), of the NSE for an initial data in II := 
H x L 2 (0, T, V) x L°°([0, T],S r ) x R+. Changing the topology on the third 
factor of the previous product to the w-relaxation one we arrive to the space 
L^ rx ([0, T], R-^) and we define the function S> wrx as the function defined in the 
product n wrx := Hx L 2 (0, T, V') x L™ rx ([0, T],M?) x R+ and taking the same 
values as §. 

Proposition 4.8.2. The map S wrx is continuous. 

Proof. As in the beginning of section 14.21 we put Y wrx := E> wrx — I - But 
here we want an equality between functions defined on Tl wrx so, we just put 
Io(wq, F, v, v) :— Iv. 

By Proposition 14.8.11 and Theorem 14.6. II the map 

Y wrx '■ n wrx — ► C([0, T], H) 

{u , F, v, v)^Y{u , F, Iv, v) = YoF(u 0) F, v, v) 
is continuous. Where I°(uo, F, v, v) := (uq, F, Iv, v). 

By the equality S wrx = Y wrx + I we conclude the continuity of S wrx . □ 

Analogously, using Proposition 14 . 8 . ll and theorems 14.6.21 E.7.31 and 14 . 7 . 41 we 

can prove the continuity on relaxation metric of the maps §2 ; §s and §2 S arriving 
to the Proposition 
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Proposition 4.8.3. The maps S wrx , S>2wrx, ^swrx ^2swrx are all continuous. 
By Remark [TBI we obtain 

Corollary 4.8.4. The TYl&pS §?-:r ; *^2rx ? ^srx ^2srx 

are all continuous. 

Remark 17. //, instead of a subspace spanned by a finite number of eigen- 
f unctions, we consider any (with either finite or infinite dimension) subspace 
F C D(A) and, instead of defining the relaxation metric using the l\-norm we 
define it on L 1 (0, T, F) using the norm induced by D{A): 

and, its corresponding weak form by 

h\\wrx[2\ ■= max 

t£[0,T] 

we still to have the continuity of the corresponding maps S rx pi, S 2rK [2], § S rx[2] 
S 2sra: [2], i-e-, the solutions of the NSE vary continuously when the control varies 
continuously in rx[2]-metric. 

The reason why we have considered a finite space is that we are interested in the 
case where the control in a space spanned by a finite number of eigenf unctions 
and, the reason to consider l\-metric in the definition of relaxation metric is 
that it will be convenient later. 

4.9 A Remark on Finite System versus Infinite 
System. 

From now we consider only strong solutions and, the external force F and the 
coefficient of viscosity v are fixed. 

Definition 4.9.1. Given T > 0, we say that the N-S system is time-T ap- 
proximately controllable in observed component if for any <fi S V and 

any finite subset of modes O C Ng ; the projection of the closure of the attainable 
set at time T from (f> S V onto span{Wk \ k £ O} is surjective. It is time-T 
controllable in observed component if the projection of the attainable set 
at time T from (f> £ V on span{Wk \ k € O} is surjective. 

In this section we ask ourselves if can we prove approximate controllability 
in observed component similarly as we have proved approximated controllability 
of Galerkin approximations. 

As we have seen in (|3.4.2|l from the N-S equation 

u t + vAu + P v Bu = F + v 5 



j(t) dr 



[21 



g(r) dr 



[21 



4 These "rx"-maps are defined similarly as the "wrx" ones, just considering the "rx"- 
topology in the factor of essentially bounded functions. 
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we can derive the infinite ODE system 

Eu m u n C m n 
= (n - m) 
o k 

m,n£N„ 
m<n 

{n(++)m) + =k 

+ f — (n - m)sign(ni - rai)sign(n 2 - m 2 ) 

m<n 

(n( jm)"*"— 

\— v u mU n C mn _ 

+ > = (n — TOjsignni — m\ 







m<n 
(n( — h)m) + = fe 



\ ~* u m u n C m n _ _ 
+ 2^ 1 (n — m)sign(n 2 — m 2 ) 

m<n 

(n(H — )m) + =fc 

+ vku k + F k +v k . (4.9.1) 

Now consider the subset JC 1 — {(n,i, n 2 ) £ Nq | ni, n 2 < 3}\ {(3, 3)} as the set 

of forced modes. Put k\ := #/C 1 . 

If we put y := u — Iv we can write u in the form 

u = y + lv (4.9.2) 

and y satisfies 

yt = ~vA(y + Iv) - P v B(y + Iv) + F. 
Writing system (|4.9.1|l (with K. 1 as set of controlled modes) in the form 

{u k = Bk{u) +vku k +F k + v k k e K} (493) 
iik = E>k(u) + vkuk + Fk k ^ K, 1 

in the new variable "y" the system becomes 

Vk=B k (y + Iv)+vk(y k +Iv k )+F k w fe = if fc ^ /C 1 . (4.9.4) 

For any time T > the closure of the attainable set at time T of system 
(|4.9.3Jl — A Uo (it) (T) — and the closure of the attainable set at time T of system 
EM — Au (y)(T) — 6 are related by 



A U0 (u)(T) = A U0 (y)(T) x R Kl = A Uo (y)(T) x R Kl . 



Indeed the inclusion A Uo (u) (T) C „4„ (y) ( r ) x RKl follows from I4.9.2|) and the 
reverse one follows from the density of the map 

L°°([0, T], M Kl ) 9DH (Iv, Iv(T)) e L 4 (0, T, R Kl ) x R Kl 



5 Recall that P v is the projection on the space of divergence free functions. 
6 Note that, by 14.9.21 . at initial time we have ug = yo- 
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and from the continuity of the map 

L 4 (0, T, M K1 ) 3 P h- Y s {u , F, P, u)(T) £ V 

for fixed uq that is a consequence of Theorem 14.7.31 Similarly if we rewrite the 
NS Equation as 

u t + vAu + P v Bu = F + v 1 + v 2 

we obtain 

where A Uo {y l ){T) stays for the attainable set at time T from uq of the system 

y t = -vA(y + lv 2 ) - P v B(y + lv 2 ) +F + V 1 . 
So Factorization works like in the finite dimensional case. 

In the step of convexification a problem arises: In section l3~5l we obtained the 
new vector fields j m ,n from the vector fields n and n and from the ex- 
pressions 

„ W + f~< » W ,/ s x 
! 5 1 = / W + A 7m : n 

and 

; 2 1 = J\ U ) - A7m,n, 

i.e., we have extracted a vector field from the convexification of the vector 
fields f± v x and f± w \ and, in the finite dimensional case we know that we 
can convex without changing the closure of attainable set at time t. From 
Factorization we know that we can follow the vector fields f± v \ and f± w x 
without changing closure of attainable set but, in the infinite dimensional case 
we do not know if we can convex these new vector fields without changing the 
closure of attainable set. 

If using the vector fields f(u) + X^ m .n and f(u) — X-y m ,n we would not change 
closure of attainable set then we could add these new directions to the old ones 
— span(K}) — without changing the closure of attainable set. 

It is known that any control vq £ span{S m n , | (m, n) £ Si, k E K, 1 } 7 can 
be approximated in relaxation metric by controls taking values on {a5 mjn , aeu \ 
(m, n) £ Si, k £ IC 1 } where a £ R is a positive constant depending on vq (as we 
will see in Lemma r4.10.5l below'l. By the continuity of §> srx we do not change the 
closure of attainable set using controls in span{5 m , n , e& | (m, n) £ Si, k £ JC 1 }. 
We have just applied a step of Convexification. Note that this last procedure of 
Convexification is quite different from that after Factorization Procedure that is 
more complicated: For example setting m := (1, 1), n := (2, 1) and considering 
the vector field f v \ ^ (u) = f(u) + Xy mt n + Vi(u, v^ n ) (see (|3.4.4H V we see that 
the candidate to new control — Xj mtn + Vi(u, „), contrary to what happens 
in the last step of Convexification depends on u and does not take values on a 
compact subset, indeed if we put u r := rW m r £ R we obtain 

A7 m ,„ + Vi(u r , u^ n ) 
=\l m .n - Xr[P v B(W n , W m ) + P v B(W m , W n )] + rvm. 



7 See 13.5.41 for the definition of S\. 
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Since A 7m ,„ - Xr[P v B(W m , W n ) + P v B(W m , W n )] takes values on 

span{W {1 ^ 2 ), W(3,2)} 

we obtain rv{X, 1) for the projection of \^ m ,n + Vi(u, „) onto span{Wn t i)}. 
So, the projection goes to oo as r does. 

Without compactness we are not able to use the Approximation Lemma f4.1 0.51 
below. 

In the case of Galerkin approximations from approximate controllability we 
could derive (exact) controllability using the bracket generating property of the 
system but, in the infinite case we do not have a property like that so, we would 
not be able to conclude controllability in observed component immediately from 
approximate controllability (in the case we could somehow prove approximate 
controllability) . 

In the next Section, using some more tools, we prove the so called solid control- 
lability in observed component for system (14.9. 311 which implies controllability in 
observed component. 



4.10 Solid Controllability in Observed Compo- 
nent 

Definition 4.10.1. Let <jp : M — > M 2 be a continuous map between two finite 
dimensional C° -manifolds, C M 1 be an open subset with compact closure 
and, S C M 2 be any subset. We say that <j> (0) covers S solidly, if for some 
C° -neighborhood Af of <jr |q there holds: S C ^>(fi). 

Let O C Nq be the finite set of modes we want to observe and, IIe> be the 
projection map from V onto span{Wk | k € O}. Define, for each T > and 
each finite subset F c N§, the "end point" map 

E T : V x L°°([0, T], K #F ) -> O 

(uo, v) h-> n o S s (u , F, v, v)(T). 

For any N G No define, also, the system 



I u. fc - B k (u) + vku k + F k + v k ; k e )C N 



that is the same as system (|4.9.3|l with K. as the finite set of controlled modes. 

Definition 4.10.2. We shall say that system [ l|4.10.1j) .iV] is time-T solidly 
controllable in observed component if for any uq G V and R > there 
exists a family 

V U0 ,r ■= {v b G L°°([0, T],E K «) | 6 G B U0 , R } 

such that Er(uo, B UQ ^) := Et(wo, V Uo .^) covers Ou(uf°) solidly. Where, by 
y#° we mean the projection of y onto M.^ = O; B Uo ji is an open relatively 
compact subset of a C° -manifold and; Or(jj) is the closed ball 

{xeO \ \\x-y\\ h < R} := {x G M*° \ \\x - y\\ h < R}. 
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Below we will need also open balls we define them by 

R {y) = {xeO\ \\x-y\\ h <R} := {x e R*° \ \\x- y\\ h < R}. 

Proposition 4.10.1. System [ (|4.1U.l|l .l] is time-T solidly controllable in ob- 
served component. 

Remark 18. Proposition \4-.l(F l\ impli es controllability in observed component 
and, it follows from Proposition \4 . 1 fO| ( with N — 1) below. Indeed given R > 
and «o € V, if T < T° it is included in the statement of Proposition \4-.10.2\ 
(with N = 1), otherwise if T > T° we apply any control v G L°°([0, T],R Kl ) 
(for example V = — no control) up to time T — T° arriving to some point 
y G V. Put R := R + \\y Kl — ^g 1 1| . Then apply first part of Proposition \l.lU.2\ 
(with N = 1 and T = T°) to the pair (y,R) G V^xjO, +oo[. The family V y ji o v 
will do. 

Proposition 4.10.2. 

1. For some T° > 0, every < T < T° and every N £ No the system 
[ H4.10.1|) .iV] is time-T solid controllable in observed component; 

2. For each pair (uq, R) G V x [0, +oo[ the family 

Vu ,R '■= { y b\ b G B Uo . R } 
can be chosen satisfying: 

• The map b i— ► Vb is (B, L 2 (0, T, R KN ))- continuous and; 

• The controls Vb(t) are uniformly (w.r.t. b and t) l\-bounded: 

IK(*)k < A = A(T,R 7 u ). 

Since K. 1 is saturating we have O C K, M from some M G N . Fix M with 
this property. We shall prove Proposition 14. 10.^1 in two steps. Prove it in the 
case N > M and prove the "back-induction" step " it holds for N implies it 
holds for N — 1" (JV = 2, . . . , M). These steps are the following subsections 

Emu and mm 

4.10.1 First Step. Proposition 14.1 U.'2l N Big 

In this subsection we shall prove that the statement of Proposition l4~10.2l holds 
for N > M. 

Decompose u G V as u — u KN + U KN where u KN :— P KN u, i.e., u KN is the 
projection of u onto R KN . So U KN = P~ KN u is the projection of u onto (R Kjv )y. 
Now write system [ (14.10. l\ .N] as 

u kn = p k n i_ vAu _ p^Bu + F + v) 
jjkn = p - KN i_ v £_ u _ pV Bu + F + v ^ 
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i.e., 

v£ N = -vAu KN - P KN P v Bu + P KN F + v=: g(u) + v 
Ut N = -vAU KN - P- Kn P w Bu + p- K "F =: G(u). 

Concisely we have 

Ut N = g{u) + v, U t KN = G(u) t £ [0, T]. (4.10.2) 

Let (uq, R) be an element of l^x]0, +oo[. Fix 7 > 1. Let T > be a positive 
real number. 

For each p £ 70_r(O) C R Rn define, on [0, T] the constant control 

v p (t) := T~ x p. 

Since 7 > 1 we have that O R (0) C 70^(0) and, v p (t) dt = f Q T T^pdt = p. 
The family V p :— {v p \ p £ 70^(0)} is parametrized continuously in i 9 -norm 
(q > 0), i.e., the map p 1— > v p is (h, L q (0, T, R#°))-continuous. Indeed for 
q > 0: 

HT-^-T-VllI^ / T T-% 2 - Pl ||« =T 1 -'||p 2 -p 1 ||f i 
Jo 

Hence 

\\t~ 1 p2 - r _1 pi||i9 = T~\\p 2 -pi\\h- 

For g = 00: \\T~ 1 p 2 - T _1 pi||l- = T^\\p 2 -Pilk- 

We can also see that 11^(^)1^! < T~ 1 jR. To prove that the family V p is 
the one we are looking for it remains to check that Et(uo, V p ) covers Or{u^°) 
solidly. 

By the (h, L 2 )-continuity of p 1— > v p and (L 2 , C([0, T], y))-continuity of u ^ 
S s (ito, F, v, v) we conclude the (7i, ^-continuity of 

p Llo o § s (it , F, v p , v)(T) = E T (w , Up). 

Rescaling time: t — T£, £ £ [0, 1]. From (|4.10.2[) we obtain the system 

u* N = T(g(u) + v) U* N =TG{u) u(0) = U , £ £ [0, 1]. (4.10.3) 

which solutions, for v = v p , will be compared with those of the following system: 

y ^ = p = 0; y(0) =y ,££ [0, 1]. (4.10.4) 

Put z = u — y. Then z satisfies 

z« N = T(g(u) +v p )-p = TG(u) £ £ [0, 1], 

i.e., 

-z* N = -vAu KN - P KN P v Bu + P KN F 

rp £ 

—Z* N = -vAU KN - P~ KN P w Bu + p- K «F. 
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which is equivalent to 



-z £ = -vAu - P v Bu + F. 
T ^ 



Multiplying by z we arrive to 

< -f((u, z)) + \(P v Bu, z)\ + \F\\z\ 

< -v\\z\\ 2 -v((y, z)) + \(Bu, z^+C^FM 

and, since 

b(u, u, z) = b(y + z, y+ z, z) = b(y + z, y, z) 
= b{y, y, z) + b(z, y, z) 

we arrive to 

< -HMI 2 - K(y, *)) + c\\ y \\ 2 \\z\\ + c\z\\\z\\\\ y \\ + c^fm 

from which we obtain 

~\z\ 2 + HMI 2 < c 2 \\ y \\ 2 + c 2 \W + c 2 \z\ 2 \\v\\ 2 + c 2 \f\ 2 . 

By Gronwall inequality: 

|z( S )| 2 < exp{TC 2 jT 1 \\y(0\\ 2 ^} (jz(0)| 2 + TC 2 jf * \\y(Q\\ 2 + ||y(0|| 4 + \F\ 
^eMTWlziO)] 2 + TD 2 ) 



2 dt 



where D\ and D 2 depend only on 7, R and 1 1 2/0 1 1 - Indeed y{£) = yo + p£, and 
\\Vo +P£\\ < \\Vo\\ + C|b£||/i and, ||Klk < 1 R - In particular we have (for fixed 
ito, 7 and R): 

Corollary 4.10.3. 



1. Ify = u Q then \u-y\< [Texp(T)]^ 



2. For "bounded" T and yo, say T < T\ and \\yo — uq\\ < (3 we have 



\u-y\< K(\u Q -y \ 2 + l) 2 

with K independent of T (K depends only on 7, R and \\yo\\). 
Note that 

EtK, v p ) = n a o $ T (u , p)(l) 
where $ T (uo, p) is the solution of system l|4.10.3jl . 

Represent by <&°(mo, p) the solution of system (|4. 10.4(1 . Then n o o$°(u 0j p)(l) = 
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Now we compute 

||n c o $ T ( Uo , - n o $> , P )(i)|| loo 
<||n e o $ T ( Uo , p )(i) - n o $°( Uo , p)(i)|| Il 

<C||$ T ( Uo , p)(l) - p)(1)||a < C||* T («o, P) ~ $°K p)llc ( [o,i],H). 

By item Q] of of corollary 14. 10. 31 

|n c o $ T ( Uo , p)(l) - n o $> , p)(l)|| Joo < CX[Texp(T)]^ 
where K is independent of T and p. Thus 

||n c o $ T ( Uo , -)(1) - H o $> , .)(l)llc( 7 o«(0),R^) ^ CX[Texp(T)]5 s 
or, defining in 70^(0): 

G t (p) := ErOcp); G (p) := n o $> , p)(l); 



"0||C( 7 C» B (0),R^) 



< CK[Texp(T)]*. 



Put T° = the unique solution of T° exp(T°) = ( {l ~ 1)R 



VT €]0, T°] ||G T - Mlc( 7 e> R (o),K#°) 



2#OCK 



Then 



< 



2#G 



(4.10.5) 



(4.10.6) 



Note that the map Go is just the restriction of a translation in 



,#e> 



Go(p) 



p — restricted to 70^(0) = C 7 _r(0). By the Degree Theory 9 we have 



that 



p £ a(G o 7fl (o)) = do lR {u*°) 

o, O 7 h(0), p) = deg(I jR , 7 h(O), p - u#°) = ll . (4.10.7) 



Where I jR is the identity function on O lR (0). 

Yet by the Degree Theory we know that for every p £ dO^ R (uf°) and every 
continuous function ip : O lR (0) — ► R# such that 

- ®o\\c { o yR (o),n*°) < \\P- d0 iR( u t°)\\i~ 

degty, O lR {Q), p) = deg(G Q , O- (R (0), p)] 10 (4.10.8) 

We claim that for T < T° we have that G T (O 7fl (0)) covers O r (u*°) solidly. 
Indeed, given T < T° and a continuous function (j> : O jR (Q) — > IR^ such that 



(7~l)fl 



HI C (o tH (0),r# o ) < 2#0 



we have, using H4.10.6(l . 



W ~ G o\\c(o yR {o),v* ) - W- Gt IIc(o 7 „(0),r# o ) + ll G r^ G o|| c(OTR(0): 

( 7 -l)i? , (7-1)^ _ (7-1)^ 



< 



2#0 



2#0 



#0 



8 Where the subscript "oo" means that we are considering the " max -norm ^oc. 
9 See, for example, 6 . 

10 Where \\p, A\\ means the distance from the element p to the set A, i.e., \\p, A\\ := inf{||p — 
a\\ | a G A}. 
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For w £ O r (u* ) we have: 

\\w - dO lR (u*°)\\ loo >^\\w- dO lR {u*°)\\ h > 

so, 

||?!.-Go|| c(0tr(0)iR # 0) || < \\w - dO^Riu* )^. 

By (|4.10.7|) and (j4.10.8jl we conclude that 

deg^, O lR (0), w) = 1. 

which means, in particular, that the equation 4>{y) = w has a solution on 
O yR (0), i.e., 0(O 7fl (O)) covers O r {u*°). 



4.10.2 Second Step. Proposition ^. 10.21 "Back-Induction". 

In this subsection we "imitate" a driving using controls on R KJV by a driving 
using controls on R KN -\ N = 2, ... , M, M is fixed and satisfies O C K M . 
Both drivings leading to the same projection onto IR"™- 1 at final time but, 
possibly going by paths with projections "far from each other" in the middle. 
The projection onto the orthogonal space (K Kjv )y of the paths will be f/-close to 
each other so, at time T the two drivings lead to points close in _ff-metric. Hence 
the end points of the projection onto the finite dimensional observed space O 
are close. Solid controllability will follow from this closeness and (again) from 
a Degree Theory argument. 

Such imitation is then the key for the prove that if the system [ (|4.1U.1|I .N] is 
solid controllable in observed component then so is system [ H4.10.lfl .N-l]. 
After we prove this a N — > N — 1" step it will be clear, from the fact that 
[ H4.10.1|I .M] is solid controllable in observed component (see subsection 14 . 1 . l|l . 
that system (|4.9.3(l is solid controllable in observed component, we just note 
that the systems [ H4.1(J.lfl .l] and (|4.9.3|l are the same system. 

To prove the "back-induction" step "N — > N— 1" we shall need some lemmas: 
For the next Lemma we may consider again the case where the external force 
depend on time: 

Lemma 4.10.4. Given: 



• A finite subset J C N§; J := span{Wk \ k £ J} 

• A function q £ W 1 '°°{[t i , tf], I), such that q(ti) = qt 



• An element Qi £ Sy, the orthogonal space to JJ in V . 

Then there exists a control v J (q,Qi) £ L°°([ti, tt], J) depending on q and Qi 
such that the projection onto J of the solution of the NSE 



-vAu - P v Bu + F + v J (q, Qi), u(U) = qi + Q. L 



equals q on [ti, tf] . 

Moreover the mapv J : (q, Qi) i— ► v J (q, Qi) is (W 1 ' 2 xjy, L 2 (ti, tf, J) -continuous. 
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Proof. Let q, qi, Qi and J be like in the statement of the Lemma. Consider the 
(non controlled) NSE: u t — ~vAu — P v Bu + F with initial condition u(ti) = 
qi + Qi ='■ Ui and split it into 

P J (u t ) = P J {-vAu - P v Bu + F) 
P- J (u t ) = p- J {-vAu - P v Bu + F) 
u{U) = u t 

with the same initial condition, where P J (resp. P~ J ) is the projection H — » JJ 
(resp. H — > If we put u J := P J u and U J := P~ J u we arrive to the 

systems 

f«t = ~vAu J - P J P v Bu + P J F Q 
\u J (U) =qi 



U t J = -vAU J - P- j P v Bu + P- J F 
U J {U) =Q t . 



(4.10.10) 



In the system 1)4.10. 10|) . for each k 6 J replace Uk(t) by qk(t) arriving, in this 
way, to the (closed) system 



U t J = -vAU J - P- J P^B{U J + q) + P- J F Minn, 
U J (t l ) = Q l . 



We can prove existence and uniqueness of a strong solution for this system as 
we prove existence and uniqueness for the "full" equation. We indicate only 
how to find some estimates: Starting from approximate solutions 



U.J..L ._ U k ' L W k 



{U t J ' L , W k ) = - is{AU J ' L , W k ) - {P- J P W B{U J ' L + q), W k ) + {P- J F, W k ) 
= - is{{U J ' L , W k )) - (B{U J ' L + q), W k ) + (F, W k ) (4.10.12) 
Vfce/C L \JJ (4.10.13) 
U J > L {U) = Q1 L ~ J = projection of Q. L onto M Ki - J , 

from which we obtain the ODE 

= -vP^U* - £ u^u J n H(W m ,W n ,W k )-A k 

m,n£}C L 

Uf' L = Q^~ J 

that has a maximal solution defined on \ti,t max [. Now we compute some esti- 
mates that, in particular, imply t max = if. 

Multiplying, for each k, the equation H4.10.12|l by U k ' L and summing up we 
obtain: 

\j t \U J ' L \ 2 ^ -v\\U J ' L \\ 2 + \{B(U J ' L + q), U JX )\ + \F\\U J ' L \ 
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and, after simplify the term with B, using some estimates from section \2. 91 and 
Young inequalities, we obtain 

j t \U J ' L \ 2 + v\\U J > L \\ 2 < C\U J ' L \ 2 \\q\\ 2 + CWqW* + C\F\ 2 
from which we can conclude that for every s G [U, tf]: 

\U J ' L (s)\ 2 <ex P ^ C\\q(t)\\ 2 dt^ 2 +cj*' \\q(t)\\ 4 + \F(t)\ 2 dt 
<£>i(||Qi|| 2 + l) 

where D\ can be taken depending only in ||<7||c([t;, t/], J) ■ So for a constant 
C\ depending only in \\q\\c([u,t f ]j) and we have 

W^L-i^t^^Ci. (4.10.14) 

We also have 

' tf ||^ i (t)|| 2 ^<i?3||^ L |||o o(tj , i/!j x ) (l + £'|k(t)|| 2 dt) 
D 3 r \\q(t)\\ 4 + \F(t)\ 2 dt. 



Hence for a constant Ci depending only in ||g||c([ti, t f ], j) & n d ||Qj|| we have 

\W J ' L \\ LHu ,t f ,^)<C^ (4.10.15) 
From H4.lU.14p and (|4.10.15|) we have that 

{(U J ' L )l remains in a bounded subset ofL°°([fj, tf], H) ^ 
(U j,l )l remains in a bounded subset oiL 2 (ti, tf, V). 

Analogously if we multiply, for each k, the equation H4.10.12f> by —kUj^' L and 
summing up we obtain: 



2 dt 
and, 

d 
dt 



' d \\U J > L \\ 2 < -v\U J >% + \(B(U J > L + q),AU J > L )\ + \F\\U J ' L \ [2] 



\\u^\\ 2 + v\u J -% } < c\u^\ 2 \\u^r + CM[ 2] + C\q\f 2] \\U^\\ 2 + C\FY 



from which, using (|4. 10. 14(1 and (|4. 10. 15(1 . we conclude that for some constants 
C 3 and C 4 depending only in || c([*^, */], J) and \\Qi\Y 

\\U J ' L \\ L ^{u,t f ^)^ C ^ (4.10.17) 
\\U^ L \\ L , {u , tsJ , {A)) <C±. (4.10.18) 



11 F is fixed. 
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The existence of a strong solution follows some classical compactness theorems. 

For the uniqueness we consider the difference w of two solutions V J = 
V J (q, Qi) and U J = (q, Qi) — w := V J - U J . Then from 

w t = -vAw - P- J P V B{V J + q) + P- J P V B{U J + q), 

multiplying by w, using some estimates from section 12.91 and appropriate in- 
equalities, we obtain 



So, \w(s)\ 2 < \w{ti)\ 2 exp fl f C\\U J (t) + q(t)\\ dt = 0. Then a weak solution is 
unique and, so is the strong one. 

We have just proved that the map (q, Qi) t— > U J (q, Qi) is well defined. 
We claim that it is {L 4 (ti, tf, J) x (JTy), X)-continuous, where X is either 
L°°([ti, tf], Sy) or L 2 (ti, tf, §jj,( A } )■ To prove these continuities we proceed 
as usually: Fix a pair (q, Qi) G (W 1 ' 00 x J v ) and consider another one (p, Pi) 
in the same product space. Define w := U J (q, Qi) — U J (p, Pi). Then we obtain 
the equation for w: 



w = -vAw - p- J P v B(U J (q, Qi) + q) + p- J P v B(U J ( Pl Pi) + p). 
To simplify the writing we put Q :— U J (q, Qi) and P := U J (p, Pi) so, 



(U J ' L )l remains in a bounded subset of L°°([i i , tf], V) 
(U j,l )l remains in a bounded subset oiL 2 (ti, tf, D(A)) 



(4.10.19) 



i\w\ 2 + v\\ W \\ 2 <C\ W \ 2 \\U J + q\\ 2 . 



w = -vAw - P- J P V B(Q + q)+ P~ J P V B(P + p) 



(4.10.20) 



and, multiplying by w: 




+ Hp - Q: w )\ + \ \ b (p - q, p, w ) + b (q, v-q, w)\ 



from which we obtain 



j t \w\ 2 + u\\wf < C\w\ 2 \\Q\\ 2 + C\\Q\\ 2 \\p -q\\ 2 + C\w\ 2 \\ P \ 
+ C(\\p\\ 2 + \\q\\ 2 )\\p-q\\ 2 . 



(4.10.21) 



Then by Gronwall Inequality 
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For \\p — <z||l4 < 1 we obtain 12 



(4.10.22) 



and conclude that the map U J is (L 4 x Jp, C([tj, i/], J^))-continuous. 
From l|4.10.21[) we obtain 



2 

II 4 



bill* + 



+ C ||p- g ||| 4 (||Q| 

^ c allHlc([o,r|,ji)+ c, 2lb-?lli*' 
for ||p — q|| L 4 < 1. Hence by <|4.1U.22[1 we arrive to 



Ib-glU* < 1 



l£ 2 (t,,t/,J^) 



^ ^CallQi-Pif + Calb-gl 



(4.10.23) 



and conclude that the map U J is (L 4 x Jy, L 2 (ti, tf, Jy))-continuous. 
Now from equation l|4.1().20)l . multiplying it by Aw, we obtain 

~\\w\\ 3 + v\w\f 2] < \(B(P + p), Aw) (B(Q + q), Aw)\ 

< \b(P, P, Aw) - b(Q, Q, Aw)\ + \b(P, p, Aw) - b{Q, q, Aw)\ 
+ \b(p, P, Aw) - b(q, Q, Aw)\ + \b(p, p, Aw) - b(q, q, Aw)\. (4.10.24) 

Note that 



b(P, P, Aw) - b(Q, Q, Aw) 
b{P, p, Aw) - b(Q, q, Aw) 
b{p, P, Aw) - b(q, Q, Aw) 
b{p, p, Aw) - b(q, q, Aw) 



b(w, w, Aw) — b(w, Q, Aw) — b(Q, w, Aw) 
-b(w, p, Aw) + b(Q, p— q, Aw) 
—b(p, w, Aw) + b(p — q, Q, Aw) 

Hp - 9, P> Aw ) + %> P-Q, Aw )- 



Hence from (14.10.241) . from the estimates of section Inl and from the continuity 
of the inclusions D(A) i— >■ V i— > H we arrive to 

i||| w || 2 + HHf 2] <c(\w\i\\w\\\w\f 2] + \\w\\^\\Q\\\w\f 2] ) 
+ c(\\w\\\p\ [2 ]\w\ [2] + \\Q\\\p-q\ m \w\ m ) + c(\p\ [2] \\w\\\w\ [2] + \p-q\ [2] \\Q\\\w\ [2 {j 

+ C(\p- q\[2]\\p\\\w\[2] + \\q\\\p - q\[2]\w\m / 

By appropriate Young inequalities: 

j t \\wf + v\w\\ 2] 

<D(\w\ 2 \\w\\ 4 + HHIQII 4 ) + D(\\wf\p\f 2] + \\Q\\ 2 \p q\f 2] ] 

+D(\p\f 2] \\w\\ 2 + \p- q\f 2] IIQH 2 ) +D(\p- q\f 2] \\pf + || g || 2 |p - g|? 2] ). 

(4.10.25) 



12 Of course we could ask for ||p — q\\ L & < D for any D > 0. What we need is a first bound 
for \\p\\ L i. 
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Thus 



H 



C{[U,t f ]Jt) 



< exp 



d / \ W (t)nw(t)r + iiqwii 4 + ip(*)if a] dt \\w{u)\ 



+d / \p(t)- q (t)\f 2] (\\Q(t)r + \\ P (t)r + \\ q (t)r)dt 



Then using gH and gj 



lb-g|U* < i 
11^ - Qih± < i 



^([ti.ty],^) 



<^i||Q 1 -P j || 2 + ^i|Ip-9|| 2 L 4 



(4.10.26) 

and conclude that the map U J is (L 4 x (ly), C([U, tf], JJy))-continuous. 
From l|4.1U.25[l we can also obtain 



U 



\w(t)\f 2] dt 



< D 2 \\w\\^„ n „ r±Jl 



tf 



c([o,t],^)v j / i«(*)l 2 IK*)ll a + IIQ(*)ll 4 + b(*)lp]<« 

D»\\p-q\\h(\\Q\\h + \Ml* + \W\\h) 



Hence by (j4.10.22j) and (|4.10.23|) we have 



lb- qWl* < 1 
\\Pi-Qi\U<l 



and, by i|4.1U.26[l we obtain 



\\p-qh* < i 



< D^Q, - P.f + D 4 \\ P - q \\ 2 Li 



^l£ 2 (*i,t/,Ji (A) ) 

(4.10.27) 

and conclude that the map U J is (L 4 x (ly), L 2 (ti 1 tf, Z)(A)))-continuous. 

Now we define another map on the product W 1 ^ {[U, tf], S) x Jy taking 
values on H: 



T J : (q, Qi) >-» -vA(U J (q, Qi) + q) - P V B(U J (q, Qi) + q) + F 



Fix (q, Qi) £ W 1 ' 2 ({t i , tj], J) x JJy and, consider another pair (p, Pi) in the same 
space. Again, as we have done before put Q — U J (q, Qi) and P = U J (p, Pi). 
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We compute the norm of the difference 

\T J (q, Q l )-T J { Pl Pi)\ 
= \vA{P+p) + P v B{P+p) - vA{Q + q)- P V B{Q + q)\ 
<\vA{P + P -Q-q)\ + \P v [B(P +p)- B(Q + q)} | 
<\uA{P +p - Q-q)\ + \B(P + p) - B(Q + q)\ 
<v\P - Q\ [2] + v\p - q\ [2] + \BP - BQ\ 
+ \B{P, p) - B(Q, q)\ + \B(p, P) ~ B(q, Q)\ + \Bp - Bq\ 

MP - Ql[ 2] + Hp - dip] + c\p - Q| [2] (||P|| + HQH) 

+c(|P-Q| [2] |M| + ||P||b-9lM) 

+c(||P|||p-«l[2] + |i > -Q|[a]Nl) 
+C\p-q\ [2] {\\p\\ + h\\). 

Therefore 

\r J ( q , Qi)-r J ( P , p)| 
<v\p - Q\m + v\v - ?|[ 2] + c\\p- q\ [2] (\\p\\ + \\Q\\ + 

+Ci|p-?Ip](||p|| + IIpII + N|) 

<c \p - 0| [2 ](i + IIPII + HQH + h\\) + c \p- q\ m (i + l|P|| + Ibll + ||g||). 

For ||p-g||i4 < 1, and ||P-Q;|| < 1, using l|4.1().26|l we obtain \\P{t)-Q{t)\\ 2 < 
2D i and then ||P(i)|| < \/2L\ + \\Q(t)\\. So we can arrive to 

\T J (q, Qi) - T J (p, P)| < C 2 \P - Q\ [2] +C 2 \p- q\ [2] (l + ||p||) (4.10.28) 

Note that we have used the fact that q £ W 1,oc and so, is bounded. But we 
can not "replace" \\p\\ by a constant in the last member because, unlike as q, p 
is not fixed and we are considering that p varying in L 4 topology. 

/ V-'fo Qi)-r J (p, m 2 dt 

<2Cf J ' |P - Q\f 2] dt + 2C| J S \p - q|f 2] (l + Ibll) 2 dt 

<2Cl\\P - Q\\l 2{t ^ tfiD{A)) + 2C* J*' b-g|f 2] 2(l + ||p|| 2 ) dt 
<C4P- Q\\l* {tutJtD{A)) + C 4 \\p- q\\ 2 Li (l + |b||| 4 ). 
Therefore, using 1(4.10. 27|) : 

f lb-g|U*<i 
\ IIP - Qih$ < 1 

\\T J (q, Qi) ~ T J (p, Pi)\\l Hti!tfiH) dt < C 5 ||P - Q t \\ 2 + C B \\p-q\\l< 

(4.10.29) 
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and we conclude the (L 4 x (Jy), L 2 (ti, tf, i/))-continuity of r J . 
Now we can indicate which is the control v appearing in the statement of 
the proposition: In fact 

v J (q, Qi) :=q-P J T J (q, Qi) 

satisfies the statement. Indeed its (W 1 ' 2 {[t i , tf], J)x(Jy), L 2 (ti, tf, J)-continuity 
follows from the (L 4 x (JJy), L 2 (ti, tf, i/)-continuity of r J and from the 
(W 1,2 ([ti, tf], J), L 2 {U, tf, J)-continuity of q i-> q. 

To prove that the projection of the solution of the system 

u t = -vAu-Bu + F + v J {q,Qi), u t = q(U) + Q 4 (4.10.30) 

coincides with q we differentiate q + U J (q, Qi) obtaining 

[q + U J (q, Qi)] t = q - vA{U J ' (q, Qi)) - p- J P V B(U J (q, Qi) + q) + P~ J F 
= -vA(q + U J (q, Qi)) - P V B(q + U J (q, Qi)) + F 
+ q- {-vAq - P J P V B{q + U J (q, Qi)) + P J F) 
= -vA{q + U J (q, Qi)) - P V B(q + U J (q, Qi)) +F + v J (q, Qi) 

showing that q + U J (q, Qi) is the (unique) solution of (|4. 10. 30(1 . 

To finish the prove remains to verify that v J £ L°°([tj, tf], J). Since q £ 
I4 7l ' 0O ([t 4 , tf], J) we have q S L°°([ti, tf], J) and by 

||r J (g, Qi)\\ v > < v\\Q + q\\ + \\Q + qf + \\F\\v> 

< C(\\Q + q\\ + \\Q + q\\ 2 + |F|) < Ci. (4.10.31) 

Hence T J (q, Qi) £ L°°([ti, t f ], V) and then, P J T J {q, Q t ) £ L°°{[t t , t f ], P J V), 
i.e.,P J T J (q, Q l )£L™({t t ,tf],J). 13 
Moreover we can see that 

\\v J (q, Qi)\\L°°([ti,t f ],j) <ll9lU°°([ti,t/],jr) + ll rJ (<7, Qi)\\L°°([ti,t f ],g) 

<ll9lli«>([ti,t/],j) +Di\\T J (q, Qi) || ioo([ ti , t/ ], v) 

and, by 14.10.31|l . we obtain 

\\v J (q, Qi)h-au,t f ],I) <D 2 (4.10.32) 

where D 2 depends only on the norms J q\\ L <*> ([ tl ,t f ], J) i h\\ L°°([u,t f ]j) and 
ll < 2II.L o °([i 4 ,i / ],j i )' Then using (|4. 10. 17(1 the constant D 2 can be chosen depend- 
ing only on ||g||vK 1 .°°([t*,t / ],J) and HQ* II- a 

Definition 4.10.3. We call (5-metric the function defined on the product space 
(L°°([0, T],R d )) 2 by 

S(u, v) := meas{t £ [0, T] | u(t) ^ v(t)}. 

13 By the equivalence of V-norm and Zi-norm in J. One can check that V' coincides with the 
domain D(A~ 2 ) of the operator A~ 2 and its Fourier characterization is V = {X]fcgN 2 u kWk I 
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Remark 19. The (5-metric is the restriction, to the space of ordinary controls, 
of the strong metric defined on the space of relaxed controls (see J^). 

In [7] (Chapter 3) we can find the so called Approximation Lemma that 
says that a strongly continuous family of relaxed controls can be weakly ap- 
proximated with arbitrary accuracy by a strongly continuous family of ordinary 
controls. At the end of that chapter we can find a "Remark on the Terminology" 
leading us to the following Lemma: 

Lemma 4.10.5 (Approximation Lemma; |7J). Let A CM. d be the convexi- 
fication of a finite set of points: 

A := Conv{pi, p 2 , p r } 

and, V := {v(t, b) £ L°°([0, T], A) \ b £ B} be a L 1 -continuous family of 
A-valued functions. Then for each e > one can construct a family V e :— 
{v e (t, b) £ L°°([0, T], {p u p r }) | b e B} of {p u p 2 , Pr}-valued func- 
tions such that 

• V e is 8-continuous, i.e., b i— > v £ (-, b) is (B, 5) -continuous; 

• V £ e- approximates, uniformly w.r.t. b, the family V in relaxation metric, 
i.e., V6 € B \\v e (-, b) — v(-, b)\\ rx < £ and; 

• The elements of V 6 are piecewise constant and the number of intervals of 
constancy is the same for all b £ B. 

Remark 20. Note that 5-metric and L q -metric (1 < q < +oo) are equivalent 
in the subset of piecewise constant functions taking values on a fixed finite set. 
Also, since the controls take values on a finite set, the S-continuity of V £ is 
equivallent to the continuity of the lengths of the intervals of constancy of the 
controls. 

In Lemma 14.10.51 is said that the intervals of constancy can be taken the 
same for all b S B but, looking at [?.; ch. 3] some of those intervals may 
degenerate to a single point. We claim that we can suppose non-degeneracy of 
the intervals, 14 . We may even suppose that there exists a lower bound 6 s for 
the lengths of the intervals of constancy of the family V e , i.e., for all b S B none 
v £ (-, b) has an interval of constancy with length less than 9 e . 

Looking at [0; ch. 3] we see that the intervals of constancy are (or can be) 
constructed in the following way: First we subdivide the interval [0, T] into n 2 
intervals — Li, i = 1, . . . , n 2 — with the same length — 5-. Then subdivide 
each one of these intervals Li into r subintervals — Ly j = 1, . . . , r — which 
lengths length(Lij) = J L( Vj(r, b)dr, v(-, b) = J2j=i v j(^ b )Pji v j(^ b ) G I ' l]i 
depend (continuously) on b and the interval L^ x lays on the left of Ly 2 if j\ < ji . 
To this partition is associated the piecewise constant control v„ 2 ('i b) defined 
by: 

t e Lij v n 2(t, b)=pj\ i = 1, n 2 , j = 1, r, 

i.e., in each interval Li we use all the controls from {pi, . . . , p r } using pj t before 
Pj 2 if ji < ]2- Note that as we have said before some control may be used for 
time zero. 



1 Note that it is not enough to eliminate the degenerate intervals because the number of 
intervals would not be the same for all b £ B. 
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Each v e (-, b) G V s have the form v e = v n 2 for some n depending only on e. 
We shall need the following strong result without non degeneracy of the 
intervals of constancy: 

Corollary 4.10.6 (Approximation Corollary). Let A C M. d be the convexi- 
ftcation of a finite set of points: 

A := Conv{pi, p 2 , p r } 

and, V := {v(t, b) G £°°([0, T], A) \ b G B} be a L 1 -continuous family of A- 
valued functions. Then for each e > there is 9 £ > and a family Z e := 
{z e (t, b) G L°°([0, T], {p u P2, Pr}) | b e B} of {pi, pz, . . . , pJ-waJued 
functions such that 

• Z £ is S -continuous; 

• Z £ e- approximates, uniformly w.r.t. b, the family V in relaxation metric, 
i.e., V6 G B \\z E (-, b) - v{; b)\\ rx < e; 

• The elements of Z 6 are piecewise constant and the number of intervals of 
constancy is the same for all b £ B and, 

• For all b G B all the intervals of constancy of z £ (-, b) have a length not 
less than 9 e > 0. 

Before the proof consider the following Lemma: 

Lemma 4.10.7. Given K > 0, 7 > and L G No- Define the sets 

L 



Vo := {(x u x L ) G R L I x t > 0, J2 x i = R } 

»=i 

L 

Ve := {(xi, . . . , x L ) G M L | Xi > 6, ^ Xi = K} 



i=l 



where 9 > 0. 

Choose n G No swc/i t/ia< < 7 and pwt 



nL 

„ K 



15 

nL 



(4.10.33) 



Then the map Po.e = {Pq£, ■ ■ ■ , Pqo)> defined on Vo by: 

is continuous, take its values on Ve and, satisfies \Pq e (x) — Xi\ < 7. 

The Approximation Corollary follows from the Approximation Lemma and 
from the last Lemma. Indeed let e > be a positive real number, put D := 
m a x ie{i,...,r}{lbi||} 5 ta kc K = 1, L = r, 7 < 57%?! G N such that £±± < 7 

15 So, 6 depends on both K, L and 7. 
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and 9 = -k 



The continuity of Po,e implies the continuity of the map 

r r 

v(t, b) = ^ Vj (t, b)p 3 ~ z(t, b) =Y, P oMt> b )) 



3=1 



3 = 1 



and then the family of the controls z(-. b) — ■_ 1 P^ g {v{t, b)) is parametrized 
in L (0, T, Conv{pi, . . . , p r })-norm. 

The family iJ |-approximates the family V in relaxation metric. Indeed 



\\z(;b)-v(;b)l 



< \z(r,b)- v(t, b)\dr = \^ Zj (t, b)pj - J2 *i (t, b)p 3 1 dr 
Jo Jo j=1 j=1 

< f \Y,\z J (t,b)-v ] (t,b)\Ddr<T 1 Dr< £ - 



3=1 

Now apply Approximation lemma to the family Z and find a family Zi that 
|-approximates Z in relaxation metric. Hence 

\\Z%(; b) - v(; b)\\ rx < \\Z%(; b) - z(; b)\\ rx + \\z(; b) - v(; b)\\ rx < S. 

Therefore the family e-approximates V in relaxation metric. 

From the Approximation Lemma the number of intervals of constancy of 
zi(-, b) is the same for all b G B and, the Lij interval of zi(-, b) has length 
f L zj(t, b)dr > ^9 for that n G No such that zi(-, b) = z n 2(- : b). Thus 
all intervals have a length not smaller than ^9 > — a positive constant 
depending only in e. 

Proof of Lemma \4-10. 1\ The continuity of each Pq g , and then of Po,e, is clear. 
From 

L L 



E^» = E[( 1 -^)(^-f) + f 



= K 



and 



1 



1. 



K K 
nL ~ nL 



we conclude that Po,e takes its values on Te- 
xt remais to estimate \P l g (x) — Xi\: 



Xi K 
n nL 



K K (L + 1)K 

< 1 < — < 7- 

n nL nL 



□ 



Now, we are ready to start the proof of the induction step. Fix uq G V. By 
"back-induction" hypothesis system [ l|4.1().l|I .N] is time-T solid controllable in 
observed component so, there is a family V := {v(-, b) G L°°([0, T], M KJV ) | b G 
B} parametrized continuously in L 2 -metric such that 

E T (u , V) covers R (uf°) solidly. (4.10.34) 
\\v(t, b)\\ h < S V(t, b) G [0, T] x B. (4.10.35) 
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Note that, by (|4.10.35|) we obtain 

V(t, b) G [0, r]xB »(t, 6) G CW{±S e fc | fc G K N } 16 
Then there exists S > such that 

V(i, 6) G [0, T] x B 

v(t, b) G C*onw{±Sefc, ±S<5 m ,„ | fc G /C^ -1 , (m, n) G SW-i} =: Cat (4.10.36) 

Now fix e > and e-approximate the family V taking values on Cn by a 
family Z e := {z £ (t, b) G| & G -B} taking values on {±Se/c, ±S<5 m) „ | A: G 
/C^- 1 , (m. n) G gjy-i} like in Corollarv l4.10.fil 17 
Choose a real number 7 such that 



0< 7 <||«f 



Remark 21. It is clear that \\u* - dE T (u , V)^ > R because of l|4.10.34j) . 

To see that the inequality is, in fact, strict we suppose not and then, there is 
p G <9Et(uo, V) such that Wu^ — p\\i ± = R. For the sequence of continuous 
functions 4> n defined in B by 



we have 



b n - E 2 



1 



1 



E T (u , v(; b)) - itj 



#0 



\C{B, 



n 



-IIErll C (B,i 



<t>n{V) = (l - -)E T K, V) + -U*° 
V 71/ rt 



1 - 



,#0 



t -(i-i)b-^°l h 



n 

< R. 



Therefore 



-)p+-u*° edMV)no R (u*o) 



and then, for all n G No, (f> n (V) does not cover On(uf°) which contradicts 

Now note that for x G R (uf°) we obtain \\uf° - <9E T (w , V)\\ h < R + 
\\x — <9Et(mo, i.e., 7 < \\x — dEr(u , V)]^ and then 



#0 



< \\x-dE T {u , V)\\ loB . 



(4.10.37) 



By the continuity of S srx there is e > such that 



\v - w\\ rx <e^> 

\\TLo o S.(« , F, w, v)(T) ~ n a o S.(« , F, v, v)(T)\\ 



C(B,R*° < 2 #0~ 



16 Note that the norm used in 14.10.351 is the l\ one. 

17 Recall that Sm-1 has been defined, in the proof of Proposition 13. 5."3l for N > 2 and in 
the proof of Proposition 13. 5,^1 for TV = 2, as the set of "extracted" pairs of modes associated 
with the new directions. 
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Fix e < e and put 

E e T {u Q , z e (b)) := n o § s ( Mo , F, z £ (6), v)(T) 
and consider </> defined on B such that 

||0(O - E§,(«o, « e (0)llo(B lH #« < 

Hence 

110 - e t\\ C (b,m*° <U~ e t\\ C (b^*° + H e t - e t|| c(B:R #« < ^ 
and, by l|4.10.37[l we conclude 

Vx G R (ufO) deg(<p, B, x) = cfeg(E T , S, x) ^ 

so, 0(B) covers R (ufO). Therefore E^(u , Z £ ) covers R (ufO) solidly for all 
positive e < e. Moreover, for all pairs (t, b) we have ||z e (£, < 3- Therefore 
for e < E system [ l|4.1().l|I .N] is is solid controllable by means of the family Z e . 

Now we fix £ g]0, e[ and define Z := Z € ; z(t, b) :— z £ (t, b). We shall com- 
pare the trajectories generated by a control z(t, b) € Z with another generated 
by some "appropriated" control taking values on R' 1 ™- 1 . 

To the control z(-, b) is associated a partition X(b) of [0, T] into m non- 
degenerated intervals of lengths Xi > 8 > 0: 

m 

X(b) = {x u x 2 , x m ) G M m , me No, = T, 

i=i 

where m and are independent of the parameter 6. We put 

A(b) = {(0 = a , ai, . . . , a m = T)} 
for the end points of the intervals in X(b). So, 

A(b) G A e := {(a , a u . . . , a m ) G M m+1 | a = 0, a m = T, 

m 

a t - a t -i > 8, )Xaj - = 
i=l 

Let w <E K, w > 3. Now in [0, T] we want to define a function </v,(-, b) associ- 
ated to X(b) with the following properties 

• <p w (-, b) vanishes at the points en,, i = 0, .. . , m; 

• ^(-j 6) G W^°°([0, T], E) with 

w(l +0) 



l!<M-, &)llc([0,T],R) < 1; \\4>w(-, &)||l«([0,T],P^ ' 

5(</> w (t, b), sin(wi)) < ^ and; 
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• For fixed w, the map A(b) <j) w (-, b) is (A$, W 1 ' 2 (0, T, M))-continuous 
(where As is endowed with the topology induced by M m+1 ). 

We proceed with the construction of <f> w (-, b): For each i e {1, . . . , m} put 
Pi = 21 . Then subdivide each interval [ctj-i, ai\ into 

[ai_i, a { ] = [<Xi-i, cti-i + pi[l)[ai-i + pi, ai - p i ]U]a l - p u on}. 18 

For each i = 1, . . . , m put 

M^-r+ Pi )) {t _ a ._ l) if te [ a ._ 1 ,a._ 1+ p.] ; 

b) = ( sin(tot) if i e + p,, a { - p^; 



-Pi 



(t - a^ if t e [a* - aj 



Then the graph of the restriction of <f> w (-, b) to an interval [a,_i, Qj] is a 
concatenation of a straight line, a piece of the graph of sm(wt) and another 
straight line. From the construction is clear that 4> w (-, b) vanishes at the points 
ai, i = 0, . . . , m and that <j) w {-, b) is continuous with ||</>k,(-, &)||c([o,t],ir) < 1- 
In the subintervals ]ai_i+pi, a^ p,[ we have ^(t, 6) = wcos(wt) so,\<p w (t, b)\ < 
w < w ( 1 f f e 1 . In the subintervals ]a,_i, + pi[ and ]«j — pi, ati[ we have 

IM*. 6)1 < j- = % < f < Hence we have II^O. 6)IUoc ([0 ,ri,R) < 

22^£i. Therefore ^(-, 6) e W 1 > oo ([0, T], R) and 

w(l + 0) 

||<Pu>('> D)||w 1 .°°([0,T],R) < H q ■ 

We see that <f> w (t, b) differs from sin(wi) only in the intervals [aj_i, aj_i + p,[ 
and ]a>i- pi, a>i] so, 

m m 2 y 

<5(^(t, 6), sin(w;i)) =V2 ft = V2^<-. 

i=l i=l 

It remains to check the continuity property. For that fix w > 3 and A e Ae- 
Let 7 > and B £ Ae such that B]]^ < 7. Let A = (a , oti, ■ ■ . , a m ) and 
B = (fa, fa, ■ ■ ■ , fan) be the coordinates of A and B. Then we have |aj— < 7 
for « = 1, . . . , m — 1, a = fa = and a m = /3 m = T. For small 7 (7 < £) 19 , 
putting ^ and 4>b for the functions <p w associated with A and B, we have that 
</iu and 4>b differ only in the the following union of subintervals 

[0, min{pf, pf }] U [min{pf, pf }, max{pf, pf }] 
U[min{a 4 - pf , - pf }, max{a, - pf , fa- pf}] 
U[max{ai - pf, fa - pf}, min{aj, fa}] 
U[min{ai, fa}, max{a i; fa}] 
U[max{a 4 , fa}, min{a 4 + p^ 1; & + pf +l }\ 
U[min{a 4 + pf +1 , fa + pf +l }, max{a 4 + pf +l , fa + pf +1 }\ 
U[min{T-p£, T-pf}, max{T - p£, T - pf}} 
U[max{T - pi , T - pf}, T]. (4.10.38) 



18 



Note that, since w > 3 we have pi < ^7- and the subdivision is well defined 



19 Note that for 7 < ^ we have ~( < pf and y < pf for all i 6 {1, ... , m}. 
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Now we prove that \\4>a — 4>b\\c([o,t],r) goes to zero as 7 does. For that is 
enough to prove that \\4>a — 4>b\\c(i,vl) goes to zero as 7 does for every interval 
/ of the union H4.10.38fl . 



For / = [0, min{pf, pf } 

\\(f>A - 4>b\\c(I 

sm(wpf) 



max{ 
te/ 



s'm(wpf) sin(u>pf ) 



pf 



m inK,pf}-^flminK,pf} 
Pi 



• For / = [min{pf , pf }, max{pf , pf }]: 

- <?M|c(/,R) <(max{pf , pf } - min{pf\ pf })||^ - <Mc([o,t],r) 
<2|pf -pf| =2-1x1-^1 <2^ 

• For i = 1, . . . , to — 1 and 

- For / = [min{a l - pf, ft - pf }, max{a l - pf, ft - pf}}: 

\\<f>A - <t>B\\c(i,M) 

<(max{a, - pf, ft - pf } - min}^ - pf, ft - pf})\\(/)A - <I>b\\c° 
<2\(a t - pf) - (ft - pf )| < 2K - A I + 2|pf - pf I 

<2 7 + 2-|x ( - < 2 7 + 2^ = 2 7 (1 + -). 
w w w 

- For / = [max{ai - pf, ft - pf }, min}^, ft}]: 

<Pb\\c(i,W) 

sin(w(ai - pf) sin(ui(ft - pf )) 



: max^ 
te/ I 



= max 



-Pt 
sin(w(ai - pf) 



-Pi 



(maxja, - pf , ft - pf } - 



sin(w(ft - pf )) 



(max( a ! - pf. A - pf } - Pi) 



sin(w(aj - pf) 



-(min{a 4 , ft} - on) 



sin( W (ft - pf)) 



(min{a i; ft} - ft) 



For / = [min}^, ft}, max{a„ ft}]: 

- <Mlc(/,R) < (max{ai, ft} - min{aj, 0i})\\4>A - 4>b\\co < 27. 



20 Note that the maximum is attained at one of the end points of the interval I, because the 
functions (f>j\ and cf>B are affme in /. 



116 CHAPTER 4. CONTROLLABILITY IN OBSERVED COMPONENT. 



- For I = [max{ai, ft}, min{a l + pf +1 , ft 4- pf +1 }\. 



\\4>A~ 4>b\\c(i, 



■ max< 
tei I 



= max 



\sm(w(ai +pf + i) 



i sin(w(a, + pf +1 ) 
I pf 



(t ~ Q.) 



shM>(ft +pf+i)) 



(t- A)|} 



(max{a„ ft} - a^) 



sinMft + p? +1 )) 



(max{a„ ft} - ft) 



sin(u>(a 4 + P-L) 



sinWft_+P, B +i)) 
pf 



(min{ai + pf +1 , ft + pf +l ) - a,) 

(min{ai + pf +1 , fa + pf +1 } - fa) 



- For J = [minjcti + p^ x , ft + pf +x }, maxfa, + pf +1 , ft + pf +1 }\. 

\\4>A - 4>b\\c(I,Wl) 

<(max{a 4 + pf +1 , ft + pf +1 ) - nim-ja, + pf +1 , ft + pf +l })\\4>A - </>s|| 
<2|(o i + pf +1 ) - (ft + pf +1 )| < 2|a« - ft| + 2|pf +1 - pf +1 | 

<2 7 + 2-|a? <+ i - < 2 7 + 2^ = 2 7 (1 + -). 

W WW 

For I = [min{T - p£, T - p£}, max{T - p£, T - p*}]: 

\\<pA - <I>b\\c(I,R) 

<(max{T - pf n , T - pf n } - min{T - pfo, T - pf n })\\4>A - <Mc([o,T],«) 

<2|p£-p£| = -\ Xm -y m \ < 

w w 

For I = [ma^{T- pi, T-pg},T]: 

\\<pA ~ <Ab||c(7,M) 



= max] 



sin( W (T - p£) 



P» 



(t-T) 



sm(w(T-pg)) 



(t~T)\} 



sin(w(T - p£) 



(max{T-p£,T-p£}-T) 



sin( W (T - p£)) 



(max{T-p£,T-p£}-T) 



Therefore since ft — > and p^ — > p^, as 7 — ► we have that for every interval 
I in the union 1|4.10.38|) 



Hence 



\\<t>A ~ <t>B\\c(i,W) as 7->0. 
0b||c(o,t,m) ->0 as 



(4.10.39) 
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For the derivatives we have 



\<j> A (t) - <j> B (t)\ 2 dt 



r iaia{pf,pf } <.max{pf,pf} 

-' - ' \<p A (t)-<p B (t)\ 2 dt+ / \Mt) 

Jmin{ P f ,pf } 



/() 

m — 1 



i B (t)| 2 ^ 



£ { Qi _pf )ft _pf} 



min{a;-pf, pi -pf} 



^t) ~ B (t)\ 2 dt 



min{a.i , (Si} 



|<^(*)-<M0l 2 * 



|^(i)-0 B (t)| 2 dt 

min{a,+pf +1 ,/3 l +pf +1 } 

max{a^, /3i} 
max{a j +p^ hl ,ft+p? hl } 



b A {t) - j) B (t)\ 2 dt 



/.max{T-p^,T-p^} 

+ / \4> A {t) - <M*)| 2 dt 

Jmin{T-p4,T-pB} 



<{T-p^,T-p^} 



b A {t) - j> B {t)\ 2 dt. 



(4.10.40) 

(4.10.41) 



Now we prove that each one of the intervals on the right-hand side of (|4. 10.41(1 
goes to zero when 7 does: 



L 



rain{ Pl , p 1 } 



\<j> A (t) - <j> B (t)\ 2 dt 



< 



sin(wp^) sin(wpf ) 



pf 



pf 



min{p! , p Y } 



/■max{pf , pf } 

/ \Mt) -<p B {t)\ 2 dt 

Jmin{pf,pf } 

<(max{pf , pf} - min{pf , p?})\\<j> A - 0s||i=( O ,T,R) 

=\p? - p?\{2^-^) =(2^-^) -\*i-*\<*>n{— ) 



Recall that out of the intervals of I4.1U.38I 4>a an d 4>b coincide. 
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For i = 1, m—1: 

pma,x{eti— pf, @i—pf} 

/ \Mt) - j> B (t)\ 2 dt 

J min{ai — pf ,Pi—pf} 



<(max{a, - Pi , Pi - Pi } - min^ - p i , Pi - p t })\\(j> A - <MIl°°(o,t,r) 

KIT) 



- 1 i - - . ' 



- (a - pf )i < ^{^-)\\ ai -Pi\+\ P f- P f\) 

Kl)s w(Y) 2 (, + 5) 



min{ai, ft} 



|^(t)-^ B (<)| 2 dt 



max{a»-pf , /3,-pf } 

sm(w(ai - pf )) _ sinQ(/3 t - pf )) 

-Pf 'Pf 

■ (min{a;, /3J - max{aj - , $ - pf }); 



6 A (t) -0s(t)| 2 dt 



<(max{a i; pi} - mm{a l7 Pi})\\<j> A - <Mll°°(o,r,R) < Aw 



min{a i +pf +1 , ft+pf +1 } 

|(/.A(t)-0B(i)| 2 * 

max{o:i , /3i } 

sm(w( ai + pf +1 )) sm{w{Pi + pf +1 )) 



Pi+i 



Pi+1 



■ (mm{ ai + pf +1 , Pi + pf +1 } - max{a l7 Pi}); 



L 



\j> A (t) - <M*)| 2 dt 



lmin{a i +pf +1 ,0 i +pf +1 } 

<(max{a 4 + pf +1 , fa + pf +l } - min^ + pf +1 , fa + pf +l })\\4>A - </>b||1oo (0iTj 



<Aw z 



a t + p? +1 )-(fa + pf +1 )\<4W( 



i+0\ 2 



l-max{T-p^,T-p^} 

/ \<P A (t) - <P B (t)\ 2 dt 

•/min{T-p£,T-pB} 



<(max{T - p£, T - p*} - min{T - p*, T - p^})\\4>A - 0b|Il~ (o ,t,i 



(i + Oy. 



m A^m I 



Aw' 



i+e\ 2 
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Finally for the last integral of ((4.10.41(1 we have 

/ \Mt) - Mt)\ 2 dt 

7max{T-p^,T-pS } 

sm(w(T-pA)) sm(w(T-p*)) " 



-p A -p B 

rrn rrn 



(T - max{T - p A , T - p*}). 



Since — ► en and p^ — > as 7 — » we have that all the intervals on the 
right-hand side of ((4.10.41(1 go to zero as 7 does, i.e., 

\\<t>A-M\ L *{Q,T,n) ^0 as 5->A (4.10.42) 

Hence from (|4.10.39() and ()4.10.42|) follows the (A, W^ 2 (0, T, R))-continuity 
of A 1—* 0^. 

Now from the (B, ^)-continuity of the map b t— > A{b) (which is equivalent 
to the (5-continuity of the family Z) and, from the (Ag, l / F 1,2 )-continuity of 
A 1 — ► 4>a we have the following corollary 

Corollary 4.10.8. The map b i-> <j> w {-, b) is (B, W 1 ' 2 ^, T, R))- continuous. 
Imitation. 

Now we "imitate" the control z(-, b) d Z taking values in {±Sefc, ±S5 m! „ | k G 
/C^ -1 , (to, n) G SV-i} by a control z w (-, 6) taking values in R KJV - 1 . 
Take the solution u°°(-, 6) of the equation 

u~(., 6) = -iMw 00 - Bm 00 + F + 2f(., 6), u(0) = u 

and, consider its projection onto R''"- 1 : 

q°°(., b) = P KN - 1 u°°(-, 6). 

Let {0 = ao < «i < ■ • • < ot m = T} be the end-points of the intervals of 
constancy of z(-, b). For w > 3 define the control z tu (-, 6) by recursion in the 
following way: 

• In the first interval of constancy [ao, a{\: 

z{; b) if z(; b) G {±Se fc I fc G JC N - 1 }- 
:•■■[■. b): { v K --Hqr(;b) + V^r(;b)(e m ±e n ),Uo) 22 
if z(-, b) G {±5<5 mj „ I (m, n) G Sjv-i}. 

where f/ is the projection of uq onto JTy = (R Kw )y — the orthogonal 
space to R KJV ^ 1 in V and, q^°(-, b) is the restriction of q°°{-, b) to [ao, a{\\ 

If the control z w (-, b) is already defined in the first p — 1 intervals of 
constancy (up to a p _i), we define it in the p th interval a p ] by: 



if z{; b) € {±^<5 mi „ I (to, n) G SV-i}- 







z™(-, 6) := < 





22 Herc is the control given by Lemma [4.1U.4l for J = K N 1 . 
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where U w is the projection onto (R Kjv )y of the solution of the equation 
<(•, b) = -vAu w - P v Bu w + F + z w (-, b), u(0) = u , t G [0, a p _i] 
and, q™(-, 6) is the restriction of q°°(-, b) to [a p _i, a p ]. 

We shall prove that at time T, u w {T) goes, uniformly w.r.t. 6, to u°°(T) in 
L 2 -norm as w; goes to oo, i.e., 

Lemma 4.10.9. For any e > £/iere exists w £ > 3 swc/i i/ia£ 

VieBVw>w e |u w (T, 6) -u°°(T, 6)| < e. 

We claim that if the statement of the previous Lemma is true and, if we put 
7 = |M|9Et(woj Z), uo||/i — Rj then, Et(uo ; -Z™ 6 ) with e = ^oc anc ^ ^ ' s a 
constant such that ||ar[|; 1 < C||x||/f (x G R# ), covers On(u n ) solidly. Indeed, 
let cj) be a continuous function defined on the closure of -B such that 

U-Eriuo, Z^)\\ ciBtRS a } < 

then 

||0-E T (w o , Z)|| C(B R # 0) 
<||^-Er(uo, ZW ')Wc(b,r*°) + HM^o, Z^O-EtCuo. z )\\c(bm*°) 

<ik +C£<2 ik- 
For a; G Or(uq) we obtain 

||0-E(u o , <^Pt(«o, 2), uollh -i?< ^I|Et(uo, Z), 

<||Er(«o, 2), z|| ioo . 

Hence deg(4>, B, x) = deg(Ex, B, x) ^ and so, <fr(B) coversO k(uo) ■ 

Therefore what remains to conclude the proof of the back-induction step and 
then, the solid controllability in Observed Component of system 1(4. 9.3(1 is the 
proof of Lemma 14 . 1 . 91 that we present in the next section. 

4.11 Proof of Lemma f4. 10. 9L 

First we note that at the times a, we have that the projections of u°°(c<i, b) 
and u w (a.i, b) onto l"™ -1 coincide. In fact, denoting the projection of u w (-, b) 
by q w {-, b) we have 

?<"(•, ft) = ?°°(-, ft) + £ CUM) 

(m, n)6Sjv-i 



where 



if v(t, b) G {±Se fc | k G /C^- 1 } 

V2Ecl) w (t, b)(e m ± e„) if z(t, 6) G {±S5 m , n | (m, n) G Sjv_i}. 
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Since <j) w (-, b) vanishes at the points a, also 0™ „(i, b) does. Hence g°°(a,, b) = 
q w (ai, b) for all i = 0, 1, . . . , m. In particular for i = m we obtain q°°{T, b) = 
q w (T, b) and so, we need only to compare the projections onto (R K ™- 1 )^. For 
that we shall need some lemmas: 

Lemma 4.11.1. The solution of the controlled NSE 

u t = —vAu — P v Bu + F + v, u(0) = u o , t <G [0, T] 

satisfies 

II m *IIl 2 (o,t, h) < c v 

where C v is a constant depending only on the norm |H|l 2 (o.t, #)- 23 

Moreover for a given family of controls {vp | (3 G B} bounded in L 2 (0, T, H)- 
norm, i.e., if there is a constant K > such that for all P € B H^Hl^o, t, h) < 
K, then we can find a constant Ck depending only on K such that 

\\uf\\mo.T,H) < C K 

where u 13 is the solution of 

uf = -vAu 13 -P v Bu p + F + Vf3 , u(0) = u , t e [0, T]. 
Proof. Multiplying the equation by u t we obtain 

-"^ ll?i||2 + q|u|l|u| [ 2 ] KI + |F|KI + |l;||Ut| 



so, 



\Wt\ 2 < -^j t \H\ 2 + Cill«ll 2 l«lf 2 ] + Ci|f| 2 + c\\v\ 2 

f T \ut(t)\ 2 dt<C 2 \\u\\l a ^ Thv) (l+ [ T \u(t)\f 2] dt) (4.11.1) 
Jo Jo 

+ 2C X T\F\ 2 + 2Ci / \v\ 2 dt. (4.11.2) 
Jo 

Multiplying the equation by u and Au, analogously as we have obtained the 
"a priori" estimates in the proofs of existence of weak and strong solutions, we 
obtain the inequalities 



IMIc([o,T],ff) < Kl +C\\\F + v\\ L2(0 T V , ) 
\\ u \\l 2 (0,T,V) < 2\\ u \\c([o,t],h) +Civ\\F + U||^ 2 ( 0)T) ^) 
\\u\\c([o,t],v) < exp(C 2 ||u||c([o,T],H)l|w||L 2 (o,T,y))(l|wo|| 2 
+ C 2 \\F + v\\ 2 LHQTH) ^ 
\ u \\l 2 (o : t,d(a)) < 2||u|| c .([ 0i T],y) + C3||w|Ic([o,t],h)IMIc([o,t], v) 

L 2 (0, T, H) ■ 



C 3 \\F + v\\ 2 L2( , T „ y (4.11.3) 



23 Recall that we have fixed u, F, T and uq. Otherwise the constant would depend on them. 



122 CHAPTER 4. CONTROLLABILITY IN OBSERVED COMPONENT. 



from which, looking at 1)4.11.2(1 . we conclude that ||itt||£2(o,T, H) < C v f° r some 
constant C v depending only on ||u||i2(o,T, H)- 

The statement relative to the family vp is also clearly true from the previous 
expressions: 



In each one of the previous estimates if we "replace ||u^||ia(o t, h) 

by the bound K" we obtain bounds independent of (3. (4.11.4) 

□ 

Corollary 4.11.2. There is a constant Ceo > such that 

VbeB \\u?(;b)\\ LHOtTiH) <C X . 
Proof. For all t S [0, T], b G B we have 

IK*, b)\\ h < SmaxjllxH^ | x € {e fe , <5 m ,„ | k e /C^ -1 , (to, n) e Sjv-i} | 
so, we have that for some constant Ki 

IK-> 6 )IIl-(o : t,h) < K 2 . 

The result follows by the last part of Lemma 14.11.11 and the continuity of the 
inclusion L°°(0, T, H) — > L 2 (0, T, H). □ 

If v(t 7 b) £ {±Sefc | k e /C^ -1 } in the interval of constancy [a%-i, a;] let us 
call this interval of the first kind. Otherwise, if z(t, b) S {±S<5 mi „ | (to, n) S 
<SV-i} on [ctj-i, «i] we call the interval of the second kind. 

Now we note that for u w (-, b) we can not find a bound for H^H^m T ff ) 
independent of the parameter w because the projection q w (-, b) in an interval 
of constancy / = [a*_i, cti] of the second kind, say v(t, b) = 5 m ,n on /, reads 



q w (t, b) = q°°(t, 6) + C,n(*. b) 



and, 



> 



^riL 2 (o,T,i 



'(t)\ 2 dt = w 2 



cos 2 (wt) dt 



~w 



I(i _ 2 -) + J- (sin(2«;(a i - -)) - sin(2u;(a l _i + -))) 
2 w 4u> V w w / 



where L is the length of /. Since w > 3 (and L > 0) we obtain 



^rili 2 (o,r, 



> w 1 



L 1 



— + — sm(2w(oti )) - sin(2w(ai_i H )) 

6 4w v 



Then we see that when w goes to oo also ||0™||l 2 (o,t,r) does and then, so does 
\\uT\\li(o,t,h)- 
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If we consider the projection U w (-, b) of u w (-,b) onto (R ftJV - 1 )^, there holds 
that we can find a bound for ^)||l 2 (o,t, (r k «-i)^) independent of both 

parameters b and w. If fact multiplying by U™(-,b) the equation 



Uf{-, b) = -vAU w {-,b) - p-^^P" 7 Bu w {-, b) + P- K > 
that is satisfied by U w we obtain 
v d 



\uT{tM\ 2 <-^ Jt \\u w (tM\ 2 



\F\\u?{t,b)\ 

+ C\\U w ^b)+ q w (t,b)\\\U w (t,b) + q w ^b)\ m \ur^b)\. (4.11.5) 
After using appropriate Young Inequalities and integrating 

' T \U?(t,b)\ 2 dt 
<Ci\\U w {-,b)\\ 2 c{[0/r]Av ^ N _ 1) ± ) + Ci\F\ 2 

+C 1 \\U w (-,b) + q w (-Mho,T],V)\\U w (-,b) + q w (-MlHO,T,D(A)) 
<C\\U w (;b)\\ ci[ ^ Tl{RKN) ^+T\F\ 2 



-C 



\^ W ('^)\\c([Q : T], (R K «)^) + lk W (")^)llc([0,T],K K W-i) 



^JV-l ) 

(4.11.6) 



Multiplying the equation by U w (-,b) and by AU w (-,b) we obtain that U w (-,b) 
satisfy the estimates (j4.10.14j) . (14.10. 15|i . Ij4.10.17jl and j4.1().18j with U J > L re- 
placed by U w (-, b) and we can easily see that there is a constant C w ^ depending 
only on \\q w {-, fr)||c([o,T],R K «-i) such that 

II U W (-,b) || c( [ 0i T ] j (K~jv)i) < C Wi b & ||£7™(-, ^)||l 2 (o,t : (k"™- 1 )^^)) ^Cw,b- 

(4.11.7) 

Moreover the family { [//"(•, 6)} is uniformly bounded, w.r.t. 6 and w, in the norm 
of L 2 (0, T, H) because, in the "a priori"-like estimates above and in (j4.11.6jl we 
can take bounds independent of the parameters depending only in the bound 



c (6, w)e.Bx[3,+oo[{||9 0) ^)llc([0 : T],R K »-i)} 



< 



Il9°°('i fr)llc([0,T],R K «-i) + 



(m, n)GSjv-i 



C([0,T],R K ™-i) 



< c 



where C does not depend neither on b nor on w. Indeed by (j4.1().36jl the 
family v(-, b) is uniformly bounded in L 2 (0, T, R^-^-norm which implies that 
q°°(-, b) is uniformly bounded in C([0, T), R ftJV - 1 )-norm. Indeed replacing u 
by u°°(-, b) in the equations Ij4.11.3jl we see that, since the family of controls 
{v(-, b)} is uniformly bounded in £ 2 (0, T, H ) norm, we can find a uniform bound 
D for the solutions u°°(-, b) in anyone of the norms C([0, T], H), L 2 (0, T, F), 
C([0, T], V) and L 2 (0, T, D(^4)). In particular 
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&)||c([o,ri,R-~-i) < D^i:, b)\\ C ([o,T],H) < D X D. 
On the other hand || E(m,™) G s JV _ 1 0m,n('' fe )llc([o,T],R"™-i)} never exceeds the 

value 2\/2S. The we can the write: 

Corollary 4.11.3. There is a constant C > suc/i f/iat /or all (w, b) in the 

product [3, +oo[x£>: 

Wti'i & )IIl 2 (0,T, (R KN - 1 )j 1 ) ^ 

Another Lemma we shall use is the following: 

Lemma 4.11.4. Let {z(-, a) e VT 1 ' 2 ([t i , t/], M) | a e £} 6e a uniformly 
bounded, w.r.t. a, family, i.e., 

3C > OVct g S IK, o-)||^i,2 ([t . it/]iK) < C. 

T/ien t/iere is a constant D\ depending only in C and in the length tf — ti of the 
interval [U, tf] and so, independent of a and w, such that 

|| sva(wt)z(t, o")|| rx < Diw^ 1 , and \\ cos(wt)z(t, a)\\ rx < Diw^ 1 . 

Proof. The Lemma follows by direct computation: Let s, r belong to [ti, tf]: 

/sin(wt)z(t, a) dt 

= — w^ 1 / — wsm(wt)z(t, a) dt — w^ 1 / (— cos(wt))z(t, a) dt 

J s J s ML 

= — w^ 1 (^\cos(wt)z(t, cr)Y s — J cos(wt)(-^z(t, a)) dt 
< — w^ 1 [cos(wr)z(r, a) — cos(ws)z(s, a)] 

cos2 ( w *)*) [J ll^(*. ^Il 2 ^ 

<2 W - 1 ||z(t, a)\\c( [u ,t f ],m+ w ^\ r ~ s ^ C - 

By the continuity of the embedding W 1 - 2 — > C° there is a constant C\ such that 
IMIc < Ci ||zi|| v^ 1 - 2 ■ Putting C2 = max{C, C\C} we obtain 

f s sm(wt)z(t, a)dt < w~ 1 C 2 (2 + (t f - U)*). 

Analogously we arrive to 



/' 

J s 



cos(wt)z(t, a) dt < w~ 1 C 2 (2 + (t f ~U)i). 

J s 

Hence 

Vo- e SVs, r e [tj, t/] 
Therefore 

Vo- e s 

Choose L»i = C 2 (2 + (t/ -ti)'). □ 



J s r sin(wt)z(t, cr)o!t < w" 1 C 2 (2 + (t/ - 
J s r cos(wt)z(t, cr)dt < W - 1 C 2 (2 + (t f - 



'||sin(u;t)z(t, o-)|| rx < vj- 1 C 2 (2 + (t f - U)*) 
|| cos(wt)z(t, a)\\ rx < w- 1 C 2 (2 + {t f - U)i). 
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Corollary 4.11.5. Let {z(; a) G W 1,2 ([0, T], R) | a € £} fee a uniformly 
bounded, w.r.t. a, family: 

3C > OVa G £ ||z(-, <x)|| W i, 2([0 , r] ,a) < C. 

T/ien there is a constant D 2 depending only in C and so, independent of a and 
of the parameter b of our controls, such that 

\\4> w (-, b)z(-, o-)\\ rx < D 2 w~ 1 . 

Proof. Let s, r be in [0, T] and suppose without loss of generality that s < r. 
Put D :={te [0, T] I 4> w (t, b) ^ sin(u;<)}. 



<fi w (t, b)z(t, a) dt = / 4> w (t, b)z(t, a) dt + / sxn(wt)z(t, a) dt 

JDn[s,r] J[s,r]\D 

< 2-dC + mw- 1 C 2 (2 + Ti) 

w ' 

< w- 1 C 2 (2T + 2m + T^m). 24 

Choose D 2 = C 2 (2T + 2m + T^m) . □ 

Nowwe compare the projections U w (-, b) and U°°(-, b) of respectively u w (-, b) 
and u°°(-, b) onto (M Kjv - 1 ) i ^. We claim that at time ctj, (i = 1, . . . , m) there 
holds 

llHa*, 6) - 6)| < aw- 1 (4.11.8) 

where Ci is a constant independent of the parameters w and b. In particular at 
time T there holds 

|Cr"(T, 6) — U°°(T, 6)| < C^w -1 . (4.11.9) 

Note that H4.11.9|) implies Lemma 14.10.91 because as we have seen at time T 
we have q w (T, b) = q°°(T, b). To prove (|4.11.9|) we shall compare U w {-,b) and 
U°°(-, b) in each interval of constancy: 

• In an interval of the first kind U w (-, b) and U°°(-, b) satisfy the same 
equation 

U t °°(., 6) = -vAU°°(; b) - p-^-ipV^oo^ ^ + ? oo^ b y^ + p-K N - ip . 
U?(; b) = -vAU w {-, b) - P- KN - 1 P V (U w (-, b)+q°°(; b)) + p- K «-i F 

• In an interval of the second kind, say v(-, b) = ±S<5 m .„ (m,n) G Sjv-ij 
[/"'(•, 6) and U°°(-, b) satisfy the equations 

C/ t °°(-, 6) = -zME7°°(-, 6) - ?- K «-^ v ([/»(. i 6) + g °°(. 7 6)) + p-"*-iF; 
f/f (-, fe) = -zM[r(-, 6) - p-««-ip v rE/*^, 6) + 6) 

+ v / 2S<T(', 6)(e m ±e n ) 



p-«w-ip ±s ^ 



24 Note that m being the number of intervals of constancy, then [s, r] \ D is a union of at 
most m intervals. 
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Define the difference 

V W (;b) :=U W (;b)-U°°(;b). 

For an interval on] of second kind we find 



T}?( S , b)=- V A V W (;b)+P 



-KN-l pV 



U°°(; b)+q°°(; b) 



-KJV-ip>V 



U w (; b) + q°°(-, b) + VZZcf> w (; b)(e m ± e„)) 



Multiplying by ?7 tu (-, 6) we obtain 
1 d 



2dt 

--v\\v w {s, b)f 



\V?(S, b)\ 2 

[p- kn - 1 p v (u°°(s, b) + q°°(s, b)), V w (s, b) 
P- KN ->P v {U w (s, b) + q°°(s, b) + y/mp°(8, b)(e m ±e n )), V w (s, b) 
± (sp- K --^ m;n , (s, 6)) 



so, 



< 



~\V W (S, b)f + v\\ V w (s, 6)f 

-kn-1 pV (tjw I 



± 



(4.11.10) 



Jp- w - ip v( tr(j) b)+q°°(s, b)), V w (s, b)) 
(P" KN - 1 P v (?7 00 (s, 6)+<z°°( S , 6)), ^(s, 6)) 

(P- KN - 1 P V (^(S, 6)+<z°°( S , 6), V2S^( S , 6)(e m ±e n )), rf (s, 6)) 
(p- kn - 1 P v (v / 2S^(s, 6)(e m ±e„), !7 w (s, + 6)), 6)) 

(P" KN - 1 P v (\/2S^(s, 6)(e m ±e n )), V w (s, 6)) 
( S p-^-^ roi „, ^(s, 6)) 
=A 1 (s) + A 2 {s) + A 3 (s) + A 4 (s) 
where we define 

A^s) := - (P- K "- 1 P V (C/^( S , b) + q°°( S , 6), v w (s, 6)) 

+ (p- KAr - 1 P v (C/ 00 (s, 6)), r/ w (s, &)) 

4>(s) := 

-(p- KAr - 1 P V (^(.s, &) + <?°°(s, b), V2S<T( S , 6)(e m ±e„)), 6)) 

Ms) ■= 

-(p-^-'P v (V2Sr(s, b)(e m ±e n ), U w (s, b) + q°°( s , b)), V w (s, 6)) 
M*) ■= - (p- Kn - 1 P V {V^4> w {s, b)(e m ± e„)), V w (s, 6)) 
± (sp- K «-^ m ,„, V w (s, 6)). 
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To estimate \rj w (z, b)\ 2 z S [a-i-i, oti\ we will integrate (|4. 1 1 . 1()|> . In the second 
member, for we have: 

\A X \ < C\ V w (z, b)\\\ V w (z, b)\\\\U™(z, b) + q°°(z, (4.11.11) 

and by Young Inequality we arrive to 

f t \rj w (z, 6)| 2 + HI^(^ b)\\ 2 

4 

<D\rj w (z, b)\ 2 \\U°°(z, b) + q°°(z, b)\\ 2 +2^1 r (z) 

and, by Gronwall Inequality we obtain 

\rf[z, b)\ 2 <\rf (oi-i, 6)| 2 exp(l) f°" \\U°°(s, b)+q°°(s, b)\\ 2 ds) 

4 ^ 

+2 Y1 / Ar(s)E(s)ds (4.11.12) 
with E(s) := exp -D\\u°°(t, b)\\ 2 dt. For r = 2 we have 
A 2 (s)E{s)ds = 

2S / 0*( a> b)(p- K »-ip v (U w (s, b)+q°°(s, b), em), V w (*, b))E(s)ds 

<t> w (s, 6)(P- K "- 1 P V (^( S , 6) +g°°( S , 6), e„), (s, bj)E(s)ds. 

(4.11.13) 

Now we estimate the derivative of the product 

(p- re -- 1 P v (C/ w (s, 6) + <T( s , b), e m ), rj w (s, by)E(s) : 



T v 2S 



(p-^- 1 P v (C/ w (s, 6)+<T( s , 6), e m ), 77 u '(s, &)) 
= (p^ k "- 1 P V (^(s, 6) + q°°(s, 6), e m ), r, w (s, 6)) 
+ (P- KN - 1 P V {U W ( S , b)+q°°(s, b), e m ), f) w (s, 6)) 

<C\U w {s, b)+q°°(s, b)\\e m \ [2] \\ V w ( Sl 6)|| 
+C\\U w (s, b)+q°°(s, b)\\\e m \ [2] \f, w (s, 6)| 

<D\U w (s, b) + f°{e, b)\\\rf(8, b)\\ + D\\U w (s, b) + q°°(s, b)\\\f, w (s, b)\ 



so. 



I — — ( P _ Kj 



-*P v (U w (s, 6)+g°°( S , 6), e m ), 7/ w (s, 6))| 2 rfs 
<A / |f/ w ( S , &)+g°°(*, 6)| 2 ||rr(s, b)\\ 2 ds 

J 0>i-l 

+D, \\U w ( S ,b)+q™( S ,b)\\ 2 \r(s,b)\ 2 d S 



<Xi-l 
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and, for 



Jp-w-tpV^^ b)+q°°(s, b), e m ), (s, 6)) 



we obtain the bound 



(p-*»-ipV(ir(s, b) + q°°(s, b), e m ), if( S , &)) 
<C\\U w (s, b) + <z°%s, 6)||||e m ||||^( Sl fe)|| < D||[/-( S , b)+f°(a, b)\\\\r) w (s, b)\\. 



Hence, by (|4.11.4|l and (|4.11.7|l and Corollaries 14.11.21 14TT3I we conclude that 
the family 

{(P- K »-1P V {U W ( S , + 6), e m ), ^(s, 6))} 

is uniformly bounded, w.r.t. w and 6, in W ' ([ati—i, ai])-norm. 
For E'(s) we find 

\E(a)\ < exp(T£>|| u °°(., b)\\ c{[0 , n v)) < C x (4.11.14) 
|^JS7(a)| < \\u°°(; 6)110(10,7!, viCi < C 2 (4.11.15) 

where C\ and C2 do not depend neither on b nor on s. In particular E is 
uniformly bounded in W 1,00 ([cti-i, oti], R)-norm and then also the family 

{(P- KN -1P V {U W ( S , b)+q°°( s , b), e m ), V w (s, b))E(s)] 

is uniformly bounded in W ' ([aj_i, aj], R)-norm. 

Therefore by (|4.11.13|) and Corollarv l4. 1 1 .51 there is a constant K\, indepen- 
dent of z and of the parameters w and b, such that 



<T( S , b)(p- K ^P v (U w (s, 6) 6), e m ), V w (s, b))E(s)ds 

1 



If we replace m by n we obtain a similar estimate and then, we conclude that 
for some constant K2 independent of z, w and b: 



A 2 (s)E(s) ds 



< K2VJ- 1 . 



(4.11.16) 



Analogously we conclude that for some constant K% independent of z, w and b: 



A 3 (s)E(s) ds 



(4.11.17) 



Recalling equation (|3.3.8|1 and supposing without loss of generality that m < n, 
we have 

P v B[VZE<b w (s, b)(e m ± e n )] = 2H(^°(«, b)) 2 P v B[(e m ± e„)] 

= ±2E(0™(«, b)) 2 6 m , n . 
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Therefore 



Ai{s)E(s) ds 



f [T2S(<T( S , 6)) 2 ±S](p- K «-^ m , n , r, w (s, bj)E(s)ds 

J aii— i 

f [T2H(^(s, 6)) 2 ±S](^„, ^( s , b))E(s)ds 
[2{<t> w {s, b)) 2 - 1] (8 m , n , rf(s, b))E(s)ds 



S y cos(2«;f)(<5 m , n , ^(s, bfjE(a)da 



|(<5 m , n , ^(s, &)||£(«)|ds. 



whe re I = [a j-i, z] n [aj_i + ^f, a* - ^f] and i, ; = cm - 
By (I4.11.14|) we have that 

/ \(6 m , n ,ri w (s,b)\\E(s)\ds 

J[o 4 _i,a]\J 

2L 

<Cl|(<W»> ^(s, &))|c([0,T],»)— 1 < ^"^oll^Cs, &))llc([0,T],(K»«-i)^) 

where the constant i^o can be taken independent of (m, n) e Sjv-Ii because 
iSjv-i is finite. 

For the derivative of the product (S m<n , rj w (s, b)) we find 



d 
dt 



(<W, ^ = (*m,n, V W (s, b)) < Ktl^is, 6)| 



so, by H4.11.4J1 . H4.ll.7jl and corollaries 14 . 1 1 . 21 and 14 . 1 1 . 31 we conclude that for 
some constant C2 independent of w and b: ||(5 m , n , rj w (-, 6))||w 12 ([o,t],r) < 
C 2 . So by the uniform boudedness of E{-) in W 1 ' 00 '^, T], K) we conclude the 
uniform boudedness of (5 m ,„, ^(s, in W x ' 2 ([0, T], M). 

By Lemma T4. 1 1 .41 there is a constant C\, independent of w and b, such that 



' cos(2tut)($ m , n , r] w {s, b)^E(s)ds 



< Ci(2w)- 1 . 



Hence 



A 4 (s)E(s) ds 



C 

< -yw^ 1 +w- 1 K \\r] w {s, &)|lc([o,T],(R' t «-i )j L ) 



and, again by Ij4.ll.4jl and Ij4.ll.7jl we arrive to 



Ai(s)E(s) ds 



< K 2 w~ 



(4.11.18) 



where Ki is independent of b and w. 
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By equation Ij4.11.12j> and estimates (|4.11. lfjji . (j4.ll.17j) and Ij4.11.18jl we 

obtain 

\v w (z, b)\ 2 

^rfia^x, b)\ 2 exp(D / \\u°°(s, b)\\ 2 ds) + 2 V / A r (s)£(s) <Zs 

In an interval of the first kind we just have A 2 = A3 = A4 = on (j4.11.1L)|) and 
then we obtain 

W{*, b)\ 2 

<\v w (ai-i, 6)| 2 exp(£»^ ||t/°°( s , h)+ 9 00 (s, fc)|| 2 rfs) < (a^ 1: b)\ 2 . 

Therefore in any interval of constancy (either of first or second kind) we have 

ItH*. b )\ 2 < D 1 \r, w {a i ^ 1 , b)\ 2 + D 2 w~ x (4.11.19) 

where D\ and D 2 are independent of z € [oti—i, oti], of the parameters w and b 
and of i € {1, 2, ... , m}, i.e., of the interval of constancy. 
We prove Ij4.11.8j) by induct ion on 1: 

• i = 1: At time en, by Ij4.11.19jl . we have 

|Cr(ai, b)-U°°(a u b)\ 2 KDaw- 1 25 
Choose C\ = D 2 . 

• Induction Step: Suppose that at time we have 

I^Ca^i, 6) - U°°(a^ u b)\ 2 < C~ lW -\ 
Then by I4.11.19jl . we have 

\U w {a i . l , b)-U x {a^ x , b)\ 2 

< D 1 {C~~ 1 w- 1 ) 2 + D2W- 1 = [DiCVT^- 1 +Z) 2 ] W - 1 

< [D 1 C~^ 2 ^+D 2 ]w-\ 

Choose C~i = D x Ci- 2 \ + D 2 . 
Therefore (j4.11.8jl holds. 



Remember that both u°° and u w start at uq at time 0. 



Chapter 5 

L 2 - Approximate 
Controllability. 



The following Proposition says that for any T > 0, system (|4.9.3|) is time-T 
approximately controllable in L 2 -norm. 

Proposition 5.0.6. For any uq € V and T > 0, the attainable set at time T 
from uq of system (|4.9.3(l is dense in H . 

Proof. Fix e > 0, uq £ V, x\ G H and T > 0. We prove that it is possible to 
drive the system from u to the ball {y g H \ \y — xi\ < e}. For that, first 
we set M e No such that \xi ~ P KM x\\ < |, where P KM is the projection map 
onto R KM and, consider the system [ (|4.10.1I) .M], As we have seen in the first 
step of the proof of Proposition 14.10.21 (section 14.10.1(1 , there is a control % 
taking values on M. KM and driving the system from uo to some point u^F such 
that P Km Ut = P KM xi. Moreover by item[T]of Corollary ETTTPI we have 



m _ p «m u M| < ^[ T exp(T)]i 



Now (see Lemma l4.10.9f) we imitate vm by another control vm-i taking values 
in I" 1 '- 1 and driving the system to a point w^f -1 £ V such that 



2M 

Repeating this procedure of imitation, at each step we find a control Vi~i taking 
values in R^- 1 and driving the system to a point u l rf l S V such that 



u l T 1 — u l T | < 



e 

2M' 



Let T <T\ where T\ equals the unique solution of T\ exp(Ti) = ( 2 £ M ■ The 

control v 1 takes its values on IR Kl and drives the system to a point u\, G V 
satisfying 



M 

l«T - *l| < I^T 1 ~ U t\ + \ U T ~ P KM u¥\ + \P KM u¥ - Xl \ 
i=2 
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Since P KM u^ = P KM x 1 we have 

14 ~xi\< ^-(M - 1) + K[Tcxp(T)]i + £ -<e. 

Hence the Proposition is proved for T < T\. If T > T\ we apply an arbitrary 
control up to time T — T\ arriving to some point yo G V and then drive the 
system in time Ti from y to the ball {y G H \ \y — x\ \ < e}. □ 
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